Exercises for Book 1 


Book 1 Chapter 1 


Exercise 1.1. What fraction of their electrons would 
a person need to lose to have a net charge of 1 wC? 
For simplicity, imagine that the person has a mass of 
60 kg and is composed entirely of water molecules. 
Each water molecule has a mass of 3.0 x 102° kg 
and contains 10 electrons. 


Exercise 1.2 Suppose that one gram of protons is 
separated from another gram of protons by 

3.84 x 10° km (the distance between the Earth and the 
Moon). Estimate the magnitude of the electrostatic 
force between these two concentrations of charge. 


Exercise 1.3. Two identical point charges, Q, are a 
distance d apart. A point charge q is placed on the line 
joining these charges, midway between them. What 
value must be chosen for g to ensure that all three 
charges experience zero electrostatic force? Is the 
system then in stable equilibrium? 


Exercise 1.4 Four identical point charges, gq, are at 
the corners of a square with length of side d. (a) What 
is the magnitude of the electric field at the centre of 
the square? (b) What is the magnitude of the electric 
field at the centre of the square when one of the corner 
charges is removed? 


Copyright © 2010 The Open University 
2.1 PT/TRx 


eS Science Supplementary Material SMT359 AddEx 
( ) Da SMT359 Electromagnetism 
— 
WO 
= 
4 = 
=< SMT359 
fo) 
QO. 
oO 
Ww e,e e 
< Additional Exercises 
Contents Exercise 1.5 Four charges are arranged in a square 
with side of length d (see Figure 1). The (x, y) 
Exercises for Book 1 1 coordinates of the charges are as follows: 4q at (0, d), 
Exercises for Book 2 4 2q at (d, d), 2q at (0,0) and —q at (d,0). What is the 
. electric field at the centre of the square? 
Exercises for Book 3 11 
Solutions for Book 1 16 Yh 
Solutions for Book 2 25 a 
Solutions for Book 3 44 


+4¢@------ = @aege* 
O 


| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
+2 -¢ = 


Y¥__» 
< Fi > 
Figure | For Exercise 1.5. 


Exercise 1.6 What is the electric field at the centre 
of a uniformly charged disk of radius R and charge Q? 


Exercise 1.7 A positive charge @ is uniformly 
distributed over a thin circular ring of radius R. (a) 
What is the electric field at the centre of the ring? (b) 
What is the electric field at the centre of the ring if a 
short segment of charge of length 1 < R is removed 
from the ring? 


Exercise 1.8 An electric field vector E has the 
value (3e, + 4e,) NC~!. What is the magnitude 
of this field and what is its unit vector? Find the 
component of the field in the direction of the unit 
vector @ = (e; — ey)/V/2. 
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Book 1 Chapter 2 


Exercise 2.1 A charge distribution is cylindrically 
symmetrical around the z-axis and is described by the 
charge density 


A 
pA(r)=— forO<z<Land0<r<R, 
iia 


where r is the distance from the z-axis and A is a 
constant. This charge distribution is independent of z 
for 0 < z < L, but vanishes outside this range and for 
r > R. What is its total charge? 


Exercise 2.2, A hydrogen atom consists of a proton, 
represented by a point charge +e at the origin, and an 
electron, represented by a spherically symmetric 
charge distribution of charge density 


p(r) = Aexp(—2r/ao), 
where r is the distance from the origin, A is a constant 
and ag = 5.29 x 107'' m is the Bohr radius. (a) Use 
the fact that a hydrogen atom is electrically neutral to 
find a formula for A. (b) Hence evaluate the charge 


density of the electron cloud at the Bohr radius, ao. 
You may use the standard integral 


| az” exp(—2x) da = 2. 
0 


Exercise 2.3. Consider the electric field 

E = Aexp(—r/a) e,, where e, is a radial unit vector 
pointing away from the origin and A and a are 
positive constants. Show that such a field could not be 
produced by an isolated point charge at the origin, 
even one moving rapidly enough to rule out any 
arguments based on Coulomb’s law. 


Exercise 2.4 Consider a spherically symmetric 
vector field F = exp(—Ar) e,, where A is a positive 
constant and e, is a radial unit vector pointing away 
from the origin. Does F obey the divergence theorem? 
Is the flux of F equal to zero through a closed surface 
that does not enclose the origin? 


Exercise 2.5 Show that there is no electric field 
inside a uniformly charged infinitely long cylindrical 
shell. 


Exercise 2.6 Consider a cylinder-shaped uniform 
distribution of positive charge. The cylinder has finite 
length and negligible radius. The electric field due to 
this distribution is measured at a point P outside the 
charge distribution and on its central axis. If this 
charge distribution were shrunk down to its central 
point, would the magnitude of the field at P increase 
or decrease? 


Exercise 2.7 The electric flux over a spherical 
surface of radius R, centred on the origin, does not 
depend on any charges that are outside this surface. 
Does this mean that the electric field at any point a 
distance R from the origin does not depend on charges 
that are further from the origin than R? 
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Exercise 2.8 Is it possible to directly use the integral 
version of Gauss’s law to find the electric field at 
every point in space due to a cube of charge with 
length of side Z and a uniform charge density p? 


Exercise 2.9 The membrane of a spherical 

animal cell has a uniform surface charge density 

of —2.5 x 10-°C m“” on its inside surface and 
+2.5 x 10°-°C m~? on its outside surface. The 
thickness of the membrane is much less than the 
radius of the cell. Ignoring any polarization effects 
and assuming there are no other concentrations of 
charge, determine the magnitude of the electric field 
(a) inside the membrane, (b) inside the cell and (c) 
outside the cell. 


Exercise 2.10 An isolated infinite slab of uniform 
charge density p > 0 lies between z = —d/2 and 

z = d/2. (a) Use the integral version of Gauss’s law 
to find the magnitude of the electric field at all points 
inside and outside the slab. (b) Confirm that your 
answer is consistent with the differential version of 
Gauss’s law. 


Exercise 2.11 An infinite plane at z = d/2 has a 
uniform surface charge density o (see Figure 2). It 
is in contact with an infinite slab of charge with 
uniform charge density p extending from z = —d/2 
to z = d/2. All the charges are fixed. Where, if 
anywhere, in the slab does the electric field vanish? 


surface 
charge 
density =o 


RY 


tz = —d/2 


volume 
charge 
density = p 


Figure 2 For Exercise 2.11. 


Exercise 2.12 A hydrogen atom can be modelled as 
a point charge +e at the centre of a spherically 
symmetric cloud of negative charge with charge 
density p = —eexp(—2r/ag)/(mag), where 

ag = 5.29 x 107" m is the Bohr radius. What is the 
radial component of the electric field at ag? 


You may use the fact that fe x? exp(—x) dx = 0.647. 


Exercise 2.13. A hemispherical shell with radius R 
has a uniform areal charge density o (see Figure 3). 
The centre of the hemisphere is at the origin, and the 
shell extends out towards values with z > 0. Consider 
now the imaginary surface S that forms the base of 
the shell, i.e. the circular, flat surface with radius R in 


the xy-plane. Show that, at all points on the surface S, 
the electric field is perpendicular to S. 


Exercise 2.14 (a) A point charge Q is at the centre 
of a cube. What is the electric flux over each face of 
the cube? (b) The charge is now moved far into a 
corner of the cube, whilst remaining wholly within the 
cube. What is the electric flux over each face of the 
cube in this situation? Hint: Embed the original cube 
in a larger cube with double the side length so that the 
point charge is close to the centre of this larger cube. 
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Figure 3 For Exercise 2.13. 


Exercise 2.15 Describe the electric field lines 
outside two concentric spherical conducting shells. 
The inner shell has a positive charge q; > O and the 
outer shell has a negative charge q. < 0. Consider the 
three cases qi + q2 < 0,q1 + q2 = Oand q + gq > 0. 


Exercise 2.16 The detailed study of ions is a topic 
of current interest in physics and chemistry. These 
studies are aided by building ion traps, which confine 
the ions in a trapping region which is a vacuum 
occupied only by trapped ions themselves. Consider a 
proposal for building a trap by distributing charges 
outside the trapping region in such a way that they 
produce a suitable electric field to trap positive ions. 
To achieve stable equilibrium, the electric field 
would have to exert a restoring force on the trapped 
ions, always pushing them back inside the trapping 
region. Use Gauss’s law to show that this proposal is 
unworkable; the required electric field cannot be 
created by charges outside the trapping region. 


Exercise 2.17 An electric field that points in the 
x-direction depends only on x and is equal to zero at 
x = 0. The charge density responsible for this field is 
p(x) = po cos(ka:). What is the electric field? 


Book 1 Chapter 3 


Exercise 3.1. The Earth’s magnetic field is similar to 
that of a magnetic dipole with a magnetic moment of 
8.0 x 102? A m?. If this were due to a loop of current 
with radius 2350 km, what would the magnitude of 
current be? 


Exercise 3.2. The current density in a region of 
space varies as 


J = Aexp(—k-r), 


where A and k are constant non-zero vectors. What 
is the direction of the current flow? Under what 
conditions does this current density lead to no 
accumulation of charge? 


Exercise 3.3. A circular loop of wire carries a steady 
current. Use the Biot—Savart law to show that the 
magnetic force on any segment of the wire is in the 
plane of the loop and points radially outwards from its 
centre. 


Exercise 3.4 Show that the electric field in an empty 
Penning trap is bound to have a radial component that 
opposes the confining effect of the magnetic field. 
(You will need to think back to an earlier chapter to 
answer this question, remembering that the trapping 
region is a vacuum with zero charge.) 


Book 1 Chapter 4 


Exercise 4.1 A current distribution is spherically 
symmetric, flowing radially away from the origin. Use 
symmetry principles to show that the magnetic field of 
this current distribution is radial. Then use the 
no-monopole equation to show that the magnetic field 
actually vanishes. 


Exercise 4.2 A thick circular ring of current is 
centred on the origin and lies in the plane z = 0. Use 
symmetry principles to express the magnetic field near 
this current distribution in the simplest possible terms. 


Exercise 4.3 A long coaxial cable consists of a 
cylinder surrounded by a cylindrical tube, both with 
the same central axis. The cylinder and the tube are 
both conductors, but the space between them is filled 
with an insulator. Steady uniform currents flow in the 
cylinder and the tube. These currents have the same 
magnitude J, but flow in opposite directions. What is 
the magnetic field outside the coaxial cable? 


Exercise 4.4 Two infinite current sheets are 
arranged so that they are parallel, and separated 

by a distance d. Sheet | is in the plane defined 

by z = —d/2, sheet 2 is in the plane defined by 

z = +d/2. Each sheet is made up of wires aligned in 
the x-direction with N// wires per unit length in the 
y-direction. Each wire carries a current of magnitude 
I, but one sheet carries it in the positive x-direction 
while the other carries it in the negative «x-direction. 
What is the magnetic field between the current sheets? 


Book 1 Chapter 5 
Exercise 5.1 Calculate the potential energy of a 


uniformly charged sphere of charge Q and radius R. 


Exercise 5.2. What is the potential at the centre of a 
uniformly charged sphere of charge Q and radius R? 


Exercise 5.3. What is the potential at the centre of a 
uniformly charged disk of charge Q and radius R? 


Exercise 5.4 In some respects an atomic nucleus 
behaves like a uniform spherical distribution of 
charge. A spherical nucleus of charge Q and radius R 
divides into two isolated spherical nuclei, each of 
charge Q/2 and radius R/2'/3. By how much does 
the electrostatic potential energy of the system 
change? 


Exercise 5.5 In a particular lightning flash, the 
potential difference between the cloud and the ground 
is 1.0 x 10° V and the quantity of charge transferred is 
30 C. What is the change in energy of the transferred 
charge? 


Exercise 5.6 Show that, in equilibrium, the electric 
field is perpendicular to the surface of any conductor. 


Book 1 Chapter 6 


Exercise 6.1 Consider a rectangular conducting 
loop ABCD in the plane z = 0 moving with constant 
velocity v = v;e,z in the x-direction (see Figure 4). 
The points ABCD defining the loop initially have the 
following (x,y) coordinates: A(0, £1), B(0,0), 
C(L2, 0), D( Le, £1). Hence sides BA and CD are 
parallel to the y-axis and have length L,. Sides AD 
and BC are parallel to the z-axis and have length 
Ly. A magnetic field B = B,(x) e, points in the 
z-direction and is a function of x, but is independent 
of y, z and t. Show that the induced emf in ABCD is 


Vemt = Uy, L1|Bz(2) | BAe Te Ly)| 


where z is the instantaneous x-coordinate of side AB 
and x + JL is the instantaneous x-coordinate of side 
DC. 


y 
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Figure 4 For Exercise 6.1. 


Exercise 6.2 A rigid loop maintains a fixed 
orientation whilst moving through a uniform static 
magnetic field. Is any emf induced in the loop? 
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Exercise 6.3 A constant uniform magnetic field of 
magnitude B points in the positive z-direction. A 
conducting rod of length L is pivoted at one end and 
rotates with angular speed w anticlockwise (as seen 
from an observer at z > 0) in the plane z = 0. Use the 
magnetic force law to calculate the voltage drop from 
the pivoted to the free end of the rod. 


Exercise 6.4 Use Faraday’s law to provide an 
alternative solution to Exercise 6.3. 


Book 1 Chapter 7 


Exercise 7.1 Inside a conducting medium that obeys 
Ohm’s law, the current density at a given point is 
related to the electric field at the same point by the 
equation J = oF, where a is the conductivity of the 
medium. Use this fact, together with the equation of 
continuity to show that a concentration of excess 
charge within the medium will decay exponentially in 
time. 


Exercise 7.2 Verify explicitly that the fields 
1 
E=f(z—cd)e, and B=-—/f(z—cd)e,, 
: , 


where c = 1/,/€pji0, obey all four of Maxwell’s 
equations in empty space. 

Exercise 7.3 Figure 7.10 in Book 1 (p. 182) shows a 
linearly-polarized monochromatic electromagnetic 
wave that propagates along the z-axis. Rank the order 
of the electric field strengths at points P, Q and R, 
which are all in a plane perpendicular to the direction 
of propagation and have spatial coordinates P(0, 0,0), 
Q(A, 0,0), R(A, A, 0), where 2 is the wavelength. 


Exercise 7.4 (a) How would you orient a long 
straight aerial to detect a plane, linearly-polarized 
electromagnetic wave? 

(b) How would you orient a circular-loop aerial to 
detect a plane, linearly-polarized electromagnetic 
wave of wavelength similar to the diameter of the 
loop? 


Exercises for Book 2 


Book 2 Chapter 1 


Exercise 1.1 Charge is distributed throughout a 
sphere of radius R with a charge density given by 
kr, r<R, 
p(r) = 
0, r>R, 
where r is the distance from the centre of the sphere 
and k is a constant. 
(a) What is the symmetry of the problem? 


(b) Use the integral form of Gauss’s law to find the 
electric field E at all points in space, both inside and 
outside the sphere. 


(c) Use the differential form of Gauss’s law to check 
your answer to part (b). 


Exercise 1.2 Consider the two vector fields 
F=se,+ (22—y)e, and 
G = —2ze, + (4y + 5)ey. 

Can either of these vector fields represent a static 


magnetic field? If appropriate, find the associated 
current density. 


Exercise 1.3. Given the electric field 
E(r,t) = Eosin [w(t - “| er, 
Cc 


where Eo, w and c are positive constants, use 
Faraday’s law to find the most general form of the 
associated magnetic field B. 


Exercise 1.4 (The derivation of a wave equation in 
this question can be done using what you have 
learned in Book 1. However, this will be discussed in 
detail in Book 3, and you may prefer to try this 
question after studying Chapter I of Book 3.) 


Use the vector identity 
curl (curl F) = grad (div F) — V°F 


to show that, in a region of empty space that contains 
no charges or currents, two of Maxwell’s equations 
can be combined to give a wave equation of the form 


1 @E 
VE=— —~. 
ce Ot 
Your answer should establish a relationship between 


the constant c and two constants appearing in 
Maxwell’s equations. 


Exercise 1.5 An electron of mass m and charge —e 
is moving in a uniform electric field E = He, anda 
uniform magnetic field B = Be,, where F and B are 
constants. Find the equations of motion for the 
electron. 


Book 2 Chapter 2 


Exercise 2.1. A sphere of radius 1.50 mm consists of 
liquid argon, atomic number Z = 18 and number of 
atoms per unit volume n = 2.17 x 107° m~?. When 
placed in a uniform electric field, it is found that the 
total dipole moment II of the sphere has a magnitude 
of 5.50 x 107-4 Cm. 


(a) Calculate the average dipole moment of an 
argon atom in this particular electric field, and the 
corresponding average separation of the centre of 
charge of the electron cloud and the argon nucleus. 


(b) Find the polarization P of the sphere. 


Exercise 2.2 A dielectric sphere of radius R has 
polarization P = ar cos? @e, — arcos@ sin 6 ep, 
where a is a constant and r, 6, ¢ are spherical 
coordinates whose origin is the centre of the sphere. 


(a) Find the bound volume charge density pp within 
the sphere. 


(b) What is the bound surface charge density oj, on the 
surface of the sphere? 


Exercise 2.3 A slab of dielectric material, with 

uniform relative permittivity ¢ = 6.0, is placed in 
air perpendicular to a uniform electric field E of 

magnitude 12000 V m“?. 


Inside the dielectric, what are the magnitudes of 
(a) the electric field, 

(b) the electric displacement, 

(c) the polarization? 

You should assume that ¢,;, = 1. 


Exercise 2.4 A parallel plate capacitor has plates 
with area A and separation d; + dz. Two slabs of 
dielectric materials of relative permittivities ¢; and €9 
and thicknesses d, and dp fill the gap between the 
plates. The free charge per unit area on plate | is 

o. There is no accumulation of free charge at the 
interface between the dielectrics. 


(a) Neglecting any edge effects, what is the free 
charge per unit area on the lower plate? Calculate the 
magnitudes of the electric displacement D, the 
electric field E and the polarization P in the dielectric 
materials between the plates. 


(b) Find expressions for the potential difference across 
each dielectric and the capacitance of the capacitor. 


Book 2 Chapter 3 


Exercise 3.1 A long cylindrical wire of radius R 
carries a steady current with uniform current density 
of magnitude J. The wire is made of a material with 
relative permeability ju. Find the magnetic intensity H 
and the magnetic field B at a point inside the wire a 
distance r from the wire’s central axis. 


Exercise 3.2 State the boundary conditions for the 
B and H fields at the plane boundary between two 
media in the absence of free currents at the interface. 


A uniform magnetic field B crosses a plane boundary 
between two LIH non-conducting media. In the first 
medium, of relative permeability 4; = 1.5, the B field 
is at an angle of 60° to the normal to the boundary. In 
the second medium, of relative permeability juz = 3.0, 
the B field is at an angle of @ to the normal to the 
boundary. Calculate the value of 6. 


Exercise 3.3 A coil has 2000 turns wound on a 
toroidal ring of soft iron alloy of relative permeability 
400. The ring has a mean radius of 12 cm (measured 
to the centre of the iron). 


(a) Assuming that the flux is restricted to the iron and 
that B has constant magnitude within the iron, 
calculate the magnitude of H and B within the iron 
when a current of 1.5 A flows in the coil. 


(b) Suppose now that a break is made in the iron so 
that there is a uniform air-gap of width 3.0mm. 
Calculate the magnitudes of H and B in this case, 
both in the iron and in the gap. 


Exercise 3.4 A sample of a linear, isotropic, 
homogeneous material is placed in a static magnetic 
field B(x). Show that the bound current density J, (r) 
is proportional to the free current density J;(r). 
Hence show that when no free currents flow through 
the sample, the bound current density vanishes, 
ensuring that all magnetization currents are confined 
to the surface. 


Book 2 Chapter 4 


Exercise 4.1 Identify each of the following 
equations, and state the circumstances in which you 
would apply each of them to solve an electrostatic 
field problem. 


(a) divE-* 
E0 
(b) divD = ps 


1 
E: ds = — 
) JES {ed 


(d) [,D-dS = f. ppdr 


Se Vre— 
EO 
() V2v=—fh 
EEQ 
i) V7V =0 


Exercise 4.2. Consider a pair of parallel conducting 
plates, 1 and 2, with plate 1 in the plane x = O and 
held at a fixed potential of zero volts, and plate 2 in 
the plane x = s at a fixed potential V;. Between the 
plates is a dielectric material with relative permittivity 
€ in which there is fixed a fixed charge density 

p(x) = A+ Ba, where A and B are constants. The 
plates are so large compared to their separation s that 
the electric field between them can be considered to 
depend only on z, the perpendicular distance from 
plate 1 towards plate 2. 


Show that the potential in the dielectric material is 


x 1 


1 Afn2 1 p(3_ 2 
A(x B 
a? (2° —sz)+¢B(a°—s°2)|), 


and determine the electrostatic field E in the material. 


Exercise 4.3 A solid metal sphere of radius a is 
concentric with a hollow metal sphere of inner radius 
2a. The space between the spheres is filled with a 
dielectric material of relative permittivity ¢ carrying a 
uniform static charge density pp. The charge density 
is zero Outside the spheres. 


(a) Write down Poisson’s equation in spherical 
coordinates for the electrostatic potential V in the 
region a < r < 2a between the spheres. 


(b) Find the general solution of Poisson’s equation in 
the region between the spheres. 


(c) Now assume that the two metal spheres are held at 
zero potential. Use the boundary conditions, together 
with your answer to part (b), to show that 


vo)~ BE fivo(2)'-7(2)] 


(a= Ps 20): 


Exercise 4.4 Two equal point charges q are located 
at points with Cartesian coordinates (d/2, d/2, 0) and 
(d/2, —d/2, 0) in front of an infinite conducting plate 
at zero potential in the plane 7 = 0. Use the method 
of images to find the electric force acting on the 
charge at (d/2, d/2, 0), and hence show that the 
magnitude of this force is 3q?/87eod?. 


Exercise 4.5 (a) Show that the potential 
Bcocos@ 


V(r, ¢, z) = Arcos¢+ 


is a solution of Laplace’s equation in cylindrical 
coordinates. 


(b) An infinitely long dielectric cylinder, with radius 
FR and relative permittivity ¢, is placed with its axis 
perpendicular to a uniform electric field Eg. Find the 
electrostatic potential inside and outside the cylinder. 


Exercise 4.6 Obtain an equation, like Equation 4.26 
on p. 88 of Book 2, for Laplace’s equation in three 
dimensions. 


Exercise 4.7 Consider the numerical solution of 
Laplace’s equation for a square region with V = 0 on 
the top and bottom sides and V = 1 (volt) on the left 
and right sides. Use the Jacobi relaxation numerical 
method to determine the values of the potential on a 
nine point square grid within the square region. The 
spacing of the grid points is i of the side length of the 
square region, and the grid is centrally located within 
the region (see Figure 5). Use the starting values for 
the potential at the internal grid points as shown in the 
figure, and work to three significant figures. Note that 
the bold values, O and 1, are the boundary potentials, 
and are therefore fixed. 


1 aw 5 405 WW 5 1 

1 40 5 40.5 40 5 1 

1 40 5 40.5 40 5 1 
0 0 0 


Figure 5 For Exercise 4.7. 
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Exercise 5.1 An infinitesimal current element J 6!’ 
is located at the point (0.02, 0, 0) m in a Cartesian 
coordinate system. It points in the direction e, 

and has magnitude 5.0 x 10-° Am. What is the 
magnetic field due to this current element at the point 
(0.08, 0.08, 0) m? 


Exercise 5.2 Derive an expression for the magnitude 
of the magnetic field at the centre of a regular 
hexagonal loop with sides of length L, carrying 
current I. 


Exercise 5.3. (a) Show that the axial field at one 
end of a long solenoid is approximately half that in the 
middle. 


(b) Show that the axial field at the middle of a 
solenoid that is ten times longer than its diameter 
differs from the value for a solenoid of infinite length 
by less than 1%. 


Exercise 5.4 (a) Find the magnetic field 
corresponding to the magnetic vector potential 


A=—5a(y’ez + 2%ey + rez). 


(b) What is the current density associated with this 
magnetic field? 


Exercise 5.5 A magnetometer is used to measure the 
radial component of the magnetic field outside the 
(perfectly spherical) head of a patient. A current 
dipole is located a distance d directly below the 
highest point of the sphere (head) and the radial field 
is measured just outside the scalp, directly above the 
highest point. What is the measured field due to the 
combination of dipole and return currents in the 
following three cases? 


(a) The dipole is perpendicular to the surface of the 
sphere at its closest point. 


(b) The dipole is parallel to the surface of the sphere at 
its closest point. 


(c) The dipole is at 45° to the surface of the sphere at 
its closest point. 


Exercise 5.6 For the situation described in Exercise 
5.5b, describe how the radial magnetic field varies 
over the scalp in the region around the highest point of 
the head, and indicate what information can be 
deduced from the measured radial field component. 


Book 2 Chapter 6 


Exercise 6.1 An electron is travelling at velocity 
v = 2.0 x 10°ms~1(e, + e,) when it passes 
through a point where the electric field is 

E = 2.5 x 10°Vm 1(e, + 2e,) and the magnetic 
field is B = 2.5T (e; + e,). What is the magnitude 
of the resultant force on the particle at this point? 


Exercise 6.2 A proton is released from rest at time 
t = 0 in a region where both the electric field and the 
magnetic field are uniform and steady. Use the 
following information to determine the magnitudes 
and directions of the fields. 


(i) The initial acceleration of the proton is in the 
+z-direction, and has magnitude 
ag = 9.6 x 10'4*ms~?. 


(ii) At a subsequent time f,, the proton is travelling at 
speed v; = 2.0 x 10’ms7! in the —y-direction 

and has an acceleration of magnitude ag in the 
—ax-direction. 


(iii) The proton’s motion takes place in a plane with 
constant z. 


Exercise 6.3. Velocity selectors are used in particle 
accelerators to remove particles from a high energy 
beam that do not have a specified velocity. These 
devices use a combination of electric and magnetic 
fields, which you may assume are uniform in the 
region of interest, to deflect particles that are not 
travelling at the required velocity, and to allow 
particles with the required velocity to pass undeflected 
through an exit aperture. If the beam is travelling in 
the e,-direction and the magnetic field is 0.50 T in the 
e,-direction, what electric field is required to select 
particles travelling at 10% of the speed of light? 


Exercise 6.4 A magnetic field of magnitude roughly 
Beal = 2 X 10—1!° T is believed to exist throughout our 
galaxy. Suppose that a cosmic ray proton travelling at 
about 10% of the speed of light is trapped in this 
uniform field. What are the period and the diameter of 
its orbit? (Assume that a particle travelling at 10% of 
the speed of light can be treated using non-relativistic 
mechanics.) 


Exercise 6.5 Mass spectrometers separate ions 
according to their velocity and charge-to-mass ratio. 
In a particular design of mass spectrometer, ions 


travelling at velocity v = ve, ina field-free region 
are injected at the origin of a Cartesian coordinate 
system into a region in which uniform electric and 
magnetic fields act in the e,-direction, that is, 

E = Ee, and B = Be,. The ions are deflected by the 
fields and fall onto a detector grid in the plane z = 0. 


(a) If the mass of an ion is m, and its charge Q, 
show that its equation of motion can be written as 
d?x QBdz dy _ QE d?z — QBdx 


dm dt’ dt  # m’ dt? m dt. 


(b) Hence show that 


(c) Show that an ion of mass m, charge Q and initial 
velocity v is collected at the point on the detector grid 
(in the z = 0 plane) with coordinates 


Imv r2mE 0 
QB’ 2QB?’ 


Exercise 6.6 An electron spirals around the 
magnetic field lines produced by a long straight wire 
carrying 100A in a large vacuum chamber. The wire 
is located on the z-axis of cylindrical coordinates, and 
the electron has a velocity component in the direction 
of the magnetic field of ug = vj = 1.0 x 10°ms-t 
and a velocity component perpendicular to the field 
direction of v, = 1.0 x 10°ms~". Starting from 
Equations 6.14 and 6.20 of Book 2, derive expressions 
for the drift velocity of the electron (i) due to the 
curvature of the magnetic field lines, and (ii) due to 
the variation of field strength with position, in terms 
of I, v1, vl and fundamental constants. Hence 
calculate the net drift speed for the given data. 


Exercise 6.7 A magnetic bottle is axially symmetric 
about the z-axis of a cylindrical coordinate system. 
Near the middle of the bottle the magnetic field is 
uniform, has a magnitude of 2.0 x 107° T and is in 
the + z-direction. At t = 0, an electron at a point on 
the plane z = 0 has velocity v = vg eg + vz ez, where 
ug = 3.0 x 10° ms~! and v, = 2.0 x 10°ms!. 


(a) What is the initial radius of the electron’s orbit? 


(b) What is the magnitude of the magnetic field at 
the mirror point for this electron? 


(c) What is the radius of the electron’s orbit at the 
mirror point? 


(d) What is the velocity of the electron when it 
returns to the plane z = 0? 


Book 2 Chapter 7 


Exercise 7.1. Two parallel-sided slabs, 1 and 2, of 
materials with low conductivities, 0; and 09, have 
thicknesses, d, and dz, and the same area. One face of 
each slab is in good electrical contact with a face of 
the other slab, and the two external faces of the 
combined slab are coated with thick layers of a 
high-conductivity metal. 


If the potential difference between the metal layers is 
AV, with the metal layer on slab | at Earth potential, 
V = 0, neglecting any edge effects, what is the 
potential at the interface between the two materials? 
Check that the expression you derive gives reasonable 
results for limiting values of the conductivities and 
thicknesses of the slabs. 


Exercise 7.2 Electric currents are set up in a steel 
tube in three different ways. 

(a) The ends of the tube are clamped between two 
copper blocks, and the current flows along the length 
of the tube. 

(b) The tube is filled with a high-conductivity fluid, its 
ends are capped with insulating material, and the tube 
immersed in a high-conductivity fluid. A current is 
passed between the two fluids through the walls of the 
tube. 

(c) A magnetic field that increases linearly with time 
is applied parallel to the axis of the tube so that a 
current is induced to flow around the wall of the tube. 


Calculate the resistance to current flow in these three 
situations for a tube that is 10 cm long, has radius 

5.0 mm and wall thickness 0.50 mm, and conductivity 
10410) 0 mm, 


Exercise 7.3. Solve the problem illustrated in Figure 
7.9 of Book 2, by deriving an expression for the 
resistance between a small sphere, radius a, embedded 
in a large block of poorly conducting material, and a 
metal coating on one face of the block. The sphere 

is at a distance d from the metal surface. Hence 
calculate the resistance for a semiconducting block, 

o =5.00-!m“|, if the ratio d/a = 10. 


Exercise 7.4 Consider a situation similar to the 
‘jumping ring’ demonstration described in Exercise 
7.6 and Figure 7.15 on p. 157 of Book 2, but with no 
iron core in the solenoid. What forces act on the 
aluminium ring when the current is switched on (a) 
when the ring is horizontal and just above the centre 
of the solenoid, and (b) when the ring is horizontal 
and just above the top end of the solenoid. 


Exercise 7.5 A circular coil, with NV turns, area A, 
rotates at frequency f about a vertical axis that lies 
along a diameter of the coil. (See Figure 6). The 
amplitude of the voltage across the terminals of 

the coil is Vo. Derive an expression for the 
magnitude Byori, of the field in the horizontal plane. 


Hence determine the value of Byoyiz if N = 50, 
A= 0.50 x 10-4 m?, f = 80 Hz and V = 25 mV. 


es 


Figure 6 For Exercise 7.5. 


Exercise 7.6 (a) Twocoils, A and B, are in 
close proximity to each other. When a current 

I, = I,gsinwt flows through coil A, an induced 
current Iz = —Ipo cos wt flows through a large 
resistance Rp connected across the terminals of coil 
B. Derive an expression for the mutual inductance 
M between the two coils in terms of parameters 
introduced above. (You should assume that 

Rp > wLp, where Lp is the self-inductance of coil 
B.) 


(b) The current source is disconnected from coil A 
and a large resistance Ra is connected between the 
terminals of the coil. If a current Iz = Ij, exp [—at] 
flows through coil B, what is the current in coil A? 
(Assume that Ra >> wl a.) 


Exercise 7.7. An inductor core is made by drilling a 
hole, radius a, through the centre of a disc that has 
thickness a, radius 2a and relative permeability pw. A 
toroidal coil with N turns of wire is tightly wound on 
the core. Derive an expression for the self-inductance 
L of the inductor, and hence calculate the value of L if 
the inductor has 100 turns, with a = 4.0 mm and 

ps = 500. (Assume the core is an LIH material.) 


Book 2 Chapter 8 


Exercise 8.1 Four equal charges, g, are located at 
the points (a, 0, 0), (0, a, 0), (—a, 0, 0), (0, —a, 0). 
(see Figure 7). What is the electrostatic energy of this 
arrangement? 


z= 0 plane 


oO © 


RY 


oO 


Figure 7 For Exercise 8.1. 


Exercise 8.2. (a) What capacitance C' is required to 
store 5.0 MJ, if the maximum voltage difference 
across the capacitor is 2.0kV? 


(b) What inductance L is required to store 5.0 MJ, if 
the maximum current through the inductor is 2.0kA? 


(c) What are the time constants for charging the 
capacitor through a resistance R of 0.2 Q and for 
establishing the current in the inductor through the 
same resistance? 


(d) What is the resonant frequency f, of a circuit in 
which components with these values of C’, L and R 
are connected in series? 


Exercise 8.3. What is the ratio of the electrical 
energy to the magnetic energy stored in unit volume 
of the atmosphere just above the Earth’s surface, 
where the electric field strength is 100 V m7! and the 
magnetic field strength is 5.0 x 107° T? 


Exercise 8.4 An electronic flash unit contains a 
1000 jxF capacitor and is charged using a 300 V power 


supply. 
(a) It takes 5s to charge the capacitor to 95% of the 


maximum charge. What resistance F is connected in 
series with the capacitor in the charging circuit? 


(b) 90% of the energy stored in the capacitor is 
discharged in 2ms. What is the average power 
dissipated during this period? 


Exercise 8.5 Show that the energy U stored in 
the electric field in the free space surrounding a 
conducting sphere, radius R, carrying charge Q, is 
given by 
2 
U= @ : 
8a eogR 

How does this expression compare with the electric 


potential energy of the charge on a conducting sphere, 
determined in Exercise 8.2 in Book 2? 


Exercise 8.6 In Worked Example 8.1 of Book 2, 
we showed that the electric potential energy of a 
uniformly charged sphere, with radius R and charge 


density p, is given by 
_ Anp?.R° 
a 15é9 


(a) Obtain an expression for U in terms of the total 
charge Q in the sphere. 


(b) Write down an expression for the energy stored 
in the electric field in the region outside the sphere. 
(Hint: you can make use of a result from the previous 
question.) 


(c) Explain why the expressions in parts (a) and 
(b) are different, and confirm your explanation 
by independently deriving an expression for the 
difference in energy. 


Exercise 8.7 A naive inventor plans to store energy 
by passing a current through a large air-filled solenoid, 
1.0m long, 10cm radius. The solenoid will be 

wound with 200 turns of copper wire, radius 2.5 mm, 
conductivity 6.5 x 10°Q-!'m!. Calculate the energy 
that would be stored by this solenoid when carrying a 
current of 100 A and the power dissipated, and hence 
comment on the feasibility of this proposal. 


Book 2 Chapter 9 


Exercise 9.1 A circular coil, with N turns of 

radius a and with self-inductance L, is wound from 
superconducting wire. The start of the first turn is 
connected to the end of the Nth turn to form a 
continuous superconducting path. The coil is cooled 
below its critical temperature in a uniform applied 
magnetic field Bo, which is perpendicular to the plane 
of the coil, and the applied magnetic field is then 
reduced to zero. 


(a) What is the magnetic flux through the coil after 
the applied field is reduced to zero? 


(b) What is the current in the coil? 
(c) What is the field at the centre of the coil? 


(d) What can you say about the magnetic field in the 
plane of the coil at distance 0.95 from its centre? 


Exercise 9.2. (a) Estimate the temperature 
below which a long thin cylinder of indium is 
superconducting in a magnetic field with strength 
B =0.014T, aligned parallel to its axis. (You will 
need to use data from Chapter 9 of Book 2.) 


(b) Estimate the critical current of a long straight 
indium wire, radius 1.0 mm, at 7 = 1.2K in the 
Earth’s magnetic field. 


Exercise 9.3. Two superconducting discs with 
identical dimensions are made from two different 
superconductors. Disc A is made from a material 
that has a penetration depth A, = 40 nm and a 
coherence length of £4 = 800 nm, whereas disc B has 
a penetration depth of Ag = 80 nm and a coherence 
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length €g = 4nm. Each disc is placed in turn in a 
solenoid, with the plane of the disc perpendicular to 
the solenoid’s axis, and the magnitude of the field in 
the solenoid is steadily increased. Describe what 
happens to the magnetic field in and around the 
discs as the magnitude of the field in the solenoid is 
increased. 


Exercise 9.4 Produce a table that briefly contrasts 
type-I and type-II superconductors. 


Exercise 9.5 A long thin superconducting cylinder 
with radius F is aligned with its axis parallel to 

the direction of a uniform applied magnetic field 
Bo = Boze:. 


(a) Using Equation 9.12 of Book 2 as a guide to 
how the magnetic field varies near to the surface 

of a superconductor, write down an expression 

in cylindrical coordinates for the field within the 
superconducting cylinder. Note any assumptions that 
you make. 


(b) Use Ampére’s law to derive an expression for the 
current density within the cylinder. 


(c) Derive an expression for the screening current 
per unit length flowing near the surface of the 
cylinder, and hence confirm that this current produces 
a uniform magnetic field within the long cylinder that 
cancels the applied field. 


Exercise 9.6 The mean magnetic field in a certain 
type-II superconductor in the mixed state is 0.50 T. 
Given that the flux within each normal core is equal to 
the flux quantum (2.07 x 10~!° T m?), what is the 
number density of cores passing through a plane 
perpendicular to the field direction? What is the 
average separation between neighbouring cores? 


Book 2 Chapter 10 


Exercise 10.1 Show that E? — c? B? is an invariant 
quantity. 


Exercise 10.2 In the laboratory frame of reference, 
F, there is an electric field of magnitude F in the 
+2-direction and a magnetic field of magnitude B 
inclined at an angle 6 to the +-direction. 


(a) What are the electric and magnetic fields in a 
frame F’ that is travelling at uniform speed v in 

the x-direction relative to the laboratory frame? 
Determine the inclination of the magnetic field to the 
+2-direction in the frame F’. (To simplify the 
derivation, assume that the magnetic field has no 
component in the z-direction.) 


(b) If @ = 30°, at what speed v must frame F’ travel 
relative to frame F for the magnetic field observed in 
this frame to be at 60° to the 2’-axis? 


Exercise 10.3. Which of the following are invariant 
quantities? 


e the rest mass of an electron; 
e the charge of an electron; 


e the speed of an electron; 


the energy of an electron; 

e the charge density of a cloud of electrons; 
e the speed of light in vacuum; 

e E-B. 


Exercise 10.4 Throughout a certain region of 

space there is a stationary distribution of positive 
charge, with uniform density p,, and a distribution of 
negative charge, with uniform density p_, that is 
travelling with speed v in the +2-direction. The net 
charge density is zero, that is p, = —p_. 


(a) What are the current densities, J, and J_, due 
to the positive charge distribution and due to the 
negative charge distribution? 


(b) What are the charge densities and current 
densities due to the positive charge distribution and 
due to the negative charge distribution measured 
by an observer who is travelling at speed v in the 
+2-direction? What are the net charge density and 
current density according to this observer? 


Exercise 10.5 An infinite cylinder with 
cross-sectional area A and uniform charge density p is 
at rest along the x-axis of an inertial frame F. A test 
charge Q is stationary at point (0, d, 0) outside the 
cylinder. 


(a) What are the electric and magnetic forces acting 
on charge Q according to an observer in frame F? 


(b) A second frame, F’, is in standard configuration 
with frame F, and travels at speed v along its x-axis. 
What are the charge density and current density 
associated with the cylinder according to an observer 
in frame F’? 


(c) What are the electric and magnetic fields 
generated by the cylinder according to the observer in 
frame F'? 


(d) What are the electric and magnetic forces on test 
charge Q according to an observer in frame F’, and 
what is the total force on Q? 


(e) Check that the net force in frame F’ determined 
in (d) is the same as that predicted by the force 
transformation relationships in Equations 10.17—10.19 
of Book 2. 


Exercises for Book 3 


Book 3 Chapter 1 


Exercise 1.1 (a) Write down an expression for the 
physical electric field of a plane wave that is travelling 
in the direction (e, +e, + e-)/V3 with amplitude 


Eo, wavenumber k, frequency f, phase shift 0 and 
polarization in the direction (e,, — ey) / V2. 
(b) Confirm that your expression represents a 
transverse wave. 
(c) Write down another possible linear polarization 
for a transverse wave travelling in the direction 
specified in part (a). 
Exercise 1.2. The expression 

E = acos (362 + 26y + Gz + yt) (ez — 3ey) 
represents the electric field of a plane monochromatic 
electromagnetic wave. If the values of the constants 
in this expression are ~ = 2.0 x 10? V m“1, 
(3 = 6.0 x 10°m~! and y = 6.7 x 10!°s—!, what are 
the wavenumber, frequency and amplitude of this 
wave? Confirm that the wave travels at speed c, the 
phase speed for electromagnetic waves in a vacuum. 


Exercise 1.3. An electromagnetic wave is travelling 

in the direction 4(2e, + 2e, — ez) and is polarized in 
the direction (e, + 2e.)/./5. In what direction is the 

magnetic field associated with this wave? 


Exercise 1.4 What types of waves are produced by 
combining the plane wave 


E = Egexp{i(kx — wt)] ey 
in turn with each of the following plane waves? 
(a) Eoexpli(kx — wt)] ey 
(b) Epo expli(kx — wt)]e, 
(c) 2£oexpli(kx — wt)] ez 
(d) Eoexpli(kx — wt + 1/2)] 
(e) Eogexpli(kx — wt — 1/2)| 
(f) 2£o expli(kx — wt + 1/2)]e, 
(g) Eoexpli(ka — wt — 7)Jey 
(h) Eo expli(ka + wt)] ey 
(i) Eoexpli(1.01lka — 1.01wt)] e, 


Exercise 1.5 A powerful research laser produces 
pulses of light, wavelength 527 nm, containing 300 J 
of energy in 1.0 x 10~°s, and concentrated into an 
area with diameter 1.0 x 10~4m. Assume that the 
pulse is uniform over its area and uniform in time over 
its nanosecond duration, and that it is polarized in the 
x-direction and is travelling in the +y-direction 
through a vacuum. 


e, 
ez 


(a) Calculate all of the following quantities, which 
specify the radiation during the pulse: w, ky, ky, kz, 
Fox, Evy, Foz, Box, Boy and Boz- 

(b) Given that air is ionized by an electric field 
strength exceeding 3 x 10° V m~!, explain why the 
beam must travel through a vacuum. 


Exercise 1.6 A linearly polarized plane wave 
travelling in the horizontal direction has its electric 
field, amplitude Eo, in the vertical direction. The 
wave passes through a polaroid filter that has its axis 
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of polarization at 30° to the vertical, and is then 
incident on a second polaroid filter. What is the 
amplitude of the electric field transmitted by the 
second filter when the axis of polarization of this filter 
is (a) vertical, and (b) horizontal? Assume that the 
filters are ‘ideal’, that is, they are perfectly transparent 
to waves polarized parallel to the axis of polarization 
and perfectly opaque to waves polarized perpendicular 
to the axis of polarization. 


Exercise 1.7. (a) In Chapter 1 of Book 3 we 
showed that the time-averaged value of the magnitude 
of the Poynting vector for a monochromatic plane 
wave is N = seoke. Show that an alternative 
expression is N = 5 Bec/ po. 


(b) Ifthe power per unit area of a radio wave in the 
region of a radio receiver is 10-'° W m2, what is the 
amplitude of the magnetic field associated with the 
wave? 


Book 3 Chapter 2 


Exercise 2.1 For each of the following oscillating 
current elements, state whether it can be regarded as a 
Hertzian dipole, and for those that can, indicate how 
far away from the dipole the point of measurement 
must be if the fields are to be well represented by the 
radiation field alone. 


(a) Dipole length 10 km, oscillating at 50 Hz. 
(b) Dipole length 10m, oscillating at 100 MHz. 
(c) Dipole length 1 cm, oscillating at 3 GHz. 


Exercise 2.2. Show that Equation 2.17 of Book 3 
reduces to Equation 2.1 in the limit when w — 0 andr 
is small. 


Exercise 2.3 The vector amplitude of the 

electric field at a point with spherical coordinates 
(1.0m, 7/2, 0) at time t = 0 due to a Hertzian dipole, 
frequency 3 GHz, located at the origin and aligned 
with the 9 = 0 direction, is 2.0 x 107° Vm~! ep. 


(a) What is the vector amplitude of the electric field 
at (2.0m, 7/4, 7/2)? 


(b) What is the vector amplitude of the magnetic 
field at (2.0m, 7/4, 7/2)? 


Exercise 2.4 A Hertzian dipole is located at the 
origin of Cartesian coordinates and is aligned with the 
y-axis. The strength of the dipole is 3.0 x 10-° Am, 
and its angular frequency is 2.0 x 10° s~. Calculate 
the time-averaged value of the Poynting vector at the 
point with coordinates (2.0, —2.0, 1.0) m. 


Exercise 2.5 The mean scattering cross-section for 
blue light with \ = 400 nm by air molecules is 

1.50 x 10~-®° m?. Determine the percentage of the 
incident power that is scattered when blue light passes 
through 1.0 km of air at sea-level (P = 1.0 x 10° Pa) 
at 0 °C, and compare this with the percentage scattered 
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for red light, \ = 600 nm, under the same conditions. 
You should assume that there is negligible scattering 
by water droplets and dust particles. 


Exercise 2.6 Assume that you are lying on your 
back on a beach under a cloudless sky with the Sun 
directly overhead and your feet pointing towards the 
south. A tiny sunshade blocks direct light from 

the Sun reaching your eyes, but otherwise all of 

the sky is visible. You use a polarizing filter that 
completely blocks the component of the incident 
light with polarization parallel to one edge of the 
filter but transmits the component with polarization 
perpendicular to this edge. Explain what you would 
expect to observe when looking at the sky through this 
filter. 


Exercise 2.7 A Hertzian dipole is located at the 
origin of spherical coordinates and is aligned with the 
direction 0 = 0. 


(a) What is the ratio of the time-average power 
flow per unit area at a point with coordinates 
(50m, 90°, 0) to that at (100 m, 60°, 90°)? 


(b) What fraction of the total power is radiated in the 
range 70° < @ < 110°? 


Book 3 Chapter 3 


Exercise 3.1 Electromagnetic radiation with angular 
frequency w = 27 x 10°s~! travels through a 
material with relative permittivity ¢ = 4.0. What are 
the values of 


(a) the frequency, 
(b) the phase speed, 
(c) the wavelength, and 


(d) the wavenumber for the radiation in this 
material? 


Exercise 3.2. An electromagnetic plane wave 
travelling in free space, with electric field given by 


E; = Eo expli(ky — wt)] ez, 


is normally incident on a boundary, in the plane y = 0, 
of a dielectric material with relative permittivity e. 


(a) Write down expressions for the magnetic field 
associated with this electric field and for the electric 
and magnetic fields of the waves transmitted and 
reflected at the boundary. You should express these 
fields in terms of Eo, k, w, < and the ratios r, and t,, of 
the amplitudes of the reflected and transmitted electric 
fields to the amplitude of the incident electric field. 


(b) By applying the boundary conditions for the 
electric and magnetic fields at the interface in the 
plane y = 0, derive expressions for the amplitude 
reflection ratio and the amplitude transmission ratio 
for normally incident radiation. 


Exercise 3.3. A plane electromagnetic wave 
travelling through free space is normally incident on 


the surface of a dielectric material with relative 
permittivity 1.7. 


(a) What is the value of the refractive index of the 
material, and what is the value of the reflectance of the 
surface for this wave? 


(b) If the amplitude of the electric field of the 
incident wave is 5.0 V m~', determine the amplitudes 
of the electric fields of the reflected and transmitted 
waves, and show that your results satisfy the boundary 
condition for E. 


Exercise 3.4 What percentage of the power of 
normally incident light is transmitted through the 
glass wall of a fish tank filled with water (a) when the 
light travels from the surrounding air into the water, 
and (b) when the light travels from the water into the 
air? Assume Ngjass = 1.50 and nwater = 1.33. 


Exercise 3.5 Linearly polarized light travelling 
through air is incident at 45° to the normal to a 
plane face of a diamond, which has refractive index 
n = 2.42. The polarization direction is at 45° to the 
scattering plane. What is the total reflectance in this 
situation? 


Exercise 3.6 Circularly polarized light is incident at 
the Brewster angle on a dielectric surface. Describe 
the polarization of the light that is reflected from the 
surface. 


Exercise 3.7 A thin parallel-sided glass plate with 
refractive index 1.55 is sandwiched between large 
blocks of glass with refractive index 1.50. An infrared 
beam travels in the thin plate at an angle 6, to the 
normal to its surface. 


(a) For what range of values of 6; will the wave be 
essentially confined within the plate? 


(b) For an infrared beam with a free space 
wavelength of 1500 nm that travels at 80° to the 
normal, what is the characteristic distance that the 
electric field penetrates into the blocks on either side 
of the plate? 
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Exercise 4.1 A plane wave travels in the 
e,-direction through a dielectric medium with 
refractive index n = Neal + 1Nimag. What is the 
ratio of the amplitudes of the magnetic fields at the 
points with Cartesian coordinates (21, y1, 21) and 
(x2, ye, 22)? 

Exercise 4.2 The physical electric field associated 
with a plane wave in a dielectric material has the form 


Epnys = Ep exp[—a1z] cos (agz — azt) ey. 


Express each of the following quantities in terms of 
symbols that appear in the expression for Epnys and 
the speed of light in free space c: 


Kereal; Kimag, W, MMreal, Timag, €real; €imag; 


Uphase, absorption length. 


Exercise 4.3. The relative permittivity of amorphous 
silica for infrared radiation with free space wavelength 
2000 nm is 1.46 + 11.0 x 107"°. Its relative 
permeability is 1.0. 


(a) Determine the values of 
(i) the refractive index, 
(ii) the wavenumber, 

(ili) the phase speed, 


(iv) the factor by which the amplitude of a plane 
wave decreases over a distance of 1.0 km. 


(b) Use the results of part (a) to write down 
expressions for the physical electric and magnetic 
fields of a plane wave with free space wavelength 
2000 nm travelling through amorphous silica in the 
x-direction, polarized in the y-direction. Assume that 
the amplitude of the electric field oscillation at the 
origin is Eo. 

Exercise 4.4 The behaviour of a particular dielectric 
material is represented by the simple classical model 
described in Section 4.2 of Book 3, with plasma 
frequency w, friction parameter y and natural 
frequency w,, related by 


Y=Wpy, and Wn = 4.0x10' s—t. 


Determine the values of the following quantities: 


Wy = 2.00wy, 


(a) the real and imaginary parts of the dielectric 
function in the limit w — 0, 


(b) the real and imaginary parts of the dielectric 
function in the limit w — oo, 


(c) the angular frequency at which é;ca) = 1, 


(d) the value of cimag at w = wy. (This is 
approximately the angular frequency at which éimag 
has its maximum value.) 


Exercise 4.5 (a) A steak has relative permittivity 
€ = 40+ 112 at 2.45 GHz, the frequency used in 
microwave ovens. Estimate the depth over which the 
amplitude of the electric field inside the steak falls to 
half of the external value. 


(b) Polystyrene has ¢ = 1.03 +i3 x 107° at 
2.45 GHz. Explain why it is a suitable container for 
use in a microwave oven. 


Exercise 4.6 Blue light is more strongly refracted by 
a glass prism than red light. 


(a) What is the sign of dn/dw for glass at optical 
frequencies? 


(b) Is the group speed greater or less than the phase 
speed for light in glass? 


(c) Fora short pulse of light travelling through glass, 
would the individual maxima and minima of the 
electric field move from the rear to the front of the 
pulse or from the front to the rear? 


Exercise 4.7 For many types of glass used to make 
lenses, the refractive index is well-represented by the 
empirical expression 


B 
=A+-—, 
n + x 
where Xo is the free space wavelength, and for hard 
crown glass (type K5) the values of the constants are 
A = 1.522 and B = 4.59 x 107° m?. 


Determine each of the following quantities for both 
blue light with free space wavelength Apiie = 400 nm 
and red light with A;eq = 600 nm: 


(a) refractive index; 
(b) phase speed; 
(c) group speed; 


(d) the transmission angle (angle of refraction) when 
light is incident on the glass surface in air at 30° to the 
normal. 


(For the speed of light, use the more precise 
value than is usually used in this course, i.e. c, 
2.998 x 108 mst, .) 


Exercise 4.8 Use data from Chapter 4 of Book 3 to 
estimate the absorption lengths in pure water and in 
amorphous silica for electromagnetic radiation with 
free space wavelengths of 500 nm, 1000 nm and 
2000 nm. 


Exercise 4.9 Pulses of light are launched into an 
optical fibre from a modulated light source that 
produces square pulses of duration 1 ys, separated by 
10 us. The source emits light in two narrow bands 
centred on free space wavelengths of 600 nm and 
650 nm and with spreads of about 5 nm around these 
wavelengths. The group speed for light in the fibre is 
given by the expression 


B —1 
Usroup = € (4 + x) ) 
0) 


where Ag is the free space wavelength, and 
A = 1.4580 and B = 1.062 x 10714 m?. 


For reliable communication, the separation of the two 
wavelength components of each pulse is required to 
be less than 10% of the 10 jus separation between 
successive pulses. What is the maximum length D 

of fibre that can be used if this condition is to be 
satisfied? 


Book 3 Chapter 5 


Exercise 5.1 A simple classical model predicts that 
the relative permittivity function for a conductor has 
the form 


2 
Wp 


a w* +iw/T. 


(a) What is the significance of the terms wp and 7;? 


14 


(b) What is the condition for e(w) to be real? 


(c) What is the value of e(w) for very high 
frequencies, and what is the physical significance of 
this result? 


(d) What can you say about the values of e(w) 
for frequencies in the range 1/T < w < wp», and 
what consequence does this have for the nature of 
electromagnetic disturbances in this range? 


(e) In the low frequency limit, derive an expression 
for ¢(w) that involves the conductivity o. 


Exercise 5.2. Compare the forms of the electric and 
magnetic fields just below the plane surfaces of a 
dielectric material and a metallic conductor when an 
electromagnetic plane wave travelling through air is 
normally incident on the surfaces. 


Exercise 5.3 A thin metallic strip on a circuit 

board has length L, width a and thickness t, with 

L > a> t. Derive an expression for the ratio of the 
resistance between the ends of the strip at frequency f 
to its DC resistance, assuming that the skin depth is 
small compared with the thickness t. 


Exercise 5.4 How many TE modes with 
w = 10 GHz can propagate between two parallel 
conducting planes separated by 10 cm? 


Exercise 5.5 Use Faraday’s law to derive 
expressions for the magnetic field components that 
accompany the electric field given by the expression 
in Equation 5.34 of Book 3. 


Exercise 5.6 Starting from the dispersion relation 
for the TE mn mode in a rectangular waveguide with 
internal cross-section a x b, with a > b, derive an 
expression for the group speed of this mode. 


Exercise 5.7 A rectangular waveguide has internal 
cross-section with dimensions 3.4cm x 1.7cm. 


(a) What are the lowest four cut-off frequencies for 
TE modes in this waveguide? 


(b) Show that the phase speed of the TE,,,,, mode is 
given by 


f2 —1/2 
Uphase = C (1 = fn . 
and hence determine the phase speeds and groups 


speeds for the four modes identified in part (a) when 
jf = 10GHz, 


Book 3 Chapter 6 


Exercise 6.1 The number density of electrons in a 
collisionless plasma is 10!° m~?. What are the 
values of the group speed and the phase speed in this 
plasma for electromagnetic waves with frequencies 
1.00 MHz, 2.00 MHz and 100 MHz? 


Exercise 6.2 A pulsar, which is 10!9 m from Earth, 
emits pulses of electromagnetic radiation across a 


broad spectrum. The pulses travel to Earth through the 
interstellar medium, which can be regarded as a dilute 
plasma within which the number density of electrons 
and ions is about 3 x 104m~?. What is the difference 
in arrival time at Earth of the pulse detected with a 
radio telescope sensitive to frequencies around 

100 MHz and the corresponding pulse detected with 
an optical telescope with a red filter? 


Exercise 6.3. An electromagnetic plane wave with 
angular frequency w is travelling in the z-direction 
and is normally incident on a dielectric window ina 
chamber containing a collisionless isotropic plasma 
that has plasma frequency wy) = 27 x 10°s~!. Derive 
an expression for the ratio of the amplitude of the 
electric field in the plasma at a distance z from the 
window to the amplitude just inside the window when 
Ww < Wp, assuming that the field is constant across a 
plane with constant z within a ‘beam’ defined by 

the area of the window. Hence evaluate this ratio 

for a distance z = 0.5m for angular frequencies 
0.99w,, 0.90w, and 0.50w,, and its limiting value 
when w/wp <1. 


Exercise 6.4 A plasma used for etching is heated by 
coupling 13.56 MHz radiowaves into the plasma in an 
applied magnetic field. The plasma frequency is 

Wp = 8.0 x 10'° s~! and the cyclotron frequency is 
i = 3.0% 10 5. 


(a) What are the values of the electron number 
density and the magnetic field strength in the plasma? 


(b) What type of waves are produced in the plasma, 
and what is their group speed? 


Exercise 6.5 Microwaves with wavelength 3.0cm 
are to be used to heat a plasma. 


(a) What steady magnetic field strength is required 
for electron cyclotron resonance to occur? 


(b) Use information from Figures 6.7 and 6.9 of 
Book 3 to explain why it would not be efficient to heat 
a plasma that is at a pressure of 104 Pa with this 
microwave source. 


Book 3 Chapter 7 


Exercise 7.1 State whether the changes listed below 
would increase, decrease or leave unchanged the 
scattering cross-section of a collagen fibril embedded 
in the transparent matrix of the cornea, and give a 
reason to support each answer. Assume that each of 
the changes, apart from (f), is small. 


(a) Increase in the fibril diameter. 

(b) Increase in the refractive index of the fibril. 
(c) Increase in the refractive index of the matrix. 
(d) Increase in the spacing of the fibrils. 

(e) Increase in the wavelength of the radiation. 


(f) Changing the polarization of the incident 
radiation from parallel to the fibril to perpendicular to 
the fibril. 


Exercise 7.2. The fundamental resonance frequency 
for a hairpin resonator in a vacuum is fp = 500 MHz. 
At what other frequencies might resonances be 
observed for this hairpin in a vacuum? 


Exercise 7.3 A (long) hairpin resonator has its 
fundamental resonance frequency in vacuum at 

fo = 100 MHz. What will the fundamental resonance 
frequency be when the hairpin is immersed in (a) 
water (€ = 80), (b) ethanol (€ = 25) and (c) a cold 
collisionless plasma with electron number density 
2.0: 10 m- 87 


Exercise 7.4 The fundamental resonance frequency 
of a hairpin resonator, length 10.0 mm, immersed in a 
dilute cold collisionless plasma is twice the value of 
the fundamental resonance frequency when the 
hairpin is in a vacuum. What is the number density of 
electrons in the plasma? 
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Solutions for Book 1 


Book 1 Chapter 1 


Solution 1.1 According to the suggested model, 
the person contains 10 x 60kg/(3.0 x 10~7° kg) 
electrons, each of charge —1.6 x 10~!°C, giving 
a total electronic charge of —3.2 x 10°C. This 
is normally exactly balanced by the positive 
charge of protons — so, to have a net charge of 

1 pC, the person would need to lose a fraction 
10-°/3.2 x 10° = 3.1 x 107° of their electrons. 


Solution 1.2 The charge-to-mass ratio of a proton is 


e  160x10-%C 
m 1.67 x 10-27kg 


= 9.58 x 10’Ckg!. 


So one gram of protons has charge q = 9.58 x 10*C. 
Coulomb’s law gives an electrostatic force of 
magnitude 

ees 

— Aneg r2 


9.58 x 10*C)? 
= 8.99 x 10°Nm?2.C-? x (9:98 x 10 CY" 
. m* — * Ra x 108m)? 


= 560N. 


This is a significant force, in spite of the enormous 
separation of the charges! Of course, it would not be 
easy to gather one gram of protons together in a small 
region — the repulsive forces would be huge. 


Solution 1.3 The charge g experiences no force 
because it is midway between two identical charges. 
The two outer charges repel one another with a force 
of magnitude Q?/47eod’. The charge q must provide 
a balancing attraction so 


1 Q 1 qQ = 
Ar eg a AT EQ (d/2)? ; 


which gives g = —Q/4. The equilibrium is not stable. 
For example, if g is moved slightly to the left, the 
attraction from the left-hand charge exceeds that from 
the right-hand charge so q tends to move even further 
to the left rather than returning to its equilibrium 
position. 


Solution 1.4 (a) Charges at opposite corners 
produce opposite electric fields at the centre of the 
square, so the total electric field is zero. (b) When one 
of the corner charges is removed, the pair of charges 
in opposite corners still produce opposite electric 
fields at the centre of the square and therefore cancel 
out. The third charge produces a field of magnitude 


Be 1 q — 1 24 
Amey (d//2)2 Ameo a?’ 
This is the total field at the centre of the square. 


Solution 1.5 The problem can be simplified by 
noting that the two 2q charges produce opposite fields 
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at the centre of the square, so their contributions 
cancel out. We are left with: 


_ 1 q Cy — Cy 4q €y — ey 
© aa (EE a 
_ 5V2 ¢ 


= eee (Gn = ey): 


Solution 1.6 By symmetry, the electric field 
vanishes at the centre of the disk. The field due to any 
element of the disk is exactly cancelled by the field 
due a corresponding element on the opposite side of 
the disk. 


Solution 1.7 (a) By symmetry, the total electric field 
vanishes at the centre of a uniformly-charged ring. 
The contribution from each small element of the ring 
is cancelled by the contribution from a similar element 
at the opposite end of the diameter passing through 
the centre of the ring. 


(b) This part of the question is most easily answered 
using the principle of superposition. Let the electric 
field due to the complete ring be E,ing, the electric 
field due to the removed segment be Exeg and the 
electric field due to the almost complete ring be Eye. 
Then the principle of superposition tells us that 


Exing = Eyseg + Erem.- 
At the centre of the ring we have seen that E,jng = 0, 
so 

Erem = —Egeg. 


The charge of the segment is Q1/27R (the 
total charge of the ring times the fraction of the 
circumference of the ring occupied by the segment). If 
the segment is short, it behaves like a point charge. At 
the centre of the ring it produces an electric field 

1 Qi/2rRk... Qt. 

r= r 

A4neg R? 8r2e9R3 
where r is a unit vector pointing radially away from 


the segment. The almost complete ring produces the 
opposite electric field 


QI. « 
7 87 2e oR? 7 
which points radially towards the gap in the ring. 
Solution 1.8 The magnitude of the field is 


Eseg a 


Brem = 


B = /(3NC4)24 (ANC)? = 5NC! 
and its unit vector is 


~  (3e; +4e,)NC7* 

E= aco = (0.6e, + 0.8e,). 
The component of the field in the direction the unit 
vector U is given by the scalar product 

x @NC ANC 1 
EB = E “uz = 
: v2 v2 v2 


Ne. 


Book 1 Chapter 2 


Solution 2.1 Because the charge distribution is 
cylindrically symmetric, we surround the z-axis by a 
nested set of thin cylindrical tubes, extending from 
z=Oto z= Land fromr = Otor = R. A thin 
cylindrical tube with radii between r and r + Ar has 
volume 27r Ar x L, where Ar is the thickness of the 
tube. The charge density in this tube is A/r, so the 
tube contributes a charge 


AQ = 2 x Qar Ar x L = 27 AL Ar. 
iy 


The total charge Q is obtained by integrating over all 
the tubes from r = 0 to r = R, giving 


R 
Q= | 2nAL dr = 27 ALR. 
0 


Solution 2.2 (a) This situation is spherically 
symmetric, so we split the volume into thin spherical 
shells centred on the origin. The total charge of the 
electron cloud is 


Qelec = | Aexp(—2r/ao) x 4nr? dr 
0 
= ana | r? exp(—2r/ag) dr. 
0 


Changing the variable of integration to x = 2r/ag and 
using the standard integral given in the question, we 
obtain 


3 oe) 
Qelec = 47 AX (S) | a? exp(—2) da = raj A. 
0 


In order for the atom to be neutral the charge of the 
electron cloud must be opposite to the charge of the 
proton. So Qelec = —e and A = —e/7a. 
(b) The charge density of the electron cloud at the 
Bohr radius is therefore 

—e X exp(—2) 

p(a0) = ——_,— 
—1.6 x 107'8C x 0.1353 
a (5.29 % 10-4 


= —4.7x 10!°Cm-?. 


Solution 2.3 If the field were produced by an 
isolated point charge @ at the origin, Gauss’s law tells 
us that the electric flux over the any sphere centred on 
the origin is Q/éo, which is independent of the radius 
of the sphere. Surround the origin by a sphere of 
radius R. The vector field E is perpendicular to the 
surface of the sphere and has a constant normal 
component A exp(—/a) over this surface. The 
total flux of E over the surface is therefore 
Aexp(—R/a) x 47R?. This clearly depends on the 
radius of the sphere, becoming negligible for very 
large spheres. This is inconsistent with Gauss’s law. 


This argument works no matter how the charge 

is moving because Gauss’s law is valid under all 
circumstances. It does not rule out the possibility that 
the field is produced by an extended distribution of 
charge. We can imagine positive and negative clouds 
of charge centred on the origin, with the negative 
cloud tailing off more gradually than the positive 
cloud. The total charge enclosed by a sphere centred 
on the origin would then depend on the radius of the 
sphere, avoiding the conflict with Gauss’s law. 


Solution 2.4 The divergence theorem is true for all 
vector fields. 


However, the flux of F is not generally equal to zero 
through a closed surface that does not enclose the 
origin because the field dies away so quickly that 
distant parts of the surface contribute virtually no 
(positive) flux, while parts of the surface close to the 
origin contribute a sizeable negative flux. As an 
example, imagine a closed surface consisting of two 
nested hemispheres, one with a small radius r, the 
other with a large radius R >> r, both centered on 
the origin, and a flat ring in the xy-plane that closes 
the gap between the hemispheres (see Figure 8). The 
ring contributes no flux, the hemisphere with radius R 
negligible flux, while the hemisphere with radius r 
contributes a sizeable negative flux. 


a 


Figure 8 For Solution 2.4. 


ry-plane 


Solution 2.5 By symmetry, any electric field inside 
the cylindrical shell can only have a radial component 
and its magnitude can only depend on the distance r 
from the axis. In cylindrical polar coordinates, 


K(r, ?, z) _ E,(r) e;. 


Choose a cylindrical Gaussian surface, radius r, 
length 61, concentric with the cylindrical shell, but of 
smaller radius. No charge is enclosed by this surface 
so Gauss’s law gives 


E,(r) x 2ar x Al=0 


Thus £,.(7r) = 0 for all r 4 0. This argument shows 
that E = 0 at all points inside the cylinder, except 
possibly along the central axis. However, the electric 
field also vanishes along this axis for symmetry 
reasons. 
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Solution 2.6 Let the point P at which the field 

is measured be a distance R from the centre 

of the cylinder (see Figure 9). Consider two 
symmetrically-placed positive charges @ in the 
distribution, each a distance r from the central point, 
where r < R. According to Coulomb’s law and the 
principle of superposition, the magnitude of the field 
at P due to these two charges is 


Q 1 Q 1 
Bi = Arreg (R— r)? a Areg (R+ 1)? 
Q (R+r)2+(R-r)? 
~ dren (R—1)2(R+r)? 
_ 2Q . R? +17? 
Aneg = (R? — r?)?” 
For comparison, if the two charges were moved to the 


centre of the cylinder, the magnitude of the field at P 
would be 


2Q 1 
2= KS 
ATED ie 
Now, it is easy to see that 
Re +r? . Rr 1 2 
(R2 = r2)2 (R2 = r2)2 _ (R?2 = r?) R2 


so EF; > Ey. Since this is true for all pairs of charges 
symmetrically-placed on either side of the centre of 
the cylinder, collapsing the entire distribution down to 
the central point would reduce the magnitude of the 
field. 


This calculation does not use Gauss’s law, but it 
shows that the results obtained from Gauss’s law for 
spheres and cylinders that are collapsed radially are 
rather special. 


Figure 9 For Solution 2.6. 


Solution 2.7 No, the statement is wrong. Obviously 
the electric field at any point depends in principle on 
all charges, no matter where they might be. Although 
charges outside the spherical surface do not contribute 
to the total flux over the surface, they do influence the 
electric field at individual points on the surface. It is 
just that they contribute a positive flux in some 
patches and a negative flux in others, and all these 
contributions cancel out when we integrate over the 
whole surface. 


Solution 2.8 No, this is not possible. For example, it 
will not help to surround the charged cube by a cubic 
Gaussian surface because the electric field produced 
by the charged cube is not perpendicular to the faces 
of this Gaussian surface and is not constant over them. 
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No other simple Gaussian surface is any better. 

We generally need spherical, cylindrical or planar 
symmetry to have any chance of obtaining complete 
knowledge of the electric field from the integral 
version of Gauss’s law. 


Solution 2.9 (a) The thickness of the membrane is 
much less than the radius of the cell, so we model 
it by a pair of plane sheets of charge. Inside the 
membrane, we can use the formula for the field 
between the plates of a parallel plate capacitor: 
p= 2 
E0 
25210 ° Cm" 
~ 8.85 x 10-12C?N-! m=? 
=28xi0 Nc, 


(b) The membrane provides two spherical shells of 
charge. A uniformly-charged spherical shell produces 
no electric field in the spherical volume it surrounds, 
so there is no electric field in the cell. 


(c) Exploiting spherical symmetry we choose a 
Gaussian surface concentric with the cell and of 
greater radius than the cell. This surface contains zero 
total charge so Gauss’s law tells us that the electric 
flux is zero over the surface. By spherical symmetry, 
the electric field is everywhere radial and has the same 
magnitude at all points on the spherical surface. We 
conclude that the electric field vanishes outside the 
cell. 


Solution 2.10 (a) By symmetry, the electric field is 
aligned with the z-axis, E = E(z)e,. For z > 0 it is 
in the positive z-direction and for z < 0, it is in the 
negative z-direction. 


To determine the electric field outside the slab, we 
choose a cylindrical Gaussian surface of cross-section 
AA with end-faces at z and —z outside the slab. Then 
Gauss’s law gives 


E(z) AA+ E(-z) AA= se 
0 
where F is the magnitude of the electric field. 
By symmetry, F(z) = E(—z), so 
_ pa 
=o 
So the electric field is 


E(2) 


NIQ 


To determine the electric field inside the slab we 
choose a similar cylindrical Gaussian surface of 
cross-section A.A, but with end-faces at z and —z 
inside the slab. This time Gauss’s law gives 


pA x 2z 


E(z) AA+ E(-z) AA= = 


Following through a similar calculation to that given 
above, we conclude that 

pz d d 

E= — f —= =. 

2 e, for 5 <z< 5 
(b) In both regions outside the slab, the electric field is 
constant and therefore has zero divergence. This 
agrees with the differential version of Gauss’s law 
because there is no charge density outside the slab. 
Evaluating the divergence inside the slab we obtain 


OE, p 


div EK = = 
ai Oz E90” 


which is consistent with the differential version of 
Gauss’s law. The surfaces of the slab create no 
problems because the electric field is continuous there. 


Solution 2.11 We solve for the plane and the slab 
separately and then use the principle of superposition. 
Using Gauss’s law we can find the electric fields due 
to the plane and slab everywhere. The calculations are 
not given here (they essentially repeat the working of 
the previous problem) but the results are as follows. 


The field due to the slab is 


pd 

ieee fe d/2 

a ‘ar 2 > d/2, 

pilab - e for —d/2 < z<d/2, 

0 

pd 

eae f —d/2. 
Jeo or z < —d/ 


The field due to the plane is 


—— for z > d/2, 
plane 2 0 
EE = a 
—_— for z < d/2. 
2€9 ° / 


pd +o d 
for z > 5, 
2€9 . A 
202 — 
B, = Es 4 pplene — for -f<z<4, 
0 
pd + d 
- for z << —§. 
2€9 ° ? 


This vanishes inside the slab at z = 0 /2p, provided 
that —d/2 < 0/2p < d/2. 


Solution 2.12 The negative charge within the Bohr 
radius is 


uo 9 4e [® 4 
q= | px4arr ar=—% [ r* exp(—2r/ao) dr. 
0 a JO 


Changing the variable of integration to x = 2r/ag we 
obtain 


_ 4e ag\? f? 5 
Q=-3x (9) [« exp(—2) dx 


_ a x 0.647 = —0.32e. 


Allowing for the positive point charge at the centre, 
the net charge within the Bohr radius is 0.68e. 


Let the radial component of the electric field at ap be 
E,.. Then Gauss’s law gives 


trea 2i= pene 
E0 
so 
_ 0.68e | 
_ Ameya? 


068 X16 10-"C 
Ar x 8.85 x 10-12 C?N7! m-? x (5.29 x 10-!! m)? 
=35x1l0UNC!. 


Solution 2.13 It is helpful to add a second 
hemispherical shell, to create a complete 
uniformly-charged spherical shell. Because the charge 
distribution is spherically symmetric, the electric field 
inside the complete sphere is zero. Let the electric 
field due to the first hemispherical shell be E, and the 
electric field due to the second hemispherical shell be 
Ez. Then, inside the sphere we have 


E, + Ey = 0. («) 


The second hemispherical shell is obtained by 
reflecting the first hemispherical shell across the 
imaginary surface S. It follows that the field pattern 
produced by the second hemispherical shell is the 
reflection of the field pattern produced by the first 
hemispherical shell. We have 


E\y = Eon and Eqy = Ey while F\, = —E», 


Combining this with Equation (*), we see that the x- 
and y-components vanish. It follows that the electric 
field due to a single hemisphere is perpendicular to S. 


Solution 2.14 (a) From Gauss’s law, the total 
electric flux over the cube is Q/e9. Because the 
charge is at the centre of the cube, each of its 6 faces 
contributes equally to the flux. The flux over a single 
face is Q/6¢e0. 


(b) Faces of the cube can be split into two groups: 
three “‘adjacent” faces (adjacent to the corner where 
the charge is located) and three “opposite” faces 
(opposite to the corner where the charge is located). 
By symmetry, the three “adjacent” faces have identical 
fluxes and the three “opposite” faces have identical 
fluxes. (See Figure 10). 


Imagine now 8 cubes, identical to the one considered, 
arranged to form a larger cube, with double the side 
length, in such a way that the point charge Q is near 
the centre of the larger cube. The flux over the entire 
larger cube is Q/eo and the flux over each of its faces 
is Q/6€0. One of the faces of the larger cube can be 
split into four equal squares with equal fluxes, one of 
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which is an “opposite” face from the original cube. 
We conclude that the flux over each “opposite” face is 


i ee 
4 6eo Q4eq’ 
leaving each “adjacent” face to have a flux 
1/Q  3Q)\_ 7Q 
3 \eg deg} deg’ 


Figure 10 For Solution 2.14b. 


Solution 2.15 Because the shells are conducting, the 
charge spreads out evenly over each shell, ensuring 

a spherically symmetric charge distribution and 
therefore a spherically symmetric field. The electric 
field outside any spherically symmetric charge 
distribution is as if all the charge were concentrated at 
the central point. Using this fact, the field outside both 
spheres points radially inwards when q, + q2 < 0, is 
zero when q; + g2 = 0 and points radially outwards 
when q; + q2 > 0. 


Solution 2.16 Suppose we wish to trap a positive 
ion. In order for a restoring force to act on this ion the 
electric field must point inwards at every point on the 
boundary of the trapping region. So the electric flux 
over the surface of the trapping region is negative. 
Gauss’s law then requires there to be a negative charge 
in the trapping region. This contradicts our initial 
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assumption that the trapping region is a vacuum 
containing no charges, apart from the ion. The 
positive charge of the ion does not help here because 
(a) it has the wrong sign and (b) it is irrelevant because 
we are looking for the origins of a field that exerts a 
force on the ion. The ion cannot exert a force on itself. 


The fact that stable equilibrium cannot be achieved 
with externally-imposed static electric fields alone 

is known as Earnshaw’s theorem. Its practical 
consequence is that an ion trap cannot be constructed 
using electrostatics alone; successful ion traps use 
magnetic fields. 


Solution 2.17 According to Gauss’s law 


OE, 2 Oy 2 OE, p 
Ox Oy Oz & 


The partial derivatives OF, /Oy and OF, /0z are equal 
to zero because the field has no y- or z-components. 
The partial derivative OE, /Ox can be written as an 
ordinary derivative because EF, depends only on z. 
Hence Gauss’s law becomes 


dE, 
dz 


= e cos(ka). (See Figure 11). 
0 


Figure I] For Solution 2.17. 


Integrating both sides with respect to x gives 


Ee = a sin(kx) + C, 
where C is an arbitrary constant. Using the initial 
value E,, = 0 at x = 0 we see that C' = 0, so we 
conclude that 


E= aE sin(kx) ex. 


Book 1 Chapter 3 


Solution 3.1 A current J flowing in a loop of 
cross-sectional area A produces a magnetic moment 
of magnitude m = |I|A. Thus 

m 8.0 x 10? A m? 


[f|= = 
7 


=46x10°A. 
R27 x (2.35 x 10°m)2 . 


Solution 3.2 The current flow is in the direction of 
J, which is in the direction of A. If the current 
density is to produce no accumulation of charge, its 


divergence must vanish. We have 


Ox 
= —k, A, exp(—k-r), 


with similar results for the y- and z-components. So, 
div J = —(k- A) exp(—k-r). 


The required condition for div J to vanish is therefore 
that k is orthogonal to A. 


Solution 3.3 Consider a current loop in the plane 
of the page, with the current in an anticlockwise 
direction viewed from above. The Biot—Savart law 
gives the magnetic force on current element J 6], due 
to current element J dl) and can be written as 


_ bo Tol, X (I dla X T12) 


Fio % 
Ar 5 


One application of the right-hand rule shows that 

I dlg X 12 points out of the page towards you. This is 
true for any current element J 6], that does not 
coincide with J 51,. A second application of the 
right-hand rule shows that J dl, X (J dlg X Ty2) is in 
the plane of the loop and points radially outwards 
from its centre. Adding up contributions from all 
current elements J 612 we conclude that the total 
magnetic force on J 5], is in the plane of the loop and 
points radially outwards. 


Solution 3.4 This question refers back to Gauss’s 
law. Consider a Gaussian surface surrounding the 
empty trapping region. There is no charge inside it, 
so the total electric flux over the surface is zero. 
Suppose the electric field points towards the origin 
along the z-axis. Then parts of the Gaussian surface 
near the z-axis contribute negative electric flux. To 
compensate for this, other parts of the Gaussian 
surface must contribute positive electric flux; this 
implies that the electric field points outwards, away 
from the origin, in at least one direction perpendicular 
to the z-axis. So, if the electric field confines along 
the axis of the magnetic field, it opposes confinement 
in a direction perpendicular to the field. In practice, 
this means that the electric field must not be too large. 


Book 1 Chapter 4 


Solution 4.1 The spherical symmetry suggests the 
use of spherical coordinates. Spherical symmetry 
shows that the spherical components of the magnetic 
field do not depend on the angular coordinates @ and 
o. So, 


B(r) = B,(r) e, + Bg(r) eg + Bo(r) eg. 


It is possible to find 180° rotations through the centre 
of the current distribution that leave the current 


distribution unchanged but reverse By and Bg. 
Symmetry then requires that 


B,(r) er + Bg(r) eg + Bo(r) eg 
= B,(r) e, — Bg(r) eg — Bo(r) €6, 


so the angular components vanish, leaving a radial 
magnetic field 


Bit) = Byer) ee 


The no-monopole equation tells us that the magnetic 
flux vanishes over any closed surface. Evaluating the 
magnetic flux over a spherical surface S of radius R, 
centred on the origin, gives 


| B- dS = B(R) x 4nR?, 
S 


so B(R) = 0 for all R > 0. The magnetic field also 
vanishes at the origin because spherical symmetry 
prevents it from pointing in one direction rather than 
another. 


Solution 4.2 It is sensible to express the field in 
cylindrical coordinates, with the z-axis perpendicular 
to the ring, and passing through its centre. Axial 
symmetry ensures that the cylindrical coordinates of 
the magnetic field do not depend on ¢, so 


Bir) = B, (7,2) & + Bg(7, 2) eg + Be(r, 2) ex: 


Consider time reversal, followed by a reflection in a 
plane that contains the z-axis. The net effect of these 
two operations is to leave the current distribution 
unchanged. Time reversal reverses all components of 
the magnetic field. Using the reflection rule, the 
reflection reverses B,(r,z) and Bz(r, z) at any point 
in the plane of reflection (since these components are 
parallel to the plane), but By(r, z) is unchanged (since 
this component is perpendicular to the plane). So the 
net effect of both operations is to convert By(r, z) to 
—Bg(r, z). Hence the ¢-component of the magnetic 
field vanishes and the magnetic field therefore takes 
the form 


Bir) = B,(r, z) er, + B,(r, z) ez. 


Note that a reflection in the plane of the ring is useless 
in this problem because it changes z to —z. This 
prevents us from comparing two expressions for the 
magnetic field at a fixed point. 


Solution 4.3 Because the situation is axially 
symmetric, we use cylindrical polar coordinates and 
express the magnetic field as 


B= Bg(r)es, 
where By(r) is a function that remains to be 
determined. 
Let R be the radius of the coaxial cable. Consider a 
circular path C' of radius r > R, centred on, and 


perpendicular to, the central axis of the cable. The line 
integral of the magnetic field around this path is 


Bg(r) x 2rr. The circle C is the perimeter of a disk 
that covers the whole cross-section of the cable. 
Because the current carried by the inner cylinder is 
equal in magnitude and opposite in direction to the 
current carried by the outer tube, the total current 
flowing through this disk is zero. Consequently, the 
integral form of Ampére’s law gives 


Bg(r) x 2nr = 0 forallr > R. 


We conclude that the magnetic field vanishes 
everywhere outside the coaxial cable. 


Solution 4.4 We can again use symmetry to argue 
that the magnetic field has no z- or z-components. We 
choose now a rectangular closed path in the plane 

x = 0, defined by the points A at (x, y, z) = (0,0, d), 
B at (0, L,d), C at (0, L,0), and D at (0,0, 0). 
Exercises 4.5 and 4.10 in Book 1 showed that the 
magnetic field of an infinite current sheet does not fall 
off with distance from the sheet. Using this fact, 
together with the right-hand rule, we conclude that 
there is no magnetic field along AB — contributions 
from the two current sheets cancel out. Let the 
magnitude of the magnetic field along CD be B. Then 
Ampeére’s law gives 


N 
Bis Hol xo Ty 


so, between the current sheets, the magnetic field has 
magnitude 


Hol N 
l 
Not surprisingly, this is twice as large as the magnetic 

field near a single current sheet. Indeed, it would be 


fine to use the principle of superposition, and the 
right-hand rule, to arrive directly at this answer. 


RS 


Book 1 Chapter 5 


Solution 5.1 Imagine building up the sphere in a 
series of spherical shells. Suppose we already have a 
sphere of uniform charge density p and radius r and 
we add a spherical shell of uniform charge density p 
and thickness dr. The thickness dr is taken to be small 
enough for us to neglect repulsion of elements within 
the spherical shell. The sphere has total charge 
Arr p/3 and the shell has total charge 47? dr p. The 
electric field of the sphere is just as if all its charge 
were concentrated at its centre, so the potential energy 
between these two charge distributions is 

1 Anrp 


6U = x 


= Anr? or p. 
Amer 3 i 


The total potential energy is found by integrating from 
0 to the full radius R of the sphere: 


2 R 
U= ae x ate) | r’ dr. 
0 


ATE 3 
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The total charge of the sphere is Q = 47 R°p/3 so 


2 
(a n 
5 4t7egR 


Solution 5.2 At the centre of the sphere, a thin 
spherical shell of radius 7, thickness dr and charge 
density p contributes a potential 


1 
V= Arr? dr. 
ieee Tr Or 


The total potential at the centre due to the whole 
sphere is 


R 2 

R 
v=2f rdr =P : 
Eo Jo 2€0 


Since Q = 47 R®p/3, we conclude that 


5 
a 2 AregR 


An alternative method can be used. We can explicitly 
calculate the work done transporting a charge q¢ 
radially from infinity to the centre of the sphere. 
Outside the sphere, the electric field is 


Q 


Mout = Teor? © 


so the work done bringing the charge from infinity to 
the surface of the sphere is 


qQ se qQ 


Anég Jong T2323» 4meQR’ 


Wout — 


The electric field inside the sphere was derived in 
Worked Example 2.2 (Book 1, pp. 48-9), and is given 
by 

Qr 
——.e 
An eo R3 J 
so the work done bringing the charge from the surface 
of the sphere to the centre of the sphere is 


qa f° qQ 


Are oR? Jp . AregR 


Ein = 


1 
Win = - eae 


The total work done is therefore 


3 
W = Wout + Win = 9 fp 
and the potential is 
WwW 3 @Q 
qs 2 AnegR’ 


Note that there is no simple relationship between the 
potential energy of the whole sphere, which was 
calculated in Exercise 5.1, and the potential at the 
centre of the sphere. In general, potential energy and 
potential are very different concepts. 


Solution 5.3. At the centre of the disk, a thin ring of 
radius 1, thickness dr and surface charge density 7 
contributes a potential 


6V= o X 2rr or. 


1 
4neor 
The total potential at the centre due to the whole disk 
is 
of ap = oR 


Since Q = 7R?<, we conclude that 


__@ 
QreoR 


Solution 5.4 The initial electrostatic energy is 


2 iF 
i=. 
‘5 AneoR 
(see Exercise 5.1). The final electrostatic energy is 
3 2? 
Uz = 2x (QR) _ 92/8, = 0.63}. 


5 Ame (R/21/3) 
The change in electrostatic energy is 


Ug — U, = —0.37 04. 


Solution 5.5 The cloud and earth can be modelled 
as a large parallel plate capacitor. The electrostatic 
energy stored in the capacitor is 


1 1 
Ce | ie 5 30C x 1.0 10° V =1.5x 101°). 


This potential energy goes into light, sound, heat and 
the potential energy of ionized atoms. 


Solution 5.6 In equilibrium, components of the 
electric field parallel to the surface of the conductor 
must vanish. If they did not vanish, free charges would 
flow in the conductor, contradicting our assumption of 
equilibrium. Since the electric field has no 
components parallel to the surface of the conductor, it 
must be perpendicular to the surface of the conductor. 


Alternatively, we can argue that, in equilibrium, the 
surface of the conductor is an equipotential surface. 
The electric field is always perpendicular to an 
equipotential surface. 


Book 1 Chapter 6 


Solution 6.1 This problem can be tackled using 
either the magnetic force law or Faraday’s law. We 
will consider both methods in turn. 


Using the magnetic force law. A charge q in side AB 
experiences a magnetic field B,(x) and a magnetic 
force with y-component F,, = —qu,B.(x). As it 
moves from A to B, the charge is displaced by 


Ay = —L and the work done by the magnetic field is 
F, Ay = quzB,(x) x Ly. 


Similarly, a charge g in side CD experiences a 
magnetic field B,(a + L2) and a magnetic force 

with y-component F;,, = —qu,B-(x + L2). As it 
moves from C to D, the charge is displaced by 

Ay = +1 and the work done by the magnetic field is 
Fy, Ay = —qv,B,(x + L2) x Th. 


These are the only electromagnetic forces that do 
work on the charge as it moves round the loop ABCD, 
so the total electromagnetic work done is 


W = qvu-l;(B,(x) — Bz(x + Le)| 


and the emf is 


Vont = — = tei [B.() ~ Be( + La) 
Using Faraday’s law. As the loop moves, the flux 
through it changes because the leading edge, CD, 
gathers more flux while the trailing edge, AB, leaves 
flux behind. If the loop moves a small distance Ax, 
the flux change associated with displacement of the 
leading edge is B,(x + L2) x L, Ax while the flux 
change associated with displacement of the trailing 
edge is —B,(x) x L; Az. Ina small time interval At, 
the loop moves a distance Ax = v, At and the total 
change in flux is 


A® = [B,(a@ + Lo) — B,(x)| x Livy At. 


The rate of change of flux is 


oe = veli[B.(n + La) ~ B-(2) 
and Faraday’s law gives 

Vane =~ = veLa[Be(a) ~ Bee + La)] 
as before. 


Solution 6.2 No emf is induced. The magnetic flux 
through the loop is constant, so Faraday’s law tells us 
that the induced emf is zero. 


Solution 6.3 Consider a positive charge q in the rod, 
at a distance r from the pivoted end. This charge has 
speed v = rw and since its motion is perpendicular to 
the magnetic field it experiences a magnetic force of 
magnitude 


F=qB=qrwB. 


Using the right-hand rule shows that this force acts 
radially outwards along the rod. If a charge q is 
transported from the pivot to the free end of the rod, 
the work done by the magnetic force is 


L L 
W= [ Fdr= we | rdr= squ BL’, 
0 0 


and the voltage drop between the two ends of the rod 
is 


WwW 
View eo 5 BwL?. 
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This is the voltage drop as we move from the pivot to 
the free end of the rod. 


Solution 6.4 As the rod rotates it sweeps out a 
circular sector OAB, where O is the pivot, OA is the 
initial position of the rod and OB is its position at a 
later time. 


If @ is the angle between OA and OB, so that d= Ww, 
the area of the sector OAB is 


Y 2_ 1972 
A= on x wL* = 5 OL". 
The magnetic flux passing through the sector is 
® = 5OLB. 
and Faraday’s law gives 
d® 
dt 


This is the emf around the circuit OAB. The entire 
voltage drop takes place within the moving rod, this is 
the change in voltage as we move from the free end of 
the rod to the pivot, so it is minus the voltage drop 
calculated in the previous exercise. 


Book 1 Chapter 7 


Vemt = ——- = —gwL?B. 


Solution 7.1 Combining the equation of continuity 
with J = oF, we obtain 
Op 
Ot 
The conductivity o is a constant, so it can be moved 
outside the divergence. We can then use Gauss’s law, 
div E = p/p, to obtain 


—div J = —div(oE). 


Op or 

SS a p 

ot E0 
At a fixed point in space, this differential equation has 
the solution 


p= poe! = poe", 


where fo is a constant (the charge density at t = 0) 
and T = €q/a is the characteristic time constant of the 
exponential decay. 


Comment: To get a feeling for the timescales 
involved, let’s insert a typical value for the 
conductivity of distilled water (~ 107 0 mm) in 
our expression for 7. This gives 

8.85 x 10°? C?N-' m? 

— 
10 0-*m-* 

So any excess charge is dispersed very rapidly. In 


metals, charge dispersal is even faster, and for most 
purposes can be taken to be instantaneous. 


~ 1077s. 


Solution 7.2 Gauss’s law and the no-monopole 
equation are satisfied because 


div E = La, =, 
Ox 


: 10 
divB = poy =), 
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Faraday’s law gives 


o=cwin+ 
_ of 1 Of 
“Ae 0 aot ® 


_ fot , ef 
= (Se cae) 


and the Ampére—Maxwell law gives 


OE 

O Scwl B <piig — 

cur Eo Lo Ot 

7 Lf, 1 oF 
~ cdz” eat” 


Loy , Lor 
~ & on |e 
Both these equations are satisfied because, letting 
u = z—ct, we have 
Of oOfodu Of 
dz Oudz du 
and 


Of  dfdu af 


at Oude du 


Hence 
of __10f 
Oz ~=—c OF 


Comment: These electric and magnetic fields remain 
undisturbed in shape and travel in the positive 
z-direction at the constant speed c. This shows that 
electromagnetic disturbances need not be sinusoidal. 
In fact, any shape travelling with constant speed 

c can be represented as a linear combination of 
sinusoidal waves (using a mathematical technique 
called Fourier analysis). The monochromatic plane 
waves considered in the main text can therefore be 
thought of as basic building blocks from which all 
other electromagnetic waves can be made. 


Solution 7.3 The electromagnetic wave is a plane 
wave, so the electric field strength is constant in any 
plane perpendicular to the direction of propagation 
and the electric field strengths at P, Q and R are 
identical. If you gave a different answer, you may 
have misinterpreted Figure 7.10. This indicates the 
z-dependence of the electric and magnetic fields but 
does not (and cannot) indicate dependence of the 
fields on x or y. The lack of x- or y-dependence 
follows from the plane-wave nature of the wave. 


Solution 7.4 (a) The time-varying electric field in an 
electromagnetic wave drives a current to and fro along 
the length of a straight aerial. The aerial is most 
effective when it is oriented parallel to the line of 
oscillation of the electric field. 


(b) The time-varying magnetic field in an 
electromagnetic wave induces an emf around a 
circular-loop aerial. The aerial is most effective when 
the plane of the loop is perpendicular to the line of 
oscillation of the magnetic field. The plane of the loop 
then contains the line of oscillation of the electric field 
and the direction of propagation of the wave. 


Solutions for Book 2 


Book 2 Chapter 1 


Solution 1.1 (a) There is spherical symmetry (and 
so we will use spherical coordinates). 


(b) Using the spherical symmetry, we can conclude 
that 


E(r) = E,(r)er. 


First consider the case r < R. We will use the integral 
form of Gauss’s law, 


[w-as== | pay, 
S £0 JV 


with a spherical Gaussian surface of radius r. For this 
surface E is everywhere parallel to dS and has 
constant magnitude E,.(r). Thus 


| E-dS = 4nr7E,(r). 
S 
The total charge enclosed by the Gaussian surface is 


[ow-=| ksdns?ds = 4nk | sds 
V 0 0 


= Ank [2 5] 


_ 4 
s<0 = tkr*. 


So, by Gauss’s law, 


We now consider the case r > R. As above, fora 
spherical Gaussian surface, 


| E-dS = 4nr7E,(r). 
S 


The total charge enclosed by the Gaussian surface is 


[ow-=| p(s) 47s? ds 
Vv 0 
R 


= p(s) 4s? ds, 


0 
since p(r) = 0 forr > R. So 


R 
[ eav =ank f sds = rkR’. 
V 0 


Hence, by Gauss’s law, 


1 
Anr?E,(r) = — wkR', 


E0 
i.e. 
k Rt 
E,(r) = 4é0 ae (r > Rs 


Combining the two cases, 
E = E,(r)e, = pe 


(c) As E(r) = E,(r)e,, 


So for r < R, 


ldf/k4 
dwE= 5 5 (or) 


k 1 
— a — - p(r). 
Similarly, for r > R, 


1d (ke 
weal o (eR 
” r? dr () 


1 

S p(r) 
Hence, the differential form of Gauss’s law is satisfied 
both inside and outside the charge distribution. 


Solution 1.2 If the vector field B(r) represents a 
magnetic field then it must satisfy the no-monopole 
law, 


divB = 0. 
For the vector field F, 


0 0 0 
div F = 5g) + ay y) + 5 (0) 


z 
=1—-1+0=0. 

Hence F can represent a magnetic field. 

Using Ampeére’s law, the associated current density is 


1 1 
J = —curlB = — curlF 
Ho Ho 


ey ey e, 


x 22—-y 0 
2 
= —e,. 
Lo 
For the vector field G, 
O O 0 
divG = — (-2 — (4 5 — (0 
VG = 5 (-22) + 5 (y+ 8) + 5 (0) 
=-24+44+0=2. 
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Hence G cannot represent a magnetic field because 
divG £ 0. 


Solution 1.3 


er ey e& 
curl E = 2 5 “ 
Epsin EC 2 )] 0 0 


E 
— au ics EC _ al e,. 
Cc Cc 


Using the differential form of Faraday’s law, 


B 
gE = —curlE 


ot 


Equ y 
= —— cos EC _ al e,. 
Cc 6 


Integrating with respect to t, we obtain 


B(r,t) = -= sin EG - *)| e, + Bo(r), 


where Bo(r) is an arbitrary field which is independent 
of time. However this field must satisfy the 
no-monopole law, so div Bo = 0. 


Solution 1.4 From Maxwell’s equations, in empty 
space, with no charges and currents, 


B 
curlE = 2 


Ot 


and 


OE 

IB= a 

cur Eo Lo aL 

Taking the curl of the first equation and then using the 
second equation, we have 


curl (curl E) = -° (curl B) 


O ( en _ 0°E 
aL EoHo aL’ EoHo ap 


But 
curl (curl E) = grad (div E) — V7E. 

Also in empty space, by Gauss’s law, 
div E = 0. 

Hence 

1 OE 

c oF” 

as required, where 


1 
Voto 


Solution 1.5 The force on the electron is given by 
the Lorentz force law 


F = q(E+ v xX B), 
26 


VE = 


where g = —e. Now 
E = Ee, 
and 
v X B= (v,e2+ uyey + 022) X (Bez) 
= vy Be, — vzBey. 
So the force acting on the electron is 
F = —(ef + eBu,)ez + eBuze, 
and the equations of motion are 
Mt, = —(eH + eBy,), 
My = eBrz 
and mv, =0. 


Comment: The motion is a combination of a circular 
motion with angular frequency w = eB/m in the 
xy-plane, a uniform motion in the negative y-direction 
with speed F/B and a uniform motion in the 
z-direction. This type of trajectory was used in the 
original determination of e/m for photoelectrons. 


Book 2 Chapter 2 


Solution 2.1 
(a) The total number of atoms in the sphere is 
N= (2.17 « 10" mi") 
x ($m x (1.50 x 107%)? m*) 
=3.07 x10", 


So the magnitude of the average dipole moment of an 
argon atom is 
I 
(p) = N 
_ 5.50 x 107-'*Cm 
3.07 x 1020 
=1.79x 10°" Cm. 


The direction of (p) is the same as that of IT. 


The separation of the centres of charge in this dipole is 
_@) _ ©) 
q Ze 
.. 199% 10-* Cm 
18 x (1.60 x 10-!9C) 
= 6.23 x 107!7 m. 


(You might like to compare this with the diameter of 
the argon atom, which is 3.4 x 10~!9 m. The charge 
separation is very small compared to the atomic 
diameter.) 


(b) The polarization P of the sphere is in the same 
direction as (p) and has magnitude 


P= n{p) 
= (2.17 x 1078 m~3) x (1.79 x 107-34 Cm) 
=5.89x 10° Cm. 


Solution 2.2 (a) The bound volume charge density is 
Pp = —divP. 
Using spherical coordinates 
1 9 1 oO 
72 Op” PH) + Tain 8 | 
1 OPs 
rsin@d 0 
1 oO 
= 5, (ar cos? 0) 
lL «2 ( 
rsin@ 06 
= 3a cos? § — —— (— sin? 6 + 2cos” 6 sin @) 
sin 
=3acos? @ + asin? @ — 2acos? 6 


=a. 


div P = sin 0 Po) 


ar cos 6 sin? 0) 


So pp=-—a. 
(b) The bound surface charge density is 

op = P-n. 
The outward-pointing unit normal to the surface of the 
sphere is Nn = e,., and on the surface of the sphere 
r= R. So 

op = (aR cos” 6 e, — aRcos 6 sin 0 eg) - e, 


= aRcos’ 6. 


Solution 2.3. By symmetry, the electric field E and 
electric displacement D inside the dielectric material 
are in the same direction as the external electric field. 
Also in an LIH material, D = ceoE. 


Furthermore, at the boundary between the dielectric 
material and air, D, is continuous. As the 
boundary is perpendicular to the direction of the 
electric displacement, the magnitude of the electric 
displacement in the dielectric slab has magnitude 
Dgiab = Dair 
= €0 Fair 
= (8.85 x 1071? x 12000) Cm~? 
=(1 105’ Cm-. 
So the electric field in the slab has magnitude 
Esiab = Dagab/EE0 
= ales © 
= 12000/6Vm~! 
=20s 10° Vm, *. 
The polarization in the slab has magnitude 
Paab = XBE0Eslab 
= (€ — DeoLEstap 
= (5.0 x 8.85 x 107? x 2000) Cm~? 
=69% 10 °Cm, 


Solution 2.4 (a) Neglecting edge effects, the field 
lines of D (and E) will be perpendicular to the plates. 


We apply Gauss’s law 


[v-as= f nav 
S V 


to a pillbox with one end in plate | and the other end 
in material 1. 


Within the conducting plate, D is zero, and so Gauss’s 
law leads to 


Di=Co 


at all positions in material 1. By the continuity of D, 
at the interface between the dielectrics, 


Dy = Di =o. 


Using the integral form of Gauss’s law for a similar 
pillbox with one end in plate 2, we find that the charge 
density at plate 2 is —c. (Alternatively we could have 
used a pillbox with one end in plate 1 and the other 
end in plate 2 and concluded that the total charge in 
this pillbox is zero.) 


Now 
D = ceoE. 
So 
1 1 
Fy = — Di = — 0 
E1E0 E1EQ 
1 1 
and Ey = —— Dy = —— 
E2EO EQEO 
Also the polarization 
P = xne0E = (€ — 1L)eoE. 
So 
—1 —1 
i a o0 and P= = Oo. 
Ej EQ 


(b) The potential difference across dielectric 1 is 


d 
db, = mel Sy 
E1E0 


and across dielectric 2 is 
2 
dg Ey = ———" 10% 
E 


So the total potential difference across the capacitor is 
1 /d d 
V=— (2 + =) 0. 
Fo (El 2 


Now the capacitance is 


Q oA 
GC = — = —— 
V Vv’ 
so that 
C= eo A 
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Solution 3.1 We will use cylindrical coordinates, 
with the z-axis along the axis of the wire and in 
the direction of the electric current. The magnetic 
intensity H and the magnetic field B are azimuthal, 
i.e. in the direction of eg. 


Using Ampeére’s law, 


gua / 3,-as, 
Cc S 


for a circular path of radius r centred on the axis of 
symmetry of the wire, 


QarH = nr? J, 
which gives 
H = 5rJ. 
Now B = plioH and so 
B= 5 Lplor J. 
Solution 3.2 In the absence of free surface currents, 
the components B, and H) are continuous. 


For LIH materials, 
H = B/pyo, so Hisin thesame direction as B. 


The continuity of B, at the interface gives 
By, cos 60° = Bo cos 8, 


whereas the continuity of Hj gives 
1 1 
—— B, sin 60° = ——Bosiné. 
LiLo H2Ho 
Dividing these last two equations leads to 
1 . 1 
— tan60° = — tan. 
1 M2 
So 


3.0 
tan = »? tan60° = “V3. 
1 1.5 


Hence 

@=tean*(2V/3) = 74°. 
Solution 3.3 (a) The magnetic field lines are circular 
lines around the toroidal ring. 


We apply Ampére’s law, 


gua / 3,-as, 
C Ss 


to a circular path C of radius r within the ring (as 
shown in Figure 3.20 in Book 2). 


If there are N turns carrying a current J, the total 
current flowing through the disc S is NI. So 
Ampére’s law gives 

27rrH = NI. 
28 


Hence 
_ NE 
- Qnr 
2000 x 1.5 4 
= ——_——Am 
2m x 0.12 
= 3980Am! = 4.0 x 10°?A m' to2 sig. figs. 
The magnitude of the magnetic field is 
B= ppoll 
— (400 x 4n x 10~" x 3980) T 
= 2.0T. 
(b) At the interface between the air in the gap and 


the iron alloy, the perpendicular component of the 
magnetic field B is continuous. So 


Bair = Biron: 
It follows that 
Hairy = LHizon- 


We again use Ampére’s law for the same circular path 
Cas used in part (a). In this case we obtain 
(2ar _ d) Hiron + dHay = NI, 
where d is the width of the air-gap. 
So 
NI 
(2rr — d) + pd 
_ NI 
~ Qar+(u—1)d 
2000 x 1. 
_ x 1.5 Agel 
2m X 0.12 + 399 x 0.003 


= 1538Am*}=15x108Am 
to 2 sig. figs. 


Hy; iron — 


Now 
Alaiy = LH iron 
= 400 x 1538 Am7! 


= 6.152 x 10°Am7! =6.2 x 10°Am7! 
to 2 sig. figs. 


Using B = poH for either the iron ring or the air in 
the gap, we obtain 
Biron = Bair = 0.773 T = 0.77 T to 2 sig. figs. 
Solution 3.4 The bound current density is given by 
J, = curlM, 


where M(r) is the magnetization. In a linear, 
isotropic, homogeneous material 


M = *2.B = pyeH, 
Lo 


where pi and yz are constants. Hence 


Jp = uxe curl H. 


But, from Ampére’s law 


curlH = Jr, 
and so 
Jb = uxeBIe. 


In other words, the bound and free current densities 
are proportional to one another. In particular when the 
material carries no free current, the magnetization 
current is zero at all points in the volume of the 
material, leaving only the possibility of bound surface 
currents ij, = M X n. 


Comments: Inhomogeneities may exist at, or near, the 
surface of a material, so that bound current density 
may exist in the surface layer. Furthermore, non-linear 
materials (such as ferromagnets) can certainly have 
non-zero bound current densities in the absence of 
free currents. 


Book 2 Chapter 4 


Solution 4.1 (a) This form of the differential 
version of Gauss’s law would be used in regions 
where the total charge density p is known. 


(b) This form of the differential version of Gauss’s 
law would be used in regions which contain dielectric 
materials, and where the free charge density pr is 
known. 


(c) Gauss’s law in its integral version would be used 
in regions where the total charge density p is known, 
and where there is a high degree of symmetry. 


(d) This form of the integral version of Gauss’s law 
would be used in regions which contain dielectric 
materials, where the free charge density p¢ is known, 
and where there is a high degree of symmetry. 


(e) This version of Poisson’s equation would be 
used in regions where the total charge density p is 
known, as in part (a). 


(f) This version of Poisson’s equation would be used 
in regions which contain an LIH dielectric material 
and where the free charge density is known. 


(g) Laplace’s equation would be used when the total 
charge density p(r) = 0 everywhere in the region. 


Solution 4.2 The potential V(x) in the dielectric 
material satisfies Poisson’s equation 


1 
V2V =—-— pr 
EE 


and because of the symmetry this reduces to 


d?v 1 
—~ =—-—(A+B : 
da = (A+ Ba) (O<e¢< 8) 


Integrating, we obtain 


dV 1 
—_=__ (4 1By? 
a seg (Att 2 x) +C 


1 


Vz) = EE 


(Ax? ~— 4 Bu”) +Ca+D, 


where Cand D are constants of integration. 


We apply the boundary conditions at x = 0 and x = s 
to determine C and D. 


At c = 0, V(0) = 0 and so D=0. 
At x = s, V(s) = Vs and so 
1 


_ 14.2, 1p.3 
Vv; = a (5As + ¢Bs°)+Cs. 


Hence 


Vs val C2 
C= A Bs*). 
ee s+ 5 Bs") 


Substituting for C and D in the expression for the 
potential, 


l Aa te = 


4 6 
+E (past $Bs" 


ue 1 
— V;, 14 2 
- = E (a — sx) 


+ £B(2? — 8? IF 


The electric field is given by 


Popa lo. 
dx 
Vs 1 |, 
A(2x 
( EE0 E eas) 


+18 (30? — )}) eL, 


Solution 4.3. (a) The electrostatic potential V in 
the dielectric material satisfies Poisson’s equation 


V7V =—— pp. 
EEQ) 


The system has spherical symmetry, and pp = po, So 
this reduces to 


Ld 9 dV 1 
Tr — 
r2 dr dr E£0 ea 
(b) Integrating, we obtain 
2dV spo 
ti ——s 


3 
= —- A 
dr 3EE0 eee 


(a<r< 2a). 


po 2 A 
V(r) = -—~—— r° -—-—-4+B, 
(7) 6E€o v = 


where A and B are constants of integration. 


(c) In order to find the constants A and B we use the 
boundary conditions V(a) = 0 and V(2a) = 0. So 
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and 


2p0 2 A 
V(2a) =0=— B. 
a 3EEQ “ 2a bi 
3 7 
Hence A= - and B= _ 
EEQ 6E€o 
Therefore 
3 2 
Po 9  poa 1 Tpoa 
V(r) = 
(r) 6E€o 7 669 6E€g 
3 2 
-- 22 f+o(4)'-7(2)] 
6EEo r r 
(a<r < 2a): 


Solution 4.4 The conducting plate can be ‘replaced’ 
by two image charges, each —gq, located at points 
(—d/2, d/2, 0) and (—d/2, —d/2, 0) (see 

Figure 12). 


—df2 d/2 © 
eee Ot¢@ 
O charge conducting 
O image charge plane « =0 


Figure 12 For Solution 4.4. 


By symmetry, the real charges and the image charges 
produce zero electrostatic potential in the plane x = 0 
midway between them, and the potential very close to 
each of the real charges will be unaffected by the 
image charges. So by the uniqueness theorem, since 
the boundary conditions are the same, the solution for 
the electrostatic fields and potentials in the region 

x > 0 for the problem with real and image charges 
must be the same as the solution for the problem with 
the real charges and the conducting plate. And since 
the force on charge gq at (d/2, d/2, 0) is just the 
charge times the electric field at that point, we can 
determine the force by adding the Coulomb forces due 
to the other real charge and due to the two image 
charges. 


The electrostatic force on g at (d/2, d/2, 0) is 
therefore (see Figure 13) 


F = F(due to charge + ¢ at (d/2, —d/2, 0)) 
+ F(due to charge — g at (—d/2, +d/2, 0)) 
+ F(due to charge — g at (—d/2, 
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d/2, 0)). 


YA 
pele 


| > 
aj2 * 
Figure 13 For Solution 4.4. 
These three contributions to the force are in 
the directions of the unit vectors e,, —e, and 
—(ez + ey) / V2, respectively, and the separations of 


the charges are d, d and \/2d, respectively. So the 
force is 


_¢ iL. I 1 (e; + ey) 
- Ameg |d2 4% dd? * Qd2 4/9 


at (+a) e+ O-aa) >] 


The magnitude of this force is 


(+3) *(-2) 
“fall ded) (se) 


_ 37 
a 87eEd2 ; 


g 


~ Ate d? 


Solution 4.5 In cylindrical coordinates, 


2,, 10 { OV Levy ov 
eae "Or geo on 
For V(r, ¢, z) = Arcos ¢+ (Bcos ¢)/r, 
OV —— 
Bp = Acoso - 
3) OV O Boos 
aG ~)= Or (Arcos = r aoe) 
= Acos 6 + =SS°, 
OV . Bsing 
a as eae a 
OV Boos 
CT a ao. 
av 
O22 
So 
VV =- = (Aces o4 Bese) 


il 
sae (— Arcos - sone) +0=0, 
r 


as required. 


(b) We use cylindrical coordinates with the z-axis 
along the axis of the cylinder and the line ¢ = 0 in 
the direction of the uniform electric field Ep. The 
potential corresponding to this uniform electric field is 
then V = —Eor cos ¢. 


We shall denote the potential outside the cylinder 
by V(r) and inside the cylinder by Vo(r). These 
potentials must satisfy Laplace’s equation with the 
following boundary conditions: 


G) V, ~—Eorcos@ for large values of r; 
Gi) Vi=VW forr=R; 

(iti) OV, /Or =eOV2/Or forr = R; 

(iv) V2 remains finite for r = 0. 


Part (a) and boundary conditions (i) and (iv) suggest 
that we look for solutions of the form 
B 
VY, = —Eprcosé+ ee 
r 


V2 = Arcos@ forr < R. 


The coefficient of r cos ¢ in V; ensures that boundary 
condition (i) is satisfied, whereas the absence of a 
term proportional to (1/7) cos ¢@ in V2 ensures that 
boundary condition (iv) is satisfied. It remains to 
choose values for the constants A and B so that 
boundary conditions (ii) and (iii) are satisfied. Using 
boundary condition (ii), we obtain 


forr > R, 


B 
—EogR+ = = AR, 
oft + R 
whereas boundary condition (iii) leads to 
B 
—Eo — R? = cA. 
Solving these two equations, 
2 e—l 
A=——"_ 8 d B= EqR? 
eer a 
So the electrostatic potential is 
—1 
VY, = —Eorcos@ + a EoR? ey forr > R, 
e+1 
2 
Le ter forr < R. 


Note that the electric field inside the cylinder is 
uniform, but reduced by a factor 2/(€ + 1) compared 
to the applied electric field. 


Solution 4.6 Amending the argument used to derive 
Equation 4.26 of Book 2, the general finite difference 
form of Laplace’s equation in three dimensions is 


v= 
V (i415 9p, 26) — 2V (as Uy 2a) HV (i215 Wy Se) 


(Az)? 
n V ig Up, 2g) — 2Y (a te) EV a, is Se) 
(Ay)? 
- V (i, Yj, 2e41) — 2V(@i, Yj, 2%) + V (Gi, Yj, Ze-1) 
(Az)? 


= 0. 


If we take the step lengths Az, Ay and Az to be 
equal, this reduces to 


1 
V (4,942) = 6 Vein wy 2) 
+ V (a5, Ups Be) PV (eis. e448) 
+ Vig Yj Zn) =i eee Yj-ls Zn) 
a V (cay te-2)] 


This equation says that the value of V at any grid 
point is the average of the values of V at the six 
closest neighbouring grid points. 


Solution 4.7 We use Equation 4.27 of Book 2 
iteratively to improve on the initial guesses for the 
potentials at the nine internal grid points. For the four 
corner points, 


V =7(0+140.5 40.5) = 0.5. 
For the midpoints at the top and bottom, 

V = 4(0+0.5 + 0.5 + 0.5) = 0.375. 
For the midpoints on left and right, 

V =4(1+0.54+ 0.5 + 0.5) = 0.625. 
For the central point, 

V = 7(0.5 + 0.5 + 0.5 + 0.5) = 0.5. 


These results are shown below. 


0 0 0 
1 0.5 0.375 0.5 1 
1 0.625 0.5 0.625 1 
1 0.5 0.375 0.5 1 


0 0 0 


A second iteration leads to the following values for V. 
For the four corner points, 


V = 7(0 + 1 + 0.625 + 0.375) = 0.5. 
For the midpoints at the top and bottom, 
V = 7(0+0.5 + 0.5 + 0.5) = 0.375. 
For the midpoints on left and right, 
V =7(1 +054 0.5 + 0.5) = 0.625. 
For the central point, 
V = 7(0.375 + 0.625 + 0.375 + 0.625) = 0.5. 


The values after the second iteration are identical to 
those after a single iteration. So the Jacobi relaxation 
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method has converged and no further iterations are 
necessary. 


Book 2 Chapter 5 
Solution 5.1 The field due to a current element is 
given by the Biot—Savart law: 
poo Idi’ X (r-1r’) 
An 
In this case, 
Tol =50% 10 * 6, Am, 
r = (0.08 e, + 0.08 e,,) m, 
r’ = 0.02e, m, 
r—r’ = (0.06 e, + 0.08 e,) m, 
lr — r’| = 0.062 + 0.082 m = 0.1m. 
Also, 
I6Y X (r-r’) 
= (5.0 x 107° ey X (0.06e, + 0.08 e,)) Am? 
= 0.3 x 10-8 (e, x e,) Am? 
= —0.3 x 10-°e, Am?. 


dB(r) 


jr —r'|3 


Substituting these results into the Biot—Savart law, 
po Tol’ X (r—1’) 
Ar 


dB(r) 


jr —r'|3 

(—0.3 x 10-®e, Am?) 
(0.1 m)3 

=—3x 10 "e«, NA ‘m’ 

=-3x107"e, T. 


=10" NA 


Solution 5.2 By symmetry, each side of the hexagon 
will produce a field of equal magnitude at the centre, 
and all of the contributions will be in the same 
direction, perpendicular to the plane of the loop. The 
magnitude of the field produced by one side is given 
by Equation 5.8 in Book 2: 


_ Hol 
~ And 
where the distance d between field point and wire in 
this case is Lsin 7/3 = Lsin 60° (the height of 
an equilateral triangle with side length L), and 


[sin ap — sin aa], 


ap = —aa = 77/6 = 30°. So 
pol : ° 
on a rp, 
fia 
I 
An Ly/3/2 
Hol 


~ 2/3nL 
So the total field at the centre of the hexagon has a 
magnitude 


Hol V3 p01 


Bhex = 6 
" 2/3rL TL 
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Solution 5.3. (a) The magnitude of the axial field 
in a solenoid of finite length, with number of turns per 
unit length n, is given by Equation 5.13 in Book 2: 


I 
Fs ees [sin ap — sin aa]. 


In the middle of a long solenoid, for which we assume 
that the diameter is very much less than the length, 
aa = —7/2 and ag = 77/2, so 

B= pon, 
which is a standard result. At one end of a long 
solenoid, a, = 0 and ag = 7/2, so 

B= pon! /2, 
which is half of the value obtained for the field in the 
middle. 


(b) In the middle of a solenoid whose length is 10 
times its diameter, 
10 10 


Vi+i02 viol 


sinap = —sina,a = 


So Equation 5.13 gives 
_ pont 20 


Be ~~ 0.995 pion, 
2/101 rs 


which is only 0.5% less than the result for a solenoid 
of infinite length. 


Solution 5.4 (a) 


B=culA=|2 2 2 


= a(zey, + rey + yez). 
(b) Using the differential version of Ampére’s law, 


1 
J = —curlB 
Ho 


Gz Cy ee, 


_ i () ie) (a) 
Ox Oy Oz 


Ho 
az ax ay 
1 
= —(a+a+a) 
HO 
_ 3a 
Lo 


Solution 5.5 (a) As discussed in Subsection 5.3.2, 
for the spherical-head approximation, radial dipoles 
and their return currents give rise to fields outside the 
scalp that are tangential to the surface of the scalp. 


A dipole that is perpendicular to the surface is a 
radial dipole, so in this case there will be no radial 
component to the magnetic field just outside the scalp. 
Since the magnetometer only measures the radial 
component, the measured field will be zero. 


(b) A dipole parallel to the surface is a tangential 
dipole, and the radial component of the magnetic field 
outside the scalp is identical to the radial field due to 
the dipole alone; there is no contribution to the radial 
field from the return currents (see discussion in 
Subsection 5.3.2). The magnetic field produced by the 
dipole is in the azimuthal direction about the dipole 
axis, that is, the field lines are circular. Directly above 
the dipole the field is parallel to the surface of the 
scalp, and so there is no radial component of the field, 
so again the magnetometer will not register a field. 


(c) A dipole J 6] at 45° to the surface can be 
resolved into component dipoles of equal magnitude, 
Tél sin 45° = I6l cos 45°, in the radial and 
tangential directions. Neither of these component 
dipoles produces a radial component of magnetic field 
at the point directly above it, so the magnetometer 
does not register a field from the inclined dipole. 


Solution 5.6 The measured radial field is that due to 
the dipole alone, with no contribution from the return 
current. As discussed in the solution to Exercise 5.5, 
the field lines are circles, centred on the axis of the 
dipole. So on one side of the highest point of the head 
the radial field will have an outward component and 
on the other side it will have an inward component. A 
line joining the position of the peak in the radially 
outward field to the position of the peak in the 
inwardly pointing field will be perpendicular to the 
direction of the current dipole, and the separation of 
the peaks will increase as the depth of the dipole 
increases. The magnitudes of the peak fields is a 
measure of the strength of the current dipole. (See 
Figure 5.25 and accompanying discussion.) 


Book 2 Chapter 6 


Solution 6.1 Problems like this require a systematic 
approach to keep track of the different components of 
the forces. The resultant force is 


F = —e(E+v X B). 


Now ex Cy &z 
vxXB=2.0x10°ms!x2.5T/0 1 1 
1 0 1 


= 5.0 x 10°Tms ‘(e, + ey — ez) 

= 5.0 x 10°Vm ‘(ez + ey — ez), 
since the unit of v X B must be the same as the unit 
of E. So 
F =—1.60 x 10-9 (2.5 x 10°V m-! (e, + 2e,) 


+5.0 x 10°V m7!(e, + ey — e2)) 
= —4.0 x 10713N (3e, + 4e, — 2e,). 


The magnitude of the vector (3e,; + 4e, — 2e,) is 


V/32 + 42 + 22 = \/29, so the magnitude of the force 
iS 


F=40x107N x 29 = 2.2 x 107!2N. 


Solution 6.2 Att = 0 the particle is at rest so it 
cannot experience a magnetic force. The initial 
acceleration is therefore due to an electric force, Fetec. 
The electric field must be in the same direction as the 
acceleration of the positively-charged proton, that is, 
in the +z-direction, and its magnitude is 


F. elec 


€ 
mag 


f= 


e€ 
1.67 x 10-2" kg x 9.6 x 10'4ms~? 
1.60 x 10-19C 
=10*10'Vm"". 


At time 1, the electric force is unchanged. 

The acceleration of the moving proton has the 
same magnitude as it had initially but the 
direction has reversed. The magnetic force must 
therefore be in the opposite direction to the 
electric force and have twice the magnitude of the 
electric force, that is, Fmag = —2maoez. Since 
Frag = qv X B = e(—v e,) X B, we can write 


—2mapez = e(—vie,) X B = —ev;(—B,e,+ B,ez). 
Equating x-components, we deduce that 


2mag 


B,= 
eV] 


_ 2% L67 x 10°?" ke x 9.6% 10" ms 
1,60 % 10-"'C % 2.0-~ 10° ims=* 
=1,07, 


and comparing z-components of the force we deduce 
that B, = 0. 


The data in (i) and (ii) do not tell us anything about 
B,,. However, if there is no motion in the z-direction, 
then there is no force in this direction, which means 
that the magnetic field must be in the z-direction, 

and so B, = 0. If B, were not equal to zero, then 
immediately after release, when the particle is moving 
in the x-direction, there would be a component of 
magnetic force in the z-direction, which conflicts with 
the observation that the proton does not move in the 
z-direction. 


Solution 6.3 Ifa particle, which has charge q and is 
travelling at velocity v, is not deflected, then the net 
force acting on it is zero, so the electric and magnetic 
forces must cancel. Thus 


qE = —qv XB, 


33 


E=-—v X B=~— ve, X Be, = —uBe, 
= —(3.0 x 10’ms7! x 0.50T) ez 
= —-15x10'Vm te. 


Note that the expression for E does not depend on the 
mass or charge of the particle, so particles of any 
mass and any charge that are travelling with velocity 
3.0% 10'ms e, will be selected. 


Solution 6.4 The magnetic field does not change the 
proton’s speed, which is 3 x 10’ ms~!. The period of 
its motion is 

27 2am 


We eB 


_ Wn x 1.67 x 107?" kg 
~ 1.60 x 10-19C x 2 x 10-19T 


~ 300s. 


The diameter of the orbit is 
2mv 


— |g/B 


_ 2X 167 «10 ke x 9% 10" mis * 
160% 10-C X 2 x 10-T 


~3x 109m. 


d= 2r. 


This orbit is about 10 times smaller than the diameter 
of the Earth’s orbit around the Sun, but it is traversed 
in about 5 minutes. Note that even a weak field can 
have a substantial effect on a high energy particle. In 
practice, the cosmic ray protons would follow a 
helical path about the magnetic field direction, rather 
than a circular path. 


Solution 6.5 (a) The force on an ion is given by 
F = Q(E+ v X B). Since 


vxB=|2 & az) — Be, + Boe.. 
0 B O 
Thus F= Q(-BSe, + Bey + Be.). 
Applying Newton’s second law to the ion, we have 
dr _ 
me) 


which in component form gives 


(G d?y d?z ) 
m + + 


dq ae ae 
dz da 
=() (- Ge. + Bey + Be.) ‘ 
Equating components on either side, 
d?x _ QBdz dy _ QE d?z — QBdx 
dt? =m dt’ dt? =m’ dt? — m dt’ 
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as required. 


(b) We solve these equations of motion using the 
initial conditions r = 0 and dr/dt = ve, when 
t =0. Solving the equation for the y-component by 
integrating twice gives 
y= QE 12 F 
2m 


Integrating the equation for the x-component of the 
acceleration, and noting that daz /dt = 0 when z = 0, 
we obtain 


Substituting this expression for da /dt into the 
equation for the z-component of the acceleration gives 


d?z QB 

— =-—(— } z. 

dt? m 
This is the equation of simple harmonic motion, with 
the solution 


—_ mv. (QB 
zZ= OB sin (Sr), 


which satisfies the boundary conditions z = 0 and 
dz/dt = v at t = 0. Solving for x, 


dx QB . (2 ) 
= z=-—v sin | —tT], 
dt m m 
and so 
~ QB m QB 


mv QB , 1 
= — |cos | —t]—-l], 
QB m 
where we have used the initial condition that x = 0 
when t = 0. 


(c) An ion hits the grid when z = 0. From 

the solution for z we see that z = 0 when 

sin QBt/m = 0. Ignoring the t = 0 solution, the next 
time that the ion passes through the z = 0 plane is 
when QBt/m = 7, or t = 1m/QB. At this instant, 


mv 2mv 
= apr l= — OB’ 
and 
QE (mm 2 mE 
Y= Om (5) ~ 9QB?’ 
as required. 


Note that the y-displacement depends on the 
charge-to-mass ratio of the ion, whereas the 
x-displacement also depends on the initial speed. 


Solution 6.6 For fields with varying direction, the 
drift velocity is given by 


where Rp is a vector drawn to a point on the field line 
from the centre of curvature of the line at that point. 
The field lines for a long straight current-carrying wire 
are circles centred on the wire. So the electron will 
spiral round a circular field line centred on the wire, 
which means that Rg will remain constant, and 
Rep = e,. The field produced by the wire is given by 
the standard expression (see Figure 5.2a of Book 2) 
Hol 


=e 
onr ? 


soRz xX B =e, X eg = e;. So 


mvj _ moj 27Rp 
qBRp 


2nmvy 


a e 
“s —eRp pol ~ puoel 
For the drift due to the varying strength of the field we 
use Equation 6.20: 
2 
muy 
Vd = gps x grad B. 

Since B = pol /2zrr, using the expression for grad in 
cylindrical coordinates, 


dB, Mol 


grad B = ap e, = 


Hence 
B 
B X grad B = Beg X | —— e; | = — ey: 
r r 
So for an electron spiralling around a field line with 
radius Rp, the drift velocity is 


_ mut Be —_ mu a 
—2eB3 Rp ~ = 2eBRp ~’ 
and substituting B = pio /27Rzg we obtain 


Vd 


mmo 


e 


Vqd=—- : 
: pioel 


Note that since vy) = UL in this case, the drift due to 
eld curvature is twice that due to the gradient o 
8 
eld strength, but both contributions are in the 
8 
—e,-direction. 


The magnitude of the drift velocity is given by 


2rmuy 


= aX 
“ Hoel 


_ ood x10 ke & 0 x 10? ms? )* 
~ 2x 10-7NA~? x 1.60 x 10-19C x 100A 
—43x10?ms—!. 


Note that the drift velocity is independent of Rp, the 
radius of the orbit about the current. 


Solution 6.7. (a) The electron’s speed 
perpendicular to the magnetic field initially is 
VL = Ug = 3.0 x 10° ms~. Its cyclotron radius is 


therefore 
MUL 


lq\B 
_ 9.11% 10-9 ke x 3.0.x 10°ms* 
~ 1.60 x 10-19C x 2.0 x 10-3 T 
= 8.5 x 1073 m. 


Te 


(b) The quantity (sin? w)/B remains constant as the 
particle moves through the bottle. The initial value ag 
of the pitch angle is 


ag = tan~'(vg/vz) = tan~1(3/2) = 56.3°, 

and the initial field strength is 2.0 x 10~° T, so 
sin?a sin? 56.3° 

B 2.0% 10-?T 

At the mirror point, @ = 90° and sin? a = 1, so 
B=1/(3.46 x 10°T™) =2.9 x 10-*T. 


— 3.46 x 107T-!. 


(c) At the mirror point, v, = 0. Since the electron’s 
energy is conserved, v, must now be equal to the 
initial speed of the electron, that is, 


v, = V3.02 + 2.02x10°ms—! = 3.6x10°ms!. 


So the radius of the electron’s orbit at the mirror point 
is 


— MvUi 
re q/B 

_ 911x107" ke x 3.6 x 10°ms* 
~ 1.60 x 10-19C x 2.9 x 10-3 T 
= 7.1x 1073 m. 


(d) The value of sin? a must be the equal to its 
initial value when the electron returns to the z = 0 
plane, that is sin? a = sin? 56.3°. This has solutions 
a = 56.3°, (180° + 56.3°), (360° — 56.3°). The 
appropriate solution, corresponding to the electron 
continuing to orbit in the same sense about the z-axis, 
but with the z-component of velocity reversed, is 

a = 123.7°, and the velocity is 


v = (3.0eg — 2.0e,) x 10°ms7?. 


Book 2 Chapter 7 


Solution 7.1 The planar symmetry indicates that the 
electric field and the current density are everywhere 
perpendicular to the plane of the slabs, and we take 
this to be the z-direction. In a steady state, the same 
current flows across any plane with constant z within 
the slabs, and therefore the current density, J,, say, is 
the same at all points in the slabs. The electric fields 
within the two slabs, F’.; e, and FH, e., are then 
given by 


J, = 0, Ey = 02 E22. 


The potential difference across the two slabs is (see 
Figure 14) 


AV = E.,d, + Exod. 
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From the first equation, F,2 = F,01/o2 and 
substituting this expression into the second equation 
leads to 


AV 
dy + (01 /02)d” 
Now the potential difference across slab | is F,,d1, 


and since the outer boundary of this slab is at potential 
V = 0, the potential at the interface is 


Vinter — Endy 
AV 
dy + (01 /02)d2 
_ AV 
1 + (o1d2/02d1) , 
If d; — 0 or dg — ~, then Vintey — 0. 
If dj — o0 or dy — O, then Vinter — AV. 


If 0, — O or ay — ~, then Vinter — AV. 
If oj — wor og — 0, then Vintey — 0. 


Ea = 


All of these limiting values are in agreement with 
what would be expected. 


aA A 


ds slab 2 AV 


dy slab 1 


Figure 14 For Solution 7.1. 


Solution 7.2 We denote the length, radius and wall 
thickness of the tube by /, r and t, respectively. For 
each pattern of current flow, resistance = ‘length’ / 
(‘area’ x conductivity), where ‘length’ is measured in 
the direction of current flow and ‘area’ is measured 
perpendicular to current flow. See Figure 15. 


(a) For current flow along the length of the tube, 
‘length’ = 1, ‘area’ = 2mrt, 


sO 
l 


2rrto 


0.1m 


Qn x 5.0 x 10-3 m x 5.0 x 10-4m x 1.0 x 1070-1 m- 


=64x1074* 9. 
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(c) 
Figure 15 For Solution 7.2. 


(b) For current flow between the inner and outer 
surfaces of the wall of the tube, ‘length’ = f, ‘area 
= 2rrl, 


’ 


5.0 x 10-47 m 
Qn x 5.0 x 107-32 m x 0.1m x 1.0 x 107Q-! m7! 
—16x 107° Q. 


(c) For current flow around the wall of the tube, 
‘length’ = 27r, ‘area’ = It, so 


27r 


lto 
Qn x 5.0 x 107-2 m 


0.1m x 5.0 x 10-4m x 1.0 x 1070Q-! m7! 
=6.3x 107° 0. 


Re = 


Solution 7.3. The equation above Figure 7.10 of 
Book 2 gives the potential along the perpendicular 
from the small sphere to the metal surface. The 
potential difference between the sphere and the 
surface is obtained by integrating this potential along 
the perpendicular line: 


I a 1 a 7 
4no Jn \e 2d-z : 


= ; [In z + In(2d — 2° 


TO 


ok int 41 d 
Ano a a 


AV = 


The resistance is therefore 
AV if d? 
R = = | . 
I dno (sais) 
For d/a = 10ando = 5.007'm"}, 
1 100 


R= 1 = 26x 10779. 
dn x5.00-!m— 19 ‘ 


Solution 7.4 As discussed in the solution to 
Exercise 7.6 of Book 2, the increasing field of the 
solenoid induces a current in the ring that creates a 
magnetic dipole in the opposite direction to the field 
of the solenoid. In the uniform field at the centre of 
the solenoid there is no net force on the dipole — just 
as there is no translational force on a bar magnet that 
is aligned with a uniform field. When the ring is at the 
end of the solenoid, the field around it is not uniform, 
and there is a net repulsive force, as described in the 
solution to the exercise, but the force is much smaller 
in the absence of the iron core. 


The effect on the ring can be understood by 
considering the magnetic forces that act on the 
currents induced in the ring. The magnetic force 
acting on a current element J 61 in a field B is 

F = 161 X B. When the ring is in the uniform axial 
field near the centre of the solenoid, the vector product 
of J 61 (in a clockwise azimuthal direction viewed 
from above) and B (vertically upwards) is radially 
inwards at all points on the ring. So there is no net 
force on the ring when it is near the centre of the 
solenoid. At the end of the solenoid the field is not 
uniform; it decreases with distance and the field lines 
spread out — they have a component in the radial 
direction, away from the axis of the solenoid, in 
addition to the axial component. The vector product 
of I 61 with this outward radial component gives an 
upwards force at each point on the ring. However, 
without the magnetic core, the force is likely to be too 
weak to propel the ring upwards. 


Solution 7.5 The flux through the coil when the 
normal to the coil is at an angle @ to the direction 


of the horizontal part of the magnetic field is 

® = NABhoriz cos @. Note that the vertical 
component of the field doesn’t contribute anything to 
the flux through the coil, since it is always parallel to 
the rotating coil. The angle @ varies as 0 = 27 ft. So 


Vent = —d®/dt = 2nfNABporiz sin 27 ft, 
which has an amplitude 

Vo = 2a fN ABhoriz, 
and hence 

Bhoria = Vo/20fNA. 
For the data quoted, 

Bhoriz = 25 x 1073 V/(2m x 80Hz 


x 50 x 0.50 x 1074 m?) 
=0.020T. 


Solution 7.6 (a) The emf, Vp, induced in coil B, is 
given by 
dl 

Va = —-M es = —wM Tyo coswt. 
Since Rp >> wp, the current in coil B does not 
depend on the self-inductance, so 
Vp _ wM Ino 
Rp Rp 
Comparing this with the expression given for the 
current in coil B, we see that 

wMIag _ 

Ra 


cos wt. 


n= 


IpoRp 
wag — 


Igo, so M= 


(b) The value of the mutual inductance between the 
two circuits is the same when a current in B induces a 
current in A as when a current in A induces a current 
in B. So in this case, 


dJp — IpoRp 
Va=—-M oi algo exp [—at]. 
Again, I, = Va/Ra, so 
alp,lpoRp 
== i]. 
= wlragoRa p= et 


Solution 7.7 To determine L we need to find the 
flux through the coil for a given current, and to 
determine the flux we need to know the magnetic field 
in the core. We assume that the magnetic flux is 
confined to the core. Symmetry indicates that the 
magnetic intensity is in the azimuthal direction for 
cylindrical coordinates with their z-axis perpendicular 
to the plane of the disc and through its centre, and that 
its magnitude is constant around a circular path with 
radius r within the core. So applying Ampére’s law to 
such a path, 


H(r)2ar = NI, and 
B(r) = peo (r) = ppoNI/2rr. 


37 


The flux through each turn is the same, so the total 
flux through the coil is N times the flux through a 
single turn. Within a single turn, the field is constant 
within a thin strip, of width dr that is parallel to the 
z-axis and at a distance r from it, and the field is 
perpendicular to the plane of the strip, so the flux 
through the strip is 


6® = Bir) adr, 


and the total flux through NV turns is 


NG rd N?I 
ee aa af dp o n2. 
20 a T 20 
d® — pipig Na 
Th L= — = —— ]n?2. 
re dI Qa . 
Substituting the data given in the question, 
L= 
500 x 4a x 10-7 NA x 10° x 4.0 x 10-7 m ao 
n 
2a 
= 2.8 x 10-°H. 
Book 2 Chapter 8 


Solution 8.1 From Equation 8.4 of Book 2, 


4 4 
_1 _ q 
U= ad ai, where Vj = | 3 wee 
i=l J=1,j#% m 
The symmetry of the problem indicates that the 
potentials V; at the location of each of the charges are 


identical. For a particular charge 7, the distances to the 
other three charges are /2a, V2a and 2a, so 


(d) Resonant frequency is 


= 
i= 5. 


1 


2rV/ LC 
1 


~~ On/DBHX LOE 


= — = 0.064 Hz. 
27 X 2.58 


Solution 8.3 The required ratio is the ratio of the 
electric and magnetic energy densities. For air, ¢ = 1 
and yw = 1, so 
Uelec = sE0E” 
= 0.5 x 8.85 x 10712 C? N71 m-?(100 Vm‘)? 
=A42 % 10°? 4, 
BP _ 60 310-"T)* 


= 9.95x10~4J. 


U = = 
ms 20) 8 X 10-7NA~? 
So the ratio is 
Uelec _ 4.43 x 107° J _ 44x 1075, 


itnse BOS 10-*] 
or 44x 107° :1. 


The magnetic field energy is clearly very much larger 
than the electric field energy. 


Solution 8.4 (a) During the charging phase, 


Vo = Ve c — exp (-z5)| 


Att = 5s, Vo/V, = 0.95, so 


and 


q q q 1 Thus 
— “Geyia ea  Anega (v2+ 3) z t ee 
Thus ~  ©lm0.05° 1.0 x 10-3 F x In 0.05 
2 = 1.7k2. 
U = 4x4xqx— (v2+3)=—*_(2v2+1), 
4rega 4rega 


where the factor 4 takes account of the fact that there 
are 4 identical charges. 


Solution 8.2. (a) U = 4CV?, so 
WU 2x5.0 x 10°F 


C = = 2.5F. 
V2 (2.0 x 10 V)2 
(b) U=4LI",so 
2 2x 5.0 x 108 
= U _ x 5.0 x OS oan 


P (2.0 x 103A)? 


(c) For the RC circuit, time constant 
= RC=0.20 x 2.5F=0.5s. 

For the RL circuit, time constant 
=[L/R=25H/020=— 1258: 
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(b) Average power P in discharge time T is given by 


p 09% 3CV? 
T 
_ 0.9 x 0.5 x 1.0 x 10-3 F x (300 V)? 
7 2x 10-35 
= 20kW. 


Solution 8.5 The charge will be distributed 
uniformly over the surface of the conducting sphere 
(since the charge density is zero inside a conductor). 
The symmetry then indicates that the electric field is 
spherically symmetric, with E = E,(r) e,. 


Using the integral version of Gauss’s law with a 
spherical Gaussian surface, radius r (> R), 
E, x 4nr? =Q/eo, so Ep = Q/Aneor?. 


Alternatively, you could justify this result by saying 
that outside any spherically symmetric charge 


distribution the field is the same as if the charge were 
all located at the centre of the distribution. 

The total energy U stored in the field is obtained by 
integrating the energy density, seo E”, over the region 
outside the sphere. Since the field strength is uniform 
within a thin spherical shell, radius r, thickness dr, 
volume 47r? dr, the energy dU in such a shell is 


Q 2 
=a 2 
6U = 5£0 (25) Arr or, 
and the total energy is 
2 oO 7 2 
U = —— | =) dr = J. 
87ée0 JR 87egoR 

as given in the question. 


In Exercise 8.2 of Book 2 we showed that the potential 
energy of a conducting sphere with surface charge 
density o is 


E0 
or, since o = Q/4rR?, 
oO 
= 8regR’ 


which is the same as derived for the energy stored in 
the electric field. 


This again indicates that there are two ways of 
thinking of the electric energy: we can think in terms 
of the energy of charges held at various potentials or 
we can think in terms of the energy stored in the 
electric fields in the regions around the charges. Both 
viewpoints lead to the same result for the energy 
stored in the system. 


Solution 8.6 (a) Since the total charge @ in the 
sphere is related to the (uniform) charge density p by 
Q= on R? p, the expression for the potential energy 
derived in Worked Example 8.1 of Book 2 can be 
written as 


7 4rR® ( 3Q.\? — 3Q? 

~ 15€9 (5) — 2repR’ 

(b) The electric field outside the sphere is the same 
as that outside the conducting sphere in the previous 
question — and the same as in the region r > R 
around a point charge @. So the energy stored in the 
electric field in this region is the same in all three 


cases, and given by the expression verified in Exercise 
8.5. 


2 


8regR 


UrsR =— 


(c) The difference between the expressions in (a) 
and (b) is 


2 2 
AU = g a g : 
TeEgR 20 8 40T7eE9R 


This must be the energy stored in the electric field 
inside the sphere, that is, in the region r < FR. In 
Worked Example 8.1 an expression was derived for 
the field inside the sphere: E,. = pr/3eo. Thus, using 
a similar integration method to the previous question, 
the energy stored in the field inside the sphere is 


R pr 2 
U,.2n = seo [ (=) Anr? dr 
) E0 


R 
= ae f r! dr 
0 


9Eo 
__ &np? R? 
ABE | 


Since p = 3Q/47R?, 

2k? ( 3Q \?_— Q? 

A5eg \47R3] = A0me9R’ 
confirming the suggestion that the difference between 


the energies in parts (a) and (b) is the energy stored in 
the electric field within the sphere. 


Ur<R = 


Solution 8.7 The self-inductance of a solenoid with 
length /, area A and number of turns per unit length n 
is given by Equation 7.26 of Book 2: 
L=on7lA 
=4n x 107-7’NA~? x (200m7!)? 
x 1.0m x 7(0.10m)? = 1.6 mH. 
Thus the energy stored is 


US iLP ]05%16410 “Ax (100A)? = 7.9). 


The power dissipated by Joule heating is given by 
P = I*R, where the resistance R is 


wire length 


o X wire area 
200 x 27 x 0.10 m 


~ 6.5 x 107Q-1m-! x 4 (2.5 x 10-3 m)2 
= 9.8 x 107-72. 


So the power dissipated is 
P= I’R= (100A)? x 9.8 x 107-7 = 980 W. 


This is clearly not a practical proposition. The energy 
dissipated per second, 980 J, is over 100 times greater 
than the energy stored! Also, the energy stored is tiny 
— the amount required to lift a 1 kg bag of sugar 
through 0.8 m. 


Book 2 Chapter 9 


Solution 9.1. (a) The coil forms a continuous 
superconducting path, so the flux through it must 
remain constant when the magnetic field is reduced to 
zero, and the value of the flux must equal the flux 
through the coil due to the uniform applied field. This 
is simply 


® = NABo = N (ra”) Bo = tNa? Bo. 
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(b) The self-inductance of the coil is defined by 
L=<d@/dI,so L = ®/I, and 

I 6 z7Na?Bo 

Ei * 

(c) The field strength at the centre of a circular coil, 
with N turns of radius R, carrying current J, is 
Beentre = poNI /2a. 
(This can be deduced from the Biot—Savart 
law: a small element of length 61 contributes 
OBeentre = Lol 61/4ra?, and since the fields due to 
all elements of the coil are in the same axial direction, 
the total field is just N2na/6l times greater than this.) 
Thus 


B WoNI poNaNa?Bo  ponN?aBo 

ca, a i |e’) aa 
(d) The symmetry indicates that the field in the 
plane of the coil must be perpendicular to this plane. 
Because a point at r = 0.95R is so close to the wire, 
the field there will be larger than at the centre. Note 
that though the field originally applied to the coil was 
uniform, it is the magnetic flux through the coil that 
must remain constant, not the magnetic field. The 
magnitude of the field near the centre will be smaller 
than that of the applied field Bp and the field close to 
the turns of the coil will have a larger magnitude. 


Solution 9.2. (a) From Table 9.1 of Book 2, 

T. = 3.4K for indium, and B,(0) = 0.028T. The 
temperature T’ that corresponds to a critical field 
strength of B.(T) = 0.014 T is given by Equation 9.1: 


BAT . 
( )_0.0MT 9. int © 
B.(0)  0.028T 


So 
(T/T. =1-05=0.5, and 
T= ¥V0.5T, = V0.5 x 34K =24K. 
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(b) The critical current produces a magnetic field at 
the surface of the wire that is equal to the critical 
magnetic field. At a temperature of 1.2 K, the critical 
field is given by 


B,(1.2K) = 0.028T [1 — (1.2/3.4)?] = 0.0245T. 


We can neglect the Earth’s magnetic field, since it is 
very much smaller than the critical magnetic field. 
The magnetic field strength B at the surface of 

a wire of radius a carrying current I is given by 

B = ol /2z7a, so the critical current corresponding to 
a critical field of 0.0245 T is 


_ 2naBe _ 20 x 1.0 x 1073 m x 0.0245T 
Lo An x 10-7NA~? 
= 120A. 


I, 


Solution 9.3 Since €, > Aq, material A is a type-I 
superconductor, whereas for material B, £8 < Ag so it 
is a type-II material. In both cases, for very small 
applied magnetic fields, the field will deviate around 
the disc, as shown in Figure 9.21b in Book 2, and a 
superconducting screening current will flow around 
the perimeter of the disc. As the field is increased it 
eventually penetrates into the discs, destroying the 
superconductivity in the regions that it penetrates. For 
disc A, the magnetic field penetrates as ‘swirling 
bands’ (Figure 9.22), typically about a millimetre 
wide, with the proportion of the material that is in 

the normal state increasing as the magnetic field 
increases. For material B, the field penetrates as thin 
cylindrical cores, within which the material is normal. 
As the applied field increases, the number of normal 
cores increases and they become more closely packed 
together. Eventually a field is reached at which the 
cores merge and all of the material is normal. 


Solution 9.4 Table contrasting type-I and type-II superconductors. 


property Type-I Type-II 
type of material elements alloys, ceramics, etc 
critical temperature low, generally < 7K alloys up to 30 K, 
ceramics (HTS) up to 140K 

refrigeration liquid helium liquid nitrogen for HTS 
critical field single critical field, lower and upper critical fields, 

<0O.1T upper B, up to 20T 
applications few many 
E/X > 1 <1 
surface energy of NS boundary _ positive, negative, 


NS boundary area minimized 


field penetration intermediate state, 


coarse layers of N and S, 


parallel to field, 


scale typically ~ mm 


NS boundary area maximized 
mixed state, 

narrow cylindrical N cores, 
parallel to field, 

scale typically < 0.1 nm 


Solution 9.5 (a) Assuming that the penetration 
depth \ is much smaller that the radius R of the 
cylinder, the field will fall exponentially with distance 
below the cylinder’s surface, with characteristic decay 
distance : 

Bx exp (==) forr < R. 
The boundary condition is that /) is continuous at the 
surface, and since ys 1 for superconductors, this is 
equivalent to By being continuous at the surface. Thus 
at r = R we require that B = Bo, e,, and this 
boundary condition is satisfied by the expression 

r 


B= Boz exp (—*) e, forr<R. 


(b) Using Ampére’s law, 


pod = curl B = curl [B,(r) ez] = — 


So 


(c) The current flowing in a ring of radius r, 
thickness dr, per unit length of cylinder is Jy(r) dr, so 
the total screening current per unit length flowing near 
the surface is given by 


R R 
Boz r—R 
J, ar=- f ex ( ) ar 
| 6 0 HOA LX 
me [o(S)), 
=— ex 
Ho " A 0 
me Poe (5) 
= 1 — exp = ; 
Lo » Lo 


ts 


since R > X. The minus sign in the result indicates 
that the surface current circulates around the z-axis in 
the —e,-direction. This circulating current produces a 
uniform field inside the cylinder, parallel to the axis 
— it’s equivalent to a solenoid with very closely 
wound turns. For a solenoid with n turns per unit 
length, carrying current J, the magnetic field inside 
the solenoid is gn e,, so replacing the current per 
unit length in this expression, nJ, by 7, for the 
superconducting cylinder, the field produced by the 
screening current is given by 


—Bo 
Bz = Lois = Uo X é == —Boz, 


which is equal in magnitude but opposite in direction 
to the applied field By = Boze,. 
Solution 9.6 The flux through unit area aligned 
perpendicular to the field is ® = BA = 0.50T m?. 
Since each core contains one flux quantum 
(2.07 x 101° T m2), the number of flux cores per 
unit area is 
0.50T m? 

2.07 x 10- Tm? 

that is, anumber density of 2.4 x 10!4m7?. 


=2.4~ 10!4, 


The cores are arranged on a hexagonal lattice 
(Figure 9.23 of Book 2). If the nearest neighbour 
separation is a, then the area of an equilateral triangle 
of side length a is 5a? sin 60°, and this is the area 
associated with one core. Thus 
; 6 1 

ay BO Sea 

so 


2 
= =9.8x107-'m. 
- laa x 10! m—? sin 60° ae 
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Solution 10.1 An invariant quantity has the same 
value in all inertial frames. Consider two frames, F 
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and F’, in standard configuration, with F’ travelling 
at speed v in the +--direction relative to frame F. Let 
the electric and magnetic fields be (E,, Ey, E,) and 
(B,, By, Bz) in frame F and (E%,, Ey, E,) and 
(B;,, Bi, Bz) in frame F’. Then using the field 
transformation equations (Equations 10.26—10.31 of 
Book 2), 

E! 2 CB’? 

= BE? +E? + Ei? — (BR? + BY + Bi?) 

= EB, +7 (Ej — 2vE,B, + vB?) 

+ (E? + 2uE,By + vB) 


UV v2 
= Ee +7 G + 2-5 By Bz + #2) 
2u vu 
+7? (2: = a BiEy + #2). 


The terms involving products of electric and magnetic 
field components all cancel, leaving 


BE’? — 2B"? = (BE? — 2B?) 


42h? (4 v 422 (4 e 


2.272 v? 2.272 v 
eB, (1-55) ex: (1-5) 


= B2+ BE? + E?-(B2+ B? +B?) 
= BE? ¢'B?. 


A similar analysis could be used whatever the relative 
velocity of the two frames of reference, so E? — c? B? 
has the same value in all inertial frames, that is, it is an 
invariant quantity — like the speed of light in free 
space, and electric charge. 


Solution 10.2. (a) The electric field E in 

frame F is (E, 0, 0) and the magnetic field B is 
(Bcos 0, Bsin@, 0). We use the field transformation 
equations (Equations 10.26—10.31 of Book 2) to find 
the fields in frame F’: 


i= c= HF, 
Ej, = (By — vBz) = 7(0 — v x 0) =0, 


E, = y(Ez + vBy) = 7(0+ vBsin 6) = yvBsin 8, 


B, = Bs =Beosd 


Bi = 7(By + vEz/c) = 7(Bsin 6 + 0) = yBsin8, 


BL = 7(B, — vE,/c?) = 7(0 — 0) =0. 


z 


So the electric field in ¥’ is E’ = (F, 0, yuBsin 6), 
and the magnetic field is B’ = (Bcos6, yBsin@, 0). 
The angle 6’ between the direction of B’ and the 
x-axis is given by 

yB sin 6 


tan 0’ = Beged. => ytan 0. 
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(b) For 6’ = 60°, we require 
tan 60° = ytan 30°, 
so y = tan 60°/ tan 30° = 3.0. 


Hence 
2 A 8 
1-5 = 5 and v= 5 C= 283x100 ms. 


Solution 10.3. An invariant quantity has the same 
value in all inertial frames. 


Rest mass is invariant; note that this mass is measured 
when the electron (or other particle) is at rest ina 
frame of reference. The relativistic mass, given by the 
ratio of the magnitude of the force acting on an 
electron to the magnitude of its acceleration, is not 
invariant but increases as the speed of the electron 
approaches the speed of light. 


Charge is an invariant quantity, so the charge of an 
electron is the same in all inertial frames of reference. 


The speed of an electron is not an invariant; the 
components of velocity transform as indicated in 
Equations 10.12—10.14 of Book 2. 


The energy of an electron is not an invariant; it is 
given by E = mc?/,/1 — U2/c2, where m is the rest 
mass of the electron and U is the electron’s speed in 
the frame of reference where the energy is measured. 


The charge density is not invariant because the 
volume occupied by a specific amount of charge will 
vary between different inertial frames, as indicated by 
the Lorentz transformation of spatial coordinates. 


The speed of light is invariant — this is one of the two 
basic postulates of special relativity. 


The quantity E - B is invariant — see Worked 
Example 10.2. 


Solution 10.4 (a) The positive charge is at rest, so 
there is no current density associated with this charge. 


The negative charge travels at speed v in the 
+<-direction. In one second, charge vp_ crosses 

unit area perpendicular to the +-direction, so 

J,— = vp_. This is the only component of the current 
density — there are no components in the y- or 
z-directions. 


(b) We will denote quantities in the frame that 
moves with the observer with primes. Then since this 
frame is in standard configuration with the frame 

in which the positive charge is at rest, the charge 

and current densities in it can be determined from 
Equations 10.22—10.25 of Book 2: 


2 
(p-- ar) 
= = Up 
p-— up 

2 

VU 
ee 
vo-( =). 


= p+ =-10-, 
Jy = V(Jx— — vp-) 
= y(vp— — vp_) 


Jy = V(Je+ — vps) 
= 7(0+ up_) = yJe-, 


i 0, 
Ty4 = Jy+ = ’ 
Je = J = 0, 
i = deel 


2 
So, net charge density = yp_ (: — =) —Yp- 
c 
2 
= ae 


and the net current density only has an x-component, 
given by 


J HJ. = lg = 0 


Solution 10.5 (a) Since charge Q is at rest in 
frame F, there is no magnetic force on it. The 
electric force is deduced by first using Gauss’s law to 
determine the electric field, which must be radially 
outward from the cylinder. For a Gaussian surface in 
the form of a cylinder with radius r, length /, coaxial 
with the cylinder of charge, the radial field at the 
surface has magnitude F’ given by 


2nrlE = pAl/e9, so E=pA/2neor. 


The force on charge Q at (0, d, 0) therefore has 
magnitude Frjece = QpA/2né0d, and is in the 
+y-direction. This is the only force on Q observed in 
frame F. 


(b) The charge density p’ and current density J’ in 
frame F’ are given by Equations 10.22-10.25 of 
Book 2. Thus 
v 

| (0 = Jz) =P; 
since the current density J, is zero in frame F. Also 

Jz = Wu — vp) 

= —~YUp; 


J, =1f= 0, d= =v. 


(c) The electric field in frame F’ is deduced in the 
same way as in part (a), but using p’ in place of p. 
Thus at distance r from the cylinder, 


E! = p' A/2neor = ypA/2re0r, 
and is directed radially outwards. 


The current J’ associated with current density J], is 

I’ = J!A; note that since dimensions transverse to the 
direction of motion are unchanged, A’ = A. At 
distance r from a long straight wire carrying current I’ 
the magnetic field strength is 


Bl = pol’ /2ar = po J). A/2ar = —poyvpA/2rr, 
with direction given by the right-hand grip rule. 


(d) The test charge Q is moving with speed v in the 
—x-direction in frame F’. The electric force on Q has 
magnitude 


Oe: _ QE’ = QypA/2ré0d, 
and is in the +y/-direction. 


The test charge moves perpendicularly to the 
azimuthal magnetic field, so the magnitude of the 
magnetic force at distance d from the cylinder is 


Frnag = Qv|B'| = Qu poyupA/2nd = Up Qyv* pA/2nd. 


Using the right-hand rule, with the velocity of Q in 
the —x-direction and the magnetic field circulating in 
the anticlockwise direction about the x-axis when 
viewed in the direction of increasing x, we deduce 
that this force is in the —y-direction. 


The net force in the +y-direction is therefore 


Bs ia aes 
_ Q7pA _ WoQyv?pA 
27Eod 2rd 


_ QYpA (1 2 

> 2rd E09 pen : 
Now c? = 1/eof10, 80 Uo = 1/egc?, and the 
expression for the force can be rewritten as 


pr — G04 (1 =) _ QA 
Qreod C2 Qreqyd 

(e) In both frames the force on charge Q is in 
the y-direction, so the equations for the z- and 
z-components are automatically satisfied and only 
Equation 10.18 of Book 2 is relevant. This equation 
indicates that the force in the frame where the particle 
is at rest is 7 times the force in the frame where the 
particle is moving at speed v. In this question, the 
particle is at rest in frame F, and so we expect 


fy = Lae Comparison of the net forces obtained in 
parts (a) and (b) show that they are related in this way. 


Note the entry for p. 227 of Book 2 in the list of errata 
for the course. 
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Solutions for Book 3 


Book 3 Chapter 1 


Solution 1.1 The general expression for the physical 
electric field of a plane wave is the real part of 


E = Epexp[i(k-r—wt+ ¢)|p, 


where Eo is the amplitude, k is the propagation 
vector, w is the angular frequency, ¢ is the phase shift 
and p is a unit vector in the polarization direction. In 
this case, 


k= Tyler ey +e2), 
k-r= (ety +2), 
w = 2rf, 
B= (ere), 
o=9, 
and so 


E= 7 ae c (je a 2nft) | (eoe), 


The physical electric field is the real part of this 
expression, that is 


Eo k 
E = — cos | —~=(x#+ y+ z) — 27 ft }(e,—e,). 
phys V2 (= 7] ) f ) ( x y) 
(b) Fora transverse wave, the electric field is 
perpendicular to the propagation direction, so 


~ 


E-k = 0. For the wave in part (a), 


a 1 
and k=-—~(e,+e,+e,), 


E=p= Va (€2—ey) V3 


a 1 

E-k= ve (2 ~ ev) -(e, +e, +e,) = 0, 
and therefore the field is transverse to the propagation 
direction. 


(c) Any polarization that satisfies p - k =Oisa 
possible solution. Examples are p = (—e, + e,)/V2 
and p = (e, — 3e, + 2e,)/V14. 


Solution 1.2 A general expression for the electric 
field of a plane monochromatic electromagnetic wave 
is 


E = Eocos(k-rtwt+ ¢)p, 


where Ep is the amplitude, |k| is the wavenumber, 
w = 27f is the angular frequency, ¢ is the phase 
shift and p is a unit vector in the polarization 
direction. Now k- r = kx + kyy + kzz, and 
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|k| = ,/k2 + k2 + k2, so the wavenumber is 


k = »/(38)? + (28)? + (8)? 
= /146 


= V14 x 6.0 x 10°m7! 


=2.2~x10'm!. 
w=Imf=y7, 
so f= /2a 
=6.7x 10" s-1/2n 
=1.1 x 10 Hz. 


The amplitude Eo of the wave is given by 
Eo = al(ez — 3e,)| 


=2.0x 10?Vm7! x 12 + 32 


=6.3x10?Vm_!. 


The phase speed of this wave is 
we 67x 10's 
ko 2.2x107m-! 
which is consistent with the value of the speed of light 
in vacuum, c = 3.00 x 108 ms~!. 

Solution 1.3 Since B = k x E/c, the direction of 
B is the direction of k X E, which is the direction of 
(2e, + 2e, —e,) X (ey + 2e,) = Se, — 4ey + 2e,. 
Expressing the direction as a unit vector, the magnetic 
field is in the direction 


(Se, — de, + 2e,)/V 45. 


os = 3.0 x 108 ms}, 


Solution 1.4 You may have recognized the types of 
some of these waves by inspecting the form of the 
expressions given and comparing them with the 
original wave. If so, a simple statement of the type of 
wave would be sufficient. However, below we include 
some explanation to justify the answers. 

(a) The plane wave Eo expli(kx — wt)] ey is 
identical to the original wave, that is, it has the same 
amplitude Eo, wavenumber k, angular frequency w, 
phase shift (zero) and polarization e, and travels in 
the +<-direction. The physical electric field of the 
combination is the wave 


Re} 2Eo exp|i(ka—wt)] ey } = 2E) cos(kx—wt) ey, 
which has twice the amplitude of the original wave, 


but otherwise is the same. 


(b) The wave Eo exp[i(ka — wt)| e, differs from the 
original wave in that it is polarized in the z-direction, 
at 90° to the polarization of the original. The 
combination is a linearly polarized wave, 


Ref Fe exp|i(ka — wt)] (ey + e.)} 


= Eo cos(kx — wt) (ey + ez), 


which has amplitude /2 Ep and is polarized in the 
direction (e, + e,) //2 (see Equation 1.25 and Figure 
1.11 in Book 3). 


(c) This case is similar to (b), but the additional 
wave has twice the amplitude of the original wave. 
The combination is again a linearly polarized wave, 


Re Eo exp|i(ka — wt)] (ey + 20.)} 


= Eo cos(kx — wt) (ey + 2e,), 


which has amplitude 5p and is polarized in the 
direction (e, + 2e.)/V5. 


(d) The wave Ep exp|i(kx — wt + 7/2)] e, is 
polarized at 90° to the original wave and its phase its 
advanced by 7/2. The physical electric field is 


Ref Eo expli(kx — wt)| ey 


+ Eo expli(kx — wt + 1/2)] e.} 


= Eo|cos(kx — wt) ey + cos(kx — wt + 7/2) e,] 
= Eo|cos(kx — wt) ey — sin(kx — wt) ez], 


which is a circularly polarized wave (see Subsection 
1.3.5 and Figures 1.12 and 1.13). 


(e) This case is similar to (d) except that the phase 
of the wave Ep expli(ka — wt — 1/2)] e, is retarded 
by 7/2 relative to the original rather than being 
advanced by 7/2. The resultant field is 


Eo|cos(kx — wt) ey + sin(kx — wt) ez], 


which is a circularly polarized wave with the electric 
field vector rotating in the opposite sense to case (d). 


(f) The wave 2Eo exp|i(kx — wt + 7/2)]e, has 
twice the amplitude of the wave in (d), but otherwise 
is the same. The combination of this wave with the 
original gives a physical electric field 


Eo|cos(ka — wt) ey — 2sin(kx — wt) e,]. 


The polarization direction rotates in the same way as 
in (d), but the electric field is twice as large when the 
polarization is in the direction e, as when it is in the 
direction e,; this is an elliptically polarized wave, 
rather than a circularly polarized wave. 


(g) The wave Eo exp|i(ka — wt — 7)] ey is identical 
to the original except that its phase is shifted by 7. At 
any position and time, the fields of the two waves are 
equal in magnitude but opposite in sign, and so the 
two waves cancel — they destructively interfere. 


(h) The wave Eo exp|i(ka + wt)] ey is identical to 
the original except that it travels in the opposite 
direction. The combination of the two waves leads to 


the field 
Ref Fe (exp[i(ka — wt)] + expfi(ka + wt)]) ey} 


= Re{ Eo exp[ikz] (exp[—iwt] + exp[+iw?]) ey} 
= 2Eo cos kx cos wt ey. 


This is a standing wave, not a travelling wave. Note 
that for kx = mz, where m is any integer, the electric 
field is always zero. Standing waves are discussed in 
Chapter 5 of Book 3 in the context of reflection of 
electromagnetic waves from conducting surfaces. 


(i) In this case the two waves differ by 1% in 
frequency and wavenumber. Their combined field is 


Ref Fe (exp[i(ka — wt) 


+ exp[i(1.01ka — 1.01wt)]) } ey 


= Eo(cos(ka — wt) + cos(1.01ka — 1.01wt)) ey 
= 2Eo cos(1.005kax — 1.005wt) x 
cos(0.005ka — 0.005wt) ey, 


where we have used the trigonometric identity 


4) (4) 
cos 
2 2 


to combine the two cosine terms. This is a wave 
that has approximately the same frequency and 
amplitude as the original (0.5% greater actually), but 
its amplitude is modulated by a more slowly varying 
sinusoidal waveform, with wavenumber 0.005 and 
angular frequency 0.005w. This is an example of the 
phenomenon of beats, which may be familiar in an 
audio context. 


cos A + 00s B= 200s ( 


Solution 1.5 (a) We first use the wavelength to 
determine w and the components of k. 


w= arf 
_ 277C 
—s 
_ 2 83.00% 10 ms 
7 527 x 10-9 m 
= 3.6 x 10'°s7}. 

27 _ 21 

\ 527 x 10-9m 


Since the pulse is travelling in the +-y-direction, 
y= 12% 10m and h, = hy = 0, 


k= =12~x10'm"!. 


To find the fields, we first deduce the power crossing 
unit area, and hence the time-average of the magnitude 
of the Poynting vector. This is directly related to the 
amplitude of the electric (and magnetic) fields. 


The mean power in the pulse is 
300J/10-° s = 3.0 x 10!! W, so the power crossing 
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unit area, which is the time-averaged value of the 
Poynting vector, N, is 
3.0 x 10’ W 


N= = 3.8x10Wm”. 
mx (0.5 x 10-4 m)2 . _ 


From Equation 1.34 of Book 3, N = se0bGC, 


NiO) 
l2N 
Eo — —— _ §_ 
Ege 
2x 3.8 x 10!9W m2 
8.85 x 10-12C2N7! m-2 x 3.00 x 108 ms! 


=1.7x10"Vm"!. 


Since the polarization is in the x-direction, 

Foe = 17 8 10"'V m-* and Boy, = Zoe = 0. 

The magnetic field is related to the electric field by 

B=k xX E/c (Equation 1.16), so 
1.7% 10" Vm" 


B= (ey Xe) 3.00 x 108ms-! 


Thus Boz = —5.7 x 10? T, and Box = Bo, = 0. 


(b) The peak value of the electric field is almost six 
orders of magnitude greater than required to ionize air. 
So the energy of the pulse would be absorbed rapidly 
in air, which is why beams with such high intensities 
are made to travel through a vacuum. When powerful 
beams must propagate in air, they are expanded with 
lenses, so that the Poynting vector has a much smaller 
magnitude. 


5.7x10?T e,. 


Solution 1.6 A polaroid filter transmits the 
component of the electric field that is parallel to the 
axis of polarization and blocks the perpendicular 
component. The vertical electric field of the incident 
wave can be resolved into component fields 

Eo cos 30° parallel to the axis of polarization of 

the first filter and Ep sin 30° perpendicular to the 
axis of polarization of this filter. Only the parallel 
component is transmitted, so the amplitude of the 
wave incident on the second filter is Ep cos 30°. 
This wave, which is polarized at 30° to the 

vertical, can be resolved into components parallel 
and perpendicular to the axis of polarization of 

the second filter. In case (a) the axis is vertical, 

and the component of the wave in that direction is 
(Eo cos 30°) x cos 30° = Ep cos? 30° = 3Ep/4, 

so this is the amplitude of the transmitted wave. 

In case (b) the axis of polarization of the second 
filter is horizontal, and the component of the wave 
in that direction is (Ep cos 30°) x sin 30° = 

Eo sin 30° cos 30° = 3/4, so this is the amplitude 
of the transmitted wave. 

Note that even though the initial wave was vertically 
polarized, after passing through the first filter the 
wave has components of polarization in both vertical 
and horizontal directions. 
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Solution 1.7. (a) For a monochromatic plane wave, 
Bo = Eo/c, 


Nie) 


W — 12,72, — 12, 772,3 — 1p? 
N = 5€0bgc = 5€0Boe’ = 5 Boc/ uo, 


where we have made use of the relationship 

a = 1 re EoLOo. 

(b) The magnitude of the Poynting vector is equal to 
the power per unit area, 


SO 
= 1/2 
2N 
By = ( He) 
Cc 
(2 x 10-19 W m-? x 4m x 10-7 NA-2\? 
7 3.00 x 108 ms~2 
=9.2x 10-* T. 
Book 3 Chapter 2 


Solution 2.1 For a Hertzian dipole, the dipole length 
dl is much less than the wavelength \ of the radiation 
it produces (Section 2.2 of Book 3). For the field to be 
well-approximated by the radiation terms alone, we 
require r >> /2m (Equation 2.18). 


@) A=c/f =3% 10° ms 50s =6 & 10° m, 
which is much greater than the dipole length, 10 km. 
So this is a Hertzian dipole. The field is a radiation 
field for distance much greater than 6 x 10° m/2r, 
that is beyond 10’ m, say. 


(b) This dipole radiates VHF radiowaves, with 

= 3x 108ms!/108s~! = 3m, which is shorter 
than the dipole length, 10 m. So this is not a Hertzian 
dipole. 


(c) This is a microwave source, with 

\=3x 108ms—!/3 x 10°s~! = 0.1m, which is 
much greater than the dipole length, 1 cm. So this is a 
Hertzian dipole. The field is a radiation field for 
distance much greater than 0.1 m/2z7, that is beyond 
20cm, say. 


Solution 2.2. When r is small the terms involving 
r—3 will be large compared with the terms involving 
r+ and r~?, so the electric field in Equation 2.17 
reduces to the expression that is printed below it in 
Book 3: 


2 
E = <" exp [i(kr — wt)] ( 


—iw 


2cos 0 sin 0 

2 ep + 73 es) , 
Now 177) = polo 61/47, and the strength of the current 
dipole Ip 61 is related to the strength of the charge 
dipole po by Ig 61 = wpo. Substituting these into the 
pre-exponential term in the equation for E we obtain 


cn _ coo 1 - c* Wowpo = c? Wopo 
—iw —i4tw —i47w An 


Then since c? 


becomes 
ipo 
Ate 0 ; 


= 1/eo,UW0, this pre-exponential term 


In the limit w — 0, the wavenumber k = w/c — 0, 
and so the exponential term tends to unity. The 
physical electric field then becomes 


_ ipo (2cosé sin 6 
Ephys = Re ro ( 73 e, + e ea) } 


D0 2c0s 6 | eee 
Ane rr i i 
which is the same as Equation 2.1 for the electrostatic 
field of a charge dipole. 


Solution 2.3. (a) The wavelength of the radiation is 
A=6/f =3 x 10°ms—'/3 * 10? s-* = 10cm, s0 

r >> /27, which means that the fields are given by 
the expressions for the radiation fields in Equations 
2.19 and 2.20. 


The electric field produced by the Hertzian dipole is 
always in the eg direction (though this direction 
depends on the location of the point where the field is 
measured). The magnitude of the amplitude is 
w7 sin 6/r so 

[amplitude at (2.0m, 7/4, 7/2)] 

__ (sin(a/4)/2.0m) 

~ (sin(a/2)/1.0m) 

1/(2.0V2 
2 OW?) 5646 Vin 
1.0m-! 
=71 x10" Vi“. 


[amplitude at (1.0m, 7/2, 0)] 


So Ep = 7.1 x 107-7 Vm ' eg. 


(b) The magnetic field at (2.0m, 7/4, 7/2) is in the 
ey direction. The magnitude of its amplitude is just 
the magnitude of the electric field amplitude divided 
by c, that is, 


7.1%107! ¥ m-*/3.00%10° ms = 2.4x10- T. 


So By = 2.4 x 10-% Teg. (Remember that the 
directions of the unit vectors e,, €9, &¢ in spherical 
coordinates depend on the location of the point that is 
being considered.) 
Solution 2.4 
From Equation 2.24 of Book 3, the time-averaged 
Poynting vector is given by 

wn? 
Qcpor? 


N= sin? Oe... 

The distance r from the origin to the point 

(2.0, —2.0, 1.0) mis 2.02 + 2.02 + 1.02 = 3.0m. 
The angle @ is the angle between the direction 

of the dipole and the direction of the position 
vector of the measurement point, and since the 


dipole is aligned with the y-axis, this is given by 
cos 6 = y/r = —2.0/3.0, so 6 = 132°. 
Also, 7 = polo 61/47. Thus 
N- w? ju (Io 51)? 
3212cr? 
(2.0 x 10°s—!)?x4a x 10-7NA~?x(3.0 x 10-8 Am)? 


sin? 132° e, = 


g277 X3.00 x 108 me! x (3.0 m)* 
x sin? 132° e, 
=29x10 Wm “e,, 


where the unit vector e, that indicates the direction of 
power flow is in the direction of the position vector of 
the measurement point, that is, 


e, = (2.0e, — 2.de, + 1.0e,)/3.0. 


Solution 2.5 


The power per unit area transmitted as radiation 
passes through a scattering medium is given by 


N = Noexp[—noz] = Noexp [—z/Lscat], 


where Np is the incident power per unit area, n is the 
number density of air molecules, o is the mean 
scattering cross-section for air molecules and Lecat 
is the scattering length. The scattered power is 

(No —N). 

The number density is given by the ideal gas law 
(Equation 2.33): 


_ 1.0 x 10° Pa 
~ 1.38 x 10-233 K7! x 273K 
= 2.65 x 107° m-?. 


Thus the transmitted power after 1.0 km is 
N= Nox 
exp [—2.65 x 107° m73 x 1.50 x 1073? m? x 10° m] 
= 0.961No. 
Therefore the percentage of the incident power that is 
scattered is 

(1 — 0.961) x 100% = 3.9%. 


The scattering cross-section is proportional to w* and 
therefore proportional to \~*. (We ignore the weak 
dependence of the polarizability on w.) For red light 
with \ = 600 nm, 


4 
= Aplue 
Ored = Pblue | \— 
Ared 


400 \ + 
= 1.50 x 10729 m? ( — 
- ™ (=) 


= 2.96 x 1073! m?, 
Hence 


N= 
No exp [—2.65 x 10°? m7? x 2.96 x 1073! m? x 108 m] 
= 0.992Np. 


47 


Therefore the percentage of the incident power that is 
scattered is 


(1 — 0.992) x 100% = 0.8%, 
which is five times smaller than for blue light. 
Solution 2.6 


Light from the sky that reaches you after travelling 
perpendicular to the direction of the incident sunlight 
is most strongly polarized, with the polarization 
perpendicular to the plane containing the Sun, the 
region of the sky that you are observing and you (see 
Figure 2.16 in Book 3). When the Sun is directly 
overhead, with direct sunlight incident in the vertical 
direction, the light reaching you from the horizon 

in any direction will be strongly polarized in the 
horizontal direction. Light from the sky in directions 
close to the Sun will be unpolarized. Thus if you are 
holding the polarizing filter so that it blocks light that 
is polarized in the north-south direction, then the 
horizon will appear bright to the east and to the west, 
but the sky to the north and south will appear dark. 
Conversely, if the polarizing filter is rotated so that 

it blocks light that is polarized in the east-west 
direction, then the horizon will appear bright to the 
north and to the south, but the sky to the east and 
west will appear dark. However, a clearer way of 
describing what you would observe is to say that when 
you look towards the horizon in any direction, the 
light will be polarized in the horizontal direction. 
When you look close to the direction of the sunshade 
overhead, the light is unpolarized. 


Solution 2.7 (a) The time-average power per unit 
area radiated by a Hertzian dipole is given by the 
magnitude of the Poynting vector (Equation 2.24 of 
Book 3): 

2,2 


= +2 
= ear sin’ 0. 


The ratio of the magnitudes of the Poynting vectors at 

the points (50m, 90°, 0) and (100m, 60°, 90°) is 
sin? 90°/(50m)?_—:16 
sin? 60°/(100m)2_— 3 

Note that the magnitude of the Poynting vector 


is axially symmetric — it doesn’t depend on the 
o-coordinate. 


= 0.3. 


(b) The average power radiated in the range 
70° < @ < 110° is given by (see Subsection 2.4.3) 


110° 27 22:2. 
w= i wT sin? @ r2 sind. dddé 


~ 2cpor? 
Qn 110° 
ae a: a sin? 6 d@ 
CHO 70° 
wn? 110° 
= (27) E cos? @ — cos q 
2ChLo0 70° 
2,2 
= 2 (om) 0.657. 
CHO 
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The total power radiated is determined in a similar 
way, except the # integral runs from zero to 180° (or 
m), and this replaces the factor of 0.657 by 1.33 (as 
shown in Subsection 2.4.3). So the fraction of the total 
power radiated in the range 70° < 6 < 110° is 
0.6577 Lae = 0.49, 


Book 3 Chapter 3 


Solution 3.1. (a) f =w/2m7 = 1.0 x 10° Hz. The 
frequency does not change when the radiation travels 
through different materials. 

(b) Phase speed is given by v = c/n, where 

the refractive index n of the material is the 


square root of the relative permittivity «. Thus 
v = 3.00 x 108 ms—!//4.0 = 1.5 x 108ms"! 


(i) Ajay =2/f = 
15% 10° ms" /10% 10 heim, 


(d) Wavenumber k = 27/\ = 27/1.5m=4.2m7! 


Solution 3.2 (a) The fields for the incident, 
reflected and transmitted waves are 


E; = Eyexpli(hy —wt)Jez, (given) 
R= “ expli (ky — wt)] ex, 

E, = rn Eo exp|i(—ky — wt)] ez, 

B,; = ns exp[i(—ky — wt)] ex, 


E; = ty Ep exp[i(nky — wt)] ez, 


tn Vek 
By = tnveEo expli(nky — wt)]e 
c 
In writing down these expressions we have used the 
relationships 


Eo nEo JeEo 
Bo ’ 
U Cc Cc 
Ee =TEqo and Ei = typ. 


We assume that all of the electric fields are in the 
same direction, and determined the directions of the 
magnetic fields from the relationship B = k X E. 


(b) The boundary conditions are that £) and H) are 
the same on both sides of the boundary. Since jz = 1, 
the second condition is equivalent to B being the 


same on both sides of the boundary. Thus 
E,+E,=E, so 14+7,= th, 


and 


1—ra= Vet: 
Solving these simultaneous equations gives 
1l— Ve 2 
ve and t, = ————~. 
1+ jé 1+ je 


These expressions are equivalent to those in Equation 
3.14 of Book 3 if we set ny; = 1 and ng = Ve. 


Tr = 


Solution 3.3. (a) The refractive index of a material 
isn = c/v = (eu)'/2. For dielectric materials, p = 1 
generally, son = el/2 — 1,71/2 = 1.304, or 1.3 to 
two significant figures. 


For normal incidence in air, the expression for the 
reflectance in Equation 3.17 of Book 3 becomes 


l—n\*  /1-1.304\? 
R= ~) =(——-—) =0017. 
l+n 1+ 1.304 
(b) The amplitudes of the reflected and transmitted 


waves are obtained from the amplitude reflection 
ratios (Equations 3.14): 


y 
= 5.0Vm7! | ——__ 
1+ 1.304 


The negative amplitude for Eo; indicates that the 
reflected wave is 7 out of phase with respect to the 
incident wave. 


The boundary condition at the interface is 

E, + E, = E;, which in this case reduces to 

Fo, + Eor = Loz. Substituting numerical values for 
the left-hand side, 


5.0Vm!—0.66Vm ! =434Vm 1), 


which is equal to the calculated amplitude of the 
transmitted wave, as required. 


Solution 3.4 The transmittance — the fraction 

of the power transmitted across a boundary — is 
independent of the direction that radiation travels 
across the boundary (see Equation 3.17 of Book 3). 
So the solutions for part (a) and part (b) are the same. 
For light travelling from air to glass, ny = 1.00 and 
n2 = 1.50, so the transmittance of the air—glass 
boundary is 


Anin2 
(n1 +n)? 

4x 1.00 x 1.50 
~ “(1.00 + 1.50)? 
= 0.960. 


Tair glass — 


For the glass—water boundary, 


4x 1.50 x 1.33 


—s a" = 0,996. 
(1.50 + 1.33)? 


Tslass—water = 


So the transmittance through the glass wall from air 
into water is 


Telass wall = 0.960 x 0.996 = 0.957, 
or 0.96 to two significant figures. 


Note that the transmittance tends to unity as the 
difference between the refractive indices on either side 
of the boundary tends to zero. 


Solution 3.5 The incident beam, which is polarized 
at 45° to the scattering plane, can be resolved into two 
components, one normal to the scattering plane and 
one in the scattering plane, with equal amplitudes and 
equal powers. For each component we determine the 
reflectance using the expressions in Equations 3.29 
and 3.33 of Book 3, which require knowledge of 6. 
From Snell’s law, 


6, = sin~1(sin 45°/2.42) = 16.7°. 
Thus 


2 
m1 cos 0 — nz cos 6; 
Risp — 
ny 1 cos & + nz cos 6; 


_ (cos 16.7° — 2.42 cos 45° \* 
~ \cos 16.7° + 2.42 cos 45° 
= 0.080. 


2 
R ny cos 6; — nz cos 6 
ns — 
: ny cos 6; + n2 cos 


_ (cos 45° — 2.42 cos 16.7° 2 
cos 45° + 2.42 cos 16.7° 
= 0.283. 


Since the two polarization components carry equal 
powers, the overall reflectance is simply the average 
of the isp and nsp values, that is, R = 0.18. 


Solution 3.6 Circularly polarized light can be 
regarded as the sum of two linearly polarized waves, 
with a phase shift of 7/2 between them, and with their 
polarizations in any pair of orthogonal directions in 
the plane perpendicular to the propagation direction. 
We can choose these directions to be the isp and nsp 
directions. For light incident at the Brewster angle, the 
reflectance is zero for the isp component, so the 
reflected wave is linearly polarized normal to the 
scattering plane. 


Of course, this result is not restricted to circularly 
polarized incident light. Any beam reflected from a 
surface at the Brewster angle will be polarized normal 
to the scattering plane. 


Solution 3.7. (a) The beam will be reflected back 
and forth in the plate if its angle of incidence to the 
interfaces is greater than the critical angle, @ri, which 
is given by 
sin Oerit = N2/M1, 
so that 
Gert = sin 1(1.50/1.55) = 75.4°. 
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So for angles of incidence between 75.4° and 90°, the 
beam will be confined within the plate (almost!). 


(b) For angles of incidence greater than the 

critical angle, the electric field is evanescent in 

the blocks either side of the plate, as discussed 

in Subsection 3.4.5 of Book 3. The field decays 
exponentially with distance z from the boundary, as 
exp[—z/6], where 6~' = kt| cos |. The wavenumber 
in the blocks is given by 


Qrn _ 27 x 1.50 
Ao «1500 x 10-9m 


ky = = 6.28 x 10°m™!. 


cos 6, = 4/ (1 — sin? 4;) 


= V¥—0.0356 = 10.189. 
So the typical penetration depth 6 for the field is 


6 = (k| cos &|)~4 
= (6.28 x 10°m™! x 0.189)~! 
= 84x 107m. 


This is approximately half of the free space 
wavelength (1.5 x 10~®m). 
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Solution 4.1 Since the wave is travelling in the 
e,-direction, the amplitude of the magnetic (and 
electric) field will decrease exponentially with the 
y-coordinate: 


Bo « exp[—Kimag y] = exp | ee | : 


So the ratio of the amplitudes at the two points is 


Nimag¥ 
exp | — Y2 
= Nimagv 
= exp (yo — y1) 
imagW 
exp | j 
Solution 4.2 
Kreal = @2, Kimag =a, W= a3, 
age Cc a,c 
Nreal = — Kreal = —, Nimag = — Kimag = 
WwW a3 Ww a3 
Since e = n? = n2,,, — n2__, + 2i MrealNimae, then 
real imag real! "mag > 


co 2 
Ereal = Nreal — “imag 
_ (a5 — az)c? 


2 o] 
a3 
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Eimag = 2Nreal Nimag 


2a1a9C? 


2 ’ 
a3 
Uphase = C/ Meal 
= a3/a2, 
absorption length = 1/kimag = 1/a1. 


Solution 4.3. (a) (i) Refractive index 
N= Nyeal + i Nimag = Ve, and when €;ea) > Eimag> 
then (as shown in Exercise 4.2 in Book 3) 


Nreal = VW Ereal = V 1.46 = 1.2, 
is _ Eimag == 1.0 x 19710 
— ay Ereal 2x V1.46 


(ii) The wavenumber is complex and is given by 


=41x107!!. 


k= Kyeal +i Kimag 


W é 
= ~ (real ae Nias) 
21 . 
— Jo (zea! +1 ina) 
_ 27 
~ 2.0 x 10-§m 
= 3.8 x 10°m-+4+i1.3 x 10-*m?. 


(1.21 +141 x 107") 


(iii) The phase speed is real and is given by 
Uphase = Ci ipasl 
= 3.00 x 108ms~1/1.2 
=25% 10? ms“; 
(iv) The amplitude falls exponentially with distance, 


as exp|—kimagy], So over a distance of 1.0km the 
amplitude falls by a factor 


exp[—1.3 x 10-4m™! x 1.0 x 10? m] = 0.88, 
which corresponds to a decrease of about 12%. 


(b) The physical electric field of a plane wave 
travelling in the x-direction and polarized in the 
y-direction, with amplitude Fo at the origin is 


Esty = Re{ Ep expli( Meat 1 iinet — wt)] } ey 

— Re{ Eo exp|—Kimae2| exp|i(Krea® — wt)]} ey 

= Eo exp|—Kimag2] cos(Krealt — wt) ey. 
Substituting the values of yea) and kKimag from part 
(a), and with w = 2mc/Ag = 9.4 x 10!4 s—!, we obtain 

Epnys = Eo exp[—1.3 x 10-4 m7 ta] x 

cos(3.8 x 10°m7*z — 9.4 x 104 st) ey. 

We determine the magnetic field using the same 
method as in Worked Example 4.2. Since 


Ephys X Bpnys must be in the propagation direction, 
e,, we deduce that Bpnys must be in the e,-direction. 


Also, the amplitude of the magnetic field is (n/c) Ep, 
so 
(Nreal +i Tima) 
€ 
exp[i(krealt — wt)] ex. 


B= 


Fo exp[—Kimag?] x 


The physical magnetic field is the real part of this, that 
is 


Fo 
Bpiye = ma exp|—Kimag 2] X 


(Tinea cos Cae a = wt) — Nimag sin ies _ wt)) e,. 


This can be rewritten as 


Eo 
Bogs = oe exp|—Kimag] @ COS(Krealt—wt+/3) ez, 


where 

_ 2 2 _ : ; 
OQ = 4/4) + Nimag = Mreal, SINCE Nreal > Nimag; 
and 


8 = tan”! (Mimag/Nreal) 
=tan™* (4.1 x 107/12) 
=tan (4x 10"), 
which is small enough to be negligible. Then since 


Nreal/C = 1/Vphase, the magnetic field is 


B 


phys => : exp|—Aimas®] cos(Kpeait — wt) ez 
phase 
x 
rae os exp[—1.3 x 1074 m1] x 


cos(3.8 x 10°m~1z — 9.4 x 10'*s~t) ex. 


Solution 4.4 The simple classical model predicts a 
complex dielectric function with real and imaginary 
parts given by Equations 4.8 and 4.9 of Book 3: 


2 2 
7 ¥ (wi — w*) 
Ereal (w) =1+ Wy (w2 _ 2p A wry?’ 
wy 


Eimag(W) = we ie 
n 


(a) Whenw — 0, 
Ww 
Ereal > 1+ 4 =1+ 2.00? = 5.0, 
WwW. 
nh 


Eimag — 0. 


(b) When w — oc, 
2 2 


wey 
and €imag — r 


Ereal 2 1-2 1, 3 
WwW Ww 


2 


(C) rea) = 1 when w = wy = 4.0 x 165-1, 
(J) Atw =w, = 40 x 10s", 
fas = w? /wn’y = 2.00? = 4.0. 


Solution 4.5 (a) We ignore the possibly 
complicated geometry of the steak, and treat it as a 


semi-infinite slab. The microwave amplitude within 
the steak will then fall exponentially with distance z 
below the surface, as exp[—kimag2], and the depth d at 
which the amplitude has dropped to half of the surface 


value is given by 
exp|—Kimagd] = 55 or d=1n2/kimag. 


Now Kimas = Rimagw/c, and 


2 2 
—Ereal + VV Sica Cimag 


Nimag = 


2 
_ —40 + V40? + 12? 
7 2 
= 0.8806. 
So 
= cln2 
Nimagv 


3.00 x 108ms~" x In2 
~ 0.8806 x 2m x 2.45 x 109 Hz 
= 1.5cm. 


(b) For polystyrene, cimag < €real and we can use 
the approximation 


Eimag 
2V/ Ereal 
3x 


2/1.03 
= 1.478 x 107°. 


Nimag ~ 


So 


is cln2 


Nimag¥ 
_ 3.00 x 108ms—! x In2 
~ 1.478 x 10-5 x Qa x 2.45 x 109 Hz 
=9.1 x 107m. 


So the microwave amplitude will hardly be affected by 
even the thickest polystyrene containers that are used 
in a microwave oven. In other words, the microwaves 
pass through polystyrene with minimal absorption. 


Solution 4.6 (a) From Snell’s law, 

sin # = sin 6;/n, so because blue light is more 
strongly refracted than red, then npiue > Nred- 
Therefore, since blue light has a higher frequency than 
red light, dn/dw > 0. 


(b) The group speed is given by (Equation 4.27 of 
Book 3) 

—— 

wdn/dw +n’ 

so since dn/dw is positive, Ugroup < c/n. But the 


phase speed is c/n, so for light in glass the group 
speed is less than the phase speed. 


Ugroup = 


However, since n ~ 1.5 for glass, both speeds are 
much less than the speed of light in a vacuum. 


51 


(c) Since Vgroup < Uphase for light in glass, the 
individual maxima and minima of the electric field 
will travel faster than the envelope of the pulse, so 
they will progress from the rear to the front of the 
pulse. 


Solution 4.7 (a) 
4.59 x 107) m? 


= 1.522 = 1.551 

Mblue 522 + (4 1-2 dol, 
4.59 x 107} m? 

Mred = 1.522 + (6x 10-7m)y2 = 1.535. 


(b) Phase speed is given by v = c/n, so 
Vbue = 2.998 x 10° ms—1/1.551 
= 1.933 x 10° ms7!, 
red = 2.998 x 108 ms! /1.535 
= 1.953 x 10°ms7!. 
Note that we have used a more precise value for c 


because the refractive index data is given to four 


significant figures. 
(c) The group speed for waves with frequencies 
around w is given by Equation 4.27 of Book 3: 
c 
wdn/dw +n 


Ugroup = 


Now 


B w2B 
m ar r e An?’ 
and so 


dn 2wB 


dw = 412¢2’ 


and 


Ugroup = 


Hence 
2.998 x 108 ms~! 
3 x 4.59 x 107) m? 
(4 x 10-* m)? 
= 1.864 x 10°ms“!, 
2.998 x 10° ms7! 
3 x 4.59 x 107) m? 
(6 x 10-7 m)? 
= 1,921 « 10° ms™, 


Ugroup blue = 


1.522 + 


Ugroup red = 


1.522 + 
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(d) The transmission angle is given by Snell’s law, 
ny, sin 6; = no sin &. With ny = 1.000 for air, and 
6 = 30°, 
A; blue = sin~} (sin 30° /npiue) 
= sin71(0.5/1.551) = 18.81°, 
req = Sin} (sin 30° /Nrea) 
gin 0541595) = 19.01". 


So the blue light is refracted 0.2° more than the red 
light, a 2% difference in the angle of deviation 


(6; — 4). 


Solution 4.8 The absorption length is the 
characteristic distance associated with the exponential 
decay of field amplitude in the factor exp [—kimag 2], 
that is, 


1 


absorption length = 
imag 


Cc Xo 


Nimag¥ 27 Nimag 


The values of Nimag for water and silica at the three 
frequencies can be estimated from Figures 4.15 and 
4.17 of Book 3, and are tabulated below. 


500 nm 1000 nm 2000 nm 


Nimag for water 10-8 ox l0. 2x 10-? 


Nimag for silica 1971 10-12 10710 


The absorption lengths calculated from these data are 
as follows 


500 nm 1000 nm 2000 nm 


absorpn length 8 5x10" 210 


for water /m 


absorpn length 8 x 10° 2x 10° 3 x 10° 
for silica /m 


Note that the absorption length is three to seven 
orders of magnitude longer in silica than in water. 


Solution 4.9 We require that the separation of the 
two wavelength components of each pulse is less than 
10% of the 10 pus interpulse separation, that is, less 
than 1 jus. The time taken for the component centred 
on the free space wavelength ., to travel distance D is 


D D B 
=—(A+S), 
Ugroup c At 


and the time for the component centred on gz to travel 
the same distance is 


D B 
—{A+—}. 
> (4458) 


The difference between these times is 


ol(4+m)- (G+) 


_ D 
~ 3.00 x 108ms-! 


x [1.062 x 107" m?] 


1 1 
x 
a x 10-"im)* (Gol. 10-? =a] 
= Dx 1450s 10- -*s: 
This difference must be less than 1 jus, so 


Dx1.455x107'' mts < 1x10-°s, or D < 69km. 
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Solution 5.1 (a) wp, = \/ne?/meép is the plasma 
frequency, which is the frequency of collective 
oscillations of the free electrons in the conductor. 7 is 
the collision time for the free electrons. 


(b) e€(w) is real when the imaginary term in the 
denominator is negligible, that is, when w >> 1/7. In 
this frequency range, e(w) = 1 — w/w. 


(c) For very high frequencies, that is, when w >> uw, 
the relative permittivity function is unity. This 
corresponds to the electrons, permanent dipoles, etc, 
being unable to respond to the rapid changes of the 
electric field because of their inertia. 


(d) In the range 1/7 < w < wp, the relative 
permittivity is real (see part (b)), and it is negative. 
This means that the refractive index (n = \/e) and 
the wavenumber are imaginary, which means that 
electromagnetic disturbances are evanescent, rather 
than travelling waves. 


(e) Forw < 1/t%, we obtain 


Hy be 
w twist 
ew) =1—-—® =14—+. 
iw/Te w 
Now ws = ne?/meo, and o = ne*7,/m, so in the low 
frequency limit, 
ig 
é(w) = 14+ — 
WEQ 


Solution 5.2 When a wave meets the 

surface of a dielectric, it is partly reflected 

and partly transmitted into the material. 

The transmitted wave has wavenumber 

k= Reet > i Aimae = (esa 1 mae) Git travels 
at speed v = C/Nyeal, and its amplitude decays 
exponentially with distance as exp|—kimagz]. For a 
material like glass the reflection is small, and the 
attenuation of the wave is also small. The electric and 
magnetic fields of the transmitted plane waves are 
transverse to the direction of propagation, orthogonal 
to each other, and they are in phase (as long as the 
attenuation is small). 


At the surface of a conductor, the wave is almost 
completely reflected. The electric field penetrates 

a very small distance into the conductor, but the 
disturbance is evanescent, dying away exponentially 
within the skin depth, which is much smaller than the 
free space wavelength of the radiation. A small 
amount of energy is dissipated in the surface layer by 
the currents set up by the oscillating electric field. As 
with dielectrics, the electric and magnetic fields 

are transverse to the direction of propagation and 
orthogonal to each other, but they are out of phase by 
m/A4 for a good conductor. 


Solution 5.3 The electric field, and therefore the 
current density, decays exponentially with depth 
inside the conductor: 


J(x£) = cE(2) = oEp exp|—2x/6], 


where o is the conductivity and the skin depth 

is 0 = \/2/90w. Since a > t, we ignore the 

field penetration at the edges, and treat this as a 
one-dimensional, planar problem, with J(a) constant 
across an element of area of width a and thickness 6x 
at distance x below the surface. The total current 
flowing across a plane perpendicular to the length is 


f= 2 | o Ep exp|—2/6] ada, 
0 


where the factor of 2 takes account of the penetration 
of field from both surfaces of the strip, and the upper 
limit of the integral can be taken as infinity as long as 
the thickness ¢ is much greater than the skin depth so 
that the field is essentially zero in the middle. Thus 


I = |-26ac Eo exp[—2/64]] > = 2dac Ep. 
Now the voltage difference between the ends of the 
strip is given by V = EL, so its resistance is 
V EoL L 
R(f)=—= - ) 
I 2aokpdo(f) 2a0d(f) 


The DC resistance is given by 


i h L 
R(0) _ engt = 
area X 0 ato 
Ne) 
R(f) t 


R00) 26(f) 2V 2. 2 
The resistance is equivalent to that of a strip that is 


two skin depths thick, because of field penetration 
from the two surfaces. 


Solution 5.4 The cut-off frequency for the mth 
mode, which has m half-wavelengths of electric field 
variation fitted between the plates, is given by 


f 3.00 x 108 ms—! 
= = —a 6 aH 
— \ a 2x 0.10m 


The m = 6 mode has f, = 9.0 GHz, but the m = 7 
mode has f, = 10.5 GHz, so six TE modes can 
propagate at 10 GHz. 


Cc m 


= 1.5m GHz. 
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Solution 5.5 The magnetic field will have the same 
time dependence, exp[—iw1t], as the electric field, so 
we can write Faraday’s law as 

OB 
Ot 
The electric field has only an x-component, so its curl 
reduces to 


curl E = — = iwB. 


OF, OE, 
curl E = = ey By e, =iwB. 
Equating vector components on left and right, 
B, =0, 
1 OF 1 
v=o 8 * = — kocos 6 Ep 2ix 
iw OZ iw 
cos(kgz cos 6) exp|i(koy sin @ — wt)| 
2h 6 
=F a cos(kgz cos 8) x 
w 
exp[i(koy sin 6 — wt)], 
1 OE, 1 
B,=--; ee —— iko sin 6 Ep 2ix 
iw Oy iw 
sin(koz cos @) exp[i(koy sin @ — wt)| 
2ko sin 0 
=-ik)— = sin(koz cos @) x 
w 


exp|i(koy sin 0 — wt)]. 


Note that B has a component in the y-direction — the 
direction of propagation — so the magnetic field is not 
transverse. 


Solution 5.6 The dispersion relation is 
— w (my ("2)" ne 
= C2 a b , 
and this can be rewritten as 


one (E+ GY] 


Then 
dw 
Usroup = dkew 
mT \ 2 nT \ 2 =e 
aie PR lide + (=) +(F)| 
_ hey 
= 


This can be written in the alternative forms 


_¢ wr mut \ 2 nn la 
wow = 5 [ae — () - Ge) 
1/2 
1 mo : nr t 
=(¢ 
2a 2b 
Note that since Uphase = “/kgw, the phase and group 
speeds are related by VphaseVgroup = C’: 
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Solution 5.7 (a) The cut-off frequency for the 
TEmn mode is given by Equation 5.51 of Book 3: 


/ m\ 2 n\2 

jos (=) = (=) , 

where a and b are the longer and shorter dimensions 
of the waveguide. In this case, b = a/2 and so 


fran = 5, Vin? + An? 
a 
8-1 
_ 3.00 x 10°ms (8 a. + Ane 


2 x 0.034m 


= 4.411 Vm? + 4n2 GHz. 


The lowest four cut-off frequencies will be those for 
the TEi9, TE20, TEo1 and TE ; modes, for 

which the values of W/m? + 4n2 are 1, 2, 2 and 

V5, respectively, leading to cut-off frequencies of 
4.4 GHz, 8.8 GHz, 8.8 GHz and 9.9 GHz, respectively. 


(b) The phase speed for the TE, mode is given by 


WwW 
Uphase = 


wu /(S-(™)- (8) 


Substituting the values of the cut-off frequencies 
determined in part (a), the phase speeds for the four 
modes are calculated to be 


TEio TE TEo. TE 
Ushaes/ 10°mis > 3.3 6.4 6.4 18. 


mode 


The group speed is most easily determined using the 
result in Equation 5.45, 


a ne 
Ugroup = © i Usheaes 


which leads to 


mode TE;0 TE29 TEo1 TE 
Cou tO me 27 1.4 1.4 0.49. 
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Solution 6.1 Fora collisionless plasma, the relative 
permittivity function is given by 


where the plasma frequency is given by 


Ne€ 
Wn = 
¥ MEQ 
7 10m? x (1.60 *. 10-!'C)? 
~ \ 9.11 x 10-81 kg x 8.85 x 10-12 Fm! 
=H heb @ 106, 
The phase speed is 


c c CW 
UV a = — 
phase A Jew) jw _ a ’ 


and (as shown in Worked Example 6.1 of Book 3) the 


group speed is 
dw (dk\* 
~ \dw 


Ugroup = de 


-1 


O89 
d esl © 
dw Cc 
= 
(: 1 
Cc 2_ 42 
wu? = Ws 
Cc 
= —4/w?—w? 
Ww 
2 
Uphase 


Solution 6.2 The interstellar medium is a very dilute 
plasma, so it can be regarded as collisionless, and we 
will ignore any thermal effects, and assume that there 
is no significant magnetic field, so that the plasma can 
be regarded as isotropic. Then the relative permittivity 
function is given by 

2 
a mt 
where the plasma frequency in the interstellar medium 
is given by 


e(w) 


_ 3 x 104m-3 x (1.60 x 10-19C)? 
~ \ 9.11 x 10-3 kg x 8.85 x 10-12 Fm-! 
= 9.76 x 10°s-*. 


As shown in Worked Example 6.1 of Book 3, the 
group speed for a material with this relative 
permittivity function is 


Ugroup = C€ 


and since w >> wy for 100 MHz radiowaves and light, 
this can be approximated by 


1 Wp 
Gig =| ag— 14 
group 27 


The difference between the arrival times for the 
radiowave and red pulses is 


Thus for a frequency of 1.00 x 10° Hz, At= 2 2 
Ugroup radio Ugroup red 
— w2 — we D 1 Ws 1 Ws 
P — 1-3-5 1-3-5 
Cc c Wradio Wred 
4/1 — w2/w? D w w? 
° ef | | Lb peo” Le 
_ 3.00 x 108 ms! radio Wred 
4/1 — (5.64 x 106 s-1/2n x 1.00 x 10 s-1)? _ Dus 1 1 | 
= 6.8 x 108 ms~}, 2¢ |Wiadio Yreal 
aad Since Wradio < Wred, We can ignore the second term 
. in brackets, so 
Cc 
Cssin= =1,5 «10? ms. 1019 m x (9.76 x 10% s—1)? ia 
Os. 


Uphase 


The speeds for frequencies of 2.00 MHz and 100 MHz 
are calculated in a similar way, and the values are 
tabulated below. 


1.00MHz 2.00MHz 100MHz 
tjese/(l0°ms + 6.8 3.4 3.0 
tino] Ome 818 27 3.0 


At 


~ 2x 3.00 x 108ms—!(2m x 108Hz)2 


Solution 6.3 The relative permittivity of a 
collisionless isotropic plasma is given by 


we 


_ Pp 
=l|- 2 

When w < wp, this is negative, which means that the 
refractive index n (= ,/e) and the wavenumber 

k (= nw/c) are imaginary. The electric field is 


therefore evanescent within the plasma. Neglecting 


e(w) 
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any spreading of the beam, the amplitude of the 
electric field will fall off as 


Re aga 
E P 
exp|ikz] = exp ee 4 = exp | i—————__ z 
c c 
ie cre 
= exp | — Zz 
c 

wpy/1 — w*/w2 

= exp 2 
c 


The value of this expression just inside the window, at 
z = 0, is 1, so the ratio of the amplitude of the electric 
field in the plasma at a distance z from the window to 
the amplitude just inside the window is just 


Wpy/1 — w?/w2 


exp z 
Cc 


For z = 0.5m and w/w, = 0.99, this ratio is 


Qn x 108 s-! V1 — 0.992 
exp 0.5m 


3.00 x 108? ms~! 
= 0.86. 


For w/w, = 0.90 the ratio is 0.63, and for 
w/Wy = 0.50 it is 0.40. 


The limiting value of this ratio when w/w, < 1 is 
exp[—7/3] = 0.35. 


Solution 6.4 (a) Since wy = \/nee?/meéo, the 
electron number density is 
_ amen _ 
Ne = 5 = 
e 


(8.0% 10% s*)? x9.11% 10-9 ke 8.85 x 10-4 Fim 


(1.60 x 10-19)? 


=2:0% 10" m-*. 
Also, since w, = eB/m, the magnetic field strength is 
pa me 
e 
_ OL 10k 43.0% 10% 5-* 
= 1.60 * 10-"'C 
=e 10" 7: 


(b) The radiowave frequency is much less than the 
plasma frequency and much less than the cyclotron 
frequency. This means that the relative permittivities 
for the RH and LH circularly polarized modes in the 
anisotropic plasma are 


2 2 
Ww Ww 

gRH_ “Po ang -tHA_—P 
WWe WW 


as shown in Subsection 6.3.2 of Book 3. The LH 
mode is evanescent, so it is low frequency RH 
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circularly polarized waves that are produced in the 
plasma. Their dispersion relation follows from 


(Ph a FE (oe a 
Cc Cc Cc We 


Then 


—1 
RH dkRH 
Ueroup a dw 
1 


— {2 
De Malle, 


4 /wu, 
= 2c =" 
Wp 


-1 


Substituting the data provided, 
up 2X 3.00 * 10° ms ** 
V2n x 13.56 x 10° Hz x 3.0 x 108s! 


8.0 x 1010 s-1 
=1.2x10°ms"!. 


Solution 6.5 (a) The microwave frequency Wyicro 
must equal the electron cyclotron frequency uw, for 
resonance to occur. Thus 


eB 27¢ 
We = — = Wyicro = ; 
m Amicro 
and so 
27cem 
B= 
€Xmicro 
_ 2x x 3.00 x 10° ms! x 9.11 x 107*' kg 
7 1.60 x 10-19C x 3.0 x 10-2m 
= 0.36T. 


(b) For efficient heating we require that the collision 
frequency, 1/7., is much less than the cyclotron 
frequency, so that the microwave electric field can 
accelerate the electrons over a number of cycles 
before the electrons are scattered. From Figure 

6.7 of Book 3 the height in the atmosphere where 
the pressure is 10* Pa is about 18 km, and from 
Figure 6.9 the collision frequency at this height (and 
pressure) is about 10!1 s~!. This is greater than ux 
(~ 2 x 101° s~!), so the heating would not be very 
efficient. 
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Solution 7.1 (a) Increasing the fibril diameter will 
increase the scattering cross-section, since there is 
more material to polarize and radiate energy. 


(b) Increasing the refractive index of the fibril will 
increase the scattering cross-section, since this change 
increases the relative permittivity and therefore 
increases the polarization. 


(c) Increasing the refractive index of the matrix will 
reduce the cross-section, because it reduces the 
difference between the polarization of the fibril and 
the matrix. If the fibrils and matrix had the same 
refractive index there would be no scattering. 


(d) Increasing the spacing of the fibrils would not 
affect the cross-section of a single fibril. The fibrils 
are weak scatterers and are not affected by their 
neighbours. 


(e) Increasing the wavelength of the radiation would 
reduce the scattering cross-section. The scattering will 
increase with frequency, just as Rayleigh scattering 
from molecules increases with frequency. 


(f) Changing the polarization of the incident 
radiation from parallel to the fibril to perpendicular to 
the fibril will reduce the scattering cross-section. 
When the polarization is perpendicular to the fibril, 
there is no scattering in the polarization direction, and 
scattering from the different elements of the fibril will 
also cancel in the direction parallel to the fibril. This 
leaves scattered radiation concentrated around the 
forward and reverse directions. 


Solution 7.2 The fundamental resonance 
corresponds to 5A of the radiation in the length L of 
the hairpin, that is, \ = 40. There is a node of 

the electric field at the short-circuited end and an 
antinode at the open end. Other resonances occur, 
each with the same boundary conditions at the ends 
of the hairpin. These all have an odd number of 
quarter-wavelengths fitted into the length of the 
hairpin, so 3\/4 = L, 5\/4=L, 7A/4 = OL, 

éte,, Or A = 4L/m, where im = 3, 5, 7, ete. 

The frequencies of the resonances are given by 

f =c/A =mc/4L = mfo, so other resonance might 
be observed at 1.5 GHz, 2.5 GHz, 3.5 GHz, etc. 


Solution 7.3. (a) The fundamental resonance 
frequency f;es in a material in which the speed of 
electromagnetic radiation is v, the refractive index is n 
and the relative permittivity is € is given by 

v Cc Cc Cc fo 


res Xr nr afE a 4/JeL = Je 
So in water, with ¢ = 80, 
tres = 100 MHz/\/80 = 11 MHz. 


(b) In ethanol, with « = 25, 
fres = 100 MHz/v 25 = 20 MHz. 


(c) The relative permittivity of a cold collisionless 
plasma is given by 


2 2 
Pa jf — P: 
Eeft = 1 2 f2’ 
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where the plasma frequency is 


_ “p 

j= Oe 
_ A. ise? 
7 Qn MEQ 


LF 20% 10" m-* x (1.60% 10-™ Cc)? 
~ Qn \ 9.11 x 10-31 kg x 8.85 x 10-12 Fm-1 


= 40.11 MHz. 


Thus 
40.11? 
=l]-— 
1002 
= 0.8391, 


fres = 100 MHz/ 0.8391 = 110 MHz. 


Solution 7.4 The fundamental resonance frequency 
in vacuum is 

re c _ 3.00 x 108 ms! 

° AL 4x 1.00 x 10-2m 


In a plasma, the fundamental resonance frequency is 
given by 


— 7.5 x 10°Hz. 


fo fo 


? 
€ 
a ae or ee 


where the plasma frequency is 


fo = (1/277) \/ nee? /meo. This simplifies to 
resp = fp + So. 


In this problem fies, p = 2 fo, so i. =3 ta Hence 


Freep = 


which can be rearranged to give 


= 127?me0 fe 


Ne 5) 


€ 
ea SOA ke 
(1.60 x 10-19)? 
x 8.85 x 107'* Fm™! x (7.5 x 10° Hz)” 


=2.1 x 10'8 m=. 
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